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1.0  INTRODUCTION 


Pressure  vessels  are  used  extensively  in  both  ground  and  spacecraft  applications. 
Explosive  failures  of  vessels  are  rare  due  to  precautions  normally  taken  including  adherence  to 
consensus  design,  fabrication  and  test  codes  and  standards.  Inservice  integrity  is  maintained 
through  monitoring  of  vessel  service  conditions  and  cyclic  history.  Yet  pressure  vessels  do 
occasionally  fail,  releasing  significant  energy  and  possible  hazardous  commodities  into  the 
surroundings.  Often  it  is  prudent  to  assess  the  damage  that  could  result  from  explosive  failure 
when  locating  pressure  vessels,  designing  nearby  structures  and  equipment,  performing  pressure 
tests,  or  considering  other  safety  precautions. 

A  considerable  body  of  data  exist  on  damage  and  injury  due  to  blast  wave  and 
fragmentation,  much  of  it  from  research  using  TNT  or  similar  high  explosives.  However 
substantially  less  is  known  about  blast  and  fragmentation  of  bursting  pressure  vessels  than  of 
chemical  explosions  such  as  TNT.  Further,  current  methods  documented  in  standards,  handbooks 
and  other  references  used  to  quantify  expected  energy  release,  blast  waves,  and  fragmentation  are 
inconsistent  and  vary  in  results*.  Accordingly,  a  pressure  vessel  burst  test  program  was  conducted 
for  the  USAF  -45th  Space  Wing  and  NASA  Headquarters.  The  program  studied  the  blast  wave 
and  fragmentation  of  bursting  gas  filled  pressure  vessels. 

1.1  Background 

An  explosive  rupture  of  a  pressure  vessel,  where  the  stored  energy  is  released 
instantaneously,  creates  a  blast  wave  (i.e.,  shockwave)  in  the  surromiding  air  and  propels 
fragments.  The  shockwave  and  fragment  characteristics  depend  on  such  things  as  vessel  contents, 
pressure,  vessel  geometry  and  mode  of  vessel  failure. 

As  the  blast  wave  advances,  the  energy  is  spread  over  the  wave’s  frontal  area,  which 
increases  with  the  square  of  the  distance  from  the  point  of  rupture.  Overpressure,  blast  wave 
velocity  and  therefore  blast  effect,  decrease  rapidly  with  distance.  After  passage  of  the 
shockwave,  the  pressure  decreases  until  a  suction  phase  follows  in  which  pressure  drops  below 
normal  atmospheric  pressure.  The  negative  pressure  is  a  result  of  the  outrush  of  gases  from  the 
center  of  the  rupture  causing  an  overexpansion.  The  pressure  above  atmospheric  at  the 
shockwave  front  is  the  peak  overpressure  and  is  used  with  impulse  to  establish  the  relative  hazard 
(i.e.,  shockwave  intensity  and  energy  in  the  shockwave,  the  impulse  being  the  area  under  the 
positive  portion  of  the  pressure  versus  time  curve)  associated  with  ruptures  and  explosions  at  a 
given  distance.  The  blast  wave  emanating  from  a  bursting  pressure  vessel  is  somewhat  similar 
to  that  caused  by  a  high  explosive  detonation.  The  pressure  close  in  (0  to  10ft  in  vessels  tested) 
due  to  vessel  burst  is  generally  lower  than  high  explosive  detonation  and  is  a  function  of  burst 
pressure.  This  is  because  the  pressure  at  the  vessel  surface  is  less  than  that  of  a  high  explosive 
blast  at  the  same  distance  from  the  explosion  center.  Other  variations  are  caused  by  vessel  and 
failure  geometry  and  distance  from  a  firm  reflecting  surface.  Blast  wave  overpressures  for  vessel 
failures  with  very  large  ratios  of  vent  area  to  vessel  volume  may  be  estimated  using  methods  of 
Baker^  or  Held  and  Jager^. 
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The  explosive  failure  of  a  pressure  vessel  not  only  generates  a  blast  wave  but  produces 
fragments,  with  very  high  velocities  possible.  Fragments  constitute  a  significant  hazard  to 
personnel,  systems,  components  and  structures  in  the  vicinity.  Primary  fragments  are  portions 
of  the  vessel  or  its  attachments  that  are  accelerated  due  to  the  internal  pressure  of  the  vessel. 
Secondary  fragments  may  also  be  produced  due  to  the  action  of  the  blast  wave  or  primary 
fragments  on  nearby  objects. 

Studies^’^’"'  of  the  characteristics  of  vessel  fragments  have  addressed  the  velocities  of 
fragments  produced,  their  trajectories  and,  as  a  result,  their  ranges  and  their  impact  velocities. 
Determination  of  the  initial  velocities  of  fragments  has  been  undertaken  by  several  researchers. 
Most  such  studies  are  based  upon  work  by  Taylor  and  Price^  which  predicted  the  velocities  of 
two  spherical  vessel  fragments  accelerated  by  an  expanding  adiabatic  ideal  gas.  Wiedermann®  has 
shown  that  real  gas  effects  can  be  expected  to  reduce  the  fragment  velocity  from  the  ideal  gas 
prediction. 

Once  the  initial  velocity  of  a  fragment  has  been  determined,  its  range  may  be  found 
through  ballistic  calculations,  generally  done  through  the  use  of  a  computer  code,  a  number  of 
which  are  available.  Code  considerations  are  drag  coefficient,  lift  coefficient  (if  any),  initial 
trajectory  angle  and  reference  area  -  either  fixed  or  varying  (tumbling  or  gradually  changing). 

The  reader  is  also  referred  to  work  by  Pittman’  for  velocities  of  fragments  burst  from 
flight  weight  vessels  at  600  and  8000  psig.  Baum'*  *  compiles  data  from  his  own  work  and  other 
researchers  on  fragments  from  vessels  burst  at  70  to  4400  psig,  mostly  at  750  psig  or  less. 

1.2  Test  Program 

A  test  program  matrix  was  developed  that  included  a  series  of  test  plans  each  with 
multiple  pneiunatic  vessel  bursts.  The  objective  of  the  program  matrix  was  to  force  vessel  bursts 
in  such  a  way  as  to  generate  worst  Ccise  blast  waves  and  fragmentation,  such  that  a  model  could 
be  developed  that  would  envelop  generally  expected  vessel  failures.  The  latter  test  plans  of  the 
matrix  would  include  such  representative  vessel  failures.  Worst  case  however  is  a  function  of 
several  variables,  including  location  and  orientation  of  failure,  pressure,  vessel  shape,  fragment 
type  and  number,  and  height  above  ground.  The  plans  and  tests  comprising  the  program  matrix 
have  been  developed  to  minimize  the  number  of  vessel  bursts  yet  meet  the  stated  objective  with 
valid  data. 

In  the  development  of  a  test  matrix,  it  was  also  recognized  that  a  pressure  vessel  burst 
may  not  produce  a  spherical  shockwave  as  does  a  TNT  explosion.  The  blast  wave  from  a 
pressure  vessel  burst  may  be  much  stronger  in  one  direction  than  another  based  on  how  the  vessel 
shell  comes  apart.  To  provide  a  direct  experimental  comparison  with  pressure  vessel  bursts, 
approximately  spherical  high  explosive  detonations  have  been  conducted  as  part  of  the  test 
program. 
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Accordingly,  a  test  program  matrix  was  developed  which  incorporated  varied  failure 
locations  and  mechanisms.  Seven  test  plans  were  envisioned  with  each  test  plan  consisting  of 
several  vessel  bursts.  The  final  and  truncated  test  matrix  includes  5  pneumatic  burst  series  and 
is  provided  in  the  Test  Summary  section,  see  Figure  2.1.  Parameters  varied  in  the  matrix  are 
burst  pressure,  height  of  burst  (HOB),  split  location,  number  of  fragments  and  vessel  geometry. 

2.0  TEST  SUMMARY 

This  summary  is  presented  early  in  this  report  because  it  will  be  referred  to  by  other 
sections.  The  section  summarizes  the  tests  and  associated  data  which  describes  the  vessels  and 
the  test  conduction.  Primary  test  output  data  are  deferred  to  other  sections,  i.e.;  Section  5  for 
blast  overpressxires  and  Section  6  for  fragment  velocities. 

2.1  High  Explosive  Tests 

The  high  explosive  tests  initially  (PHE-1  and  PHE-2)  had  a  twofold  purpose:  1.)  as  a 
checkout  device  to  ensure  all  pressure  data  chaimels  were  operating  and  2.)  as  a  calibrating  tool 
for  comparing  high  explosive  overpressure  to  pressure  vessel  burst  overpressure  when  both  vessel 
and  high  explosive  had  approximately  the  same  TNT  equivalence.  (See  Section  4.1  for  vessel 
TNT  equivalence.)  All  the  high  explosive  shots  are  summarized  in  Table  2.1.  Shot  PHE-2  was 
detonated  first  because  pressure  vessel  P-2  was  to  be  burst  first.  The  charge  weight  was  intended 
to  be  approximately  60  lbs  TNT  equivalent  to  match  the  energy  of  vessel  P-2  at  its  intended  burst 
pressure  of  6000  psig.  However,  65  pounds  of  Composition  C-4  was  furnished  (instead  of  50) 
which  is  very  energetic  and  has  a  TNT  equivalence  of  85.5  lbs.  Similarly  the  charge  weight  for 
shot  PHE-1  is  an  approximation  to  the  TNT  equivalence  of  vessel  P-1.  All  the  high  explosive 
TNT  equivalence  conversions  were  made  using  the  equivalent  weight  in  relation  to  TNT  for  the 
average  pressure  at  the  mid-arena  distance  of  30'  using  Swisdak’s^  report. 

For  Test  Plan  #1,  high  explosive  shots  HEl-1  through  HEl-3,  charge  weights  of  pentolite 
were  used  which  would  provide  an  overpressure  variation  but  would  not  exceed  the  calibrated 
range  of  the  pressure  transducers  at  the  ten  foot  distance.  The  use  of  high  explosive  shots  as  a 
calibration  check  tool  had  been  abandoned  because  of  the  greater  output  at  close  distances  for 
high  explosives.  Resetting  gains  defeated  the  pxirpose  and  using  calibrations  for  high  explosives 
provided  resolution  problems  for  vessel  burst.  The  TNT  equivalence  of  shot  HE  1-2  was  expected 
to  be  10  lbs  instead  of  14.6  and  caused  all  outputs  at  10  ft  distance  to  be  off  scale. 

The  geometry  of  the  high  explosive  shots  is  as  follows:  pentolite  -  all  1:1  diameter  to 
length  cylinders,  composition  B  - 1 : 1  diameter  to  length  cylinder,  C-4-close  to  1 : 1 : 1  square  (small 
blocks  stacked  together). 
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Table  2.1 

Summary  of  High  Explosive  Shots 


Shot# 

Date 

NSWC 

Shot# 

lbs.  Pent 
(other) 

Of&et, 

ft 

TNT  Equiv., 
lbs. 

HOB 

ft 

PHE-2 

6-21-90 

3201 

65,  (C-4) 

0.0 

85.5 

3.5 

PHE-1 

6-26-90 

3249 

33.3,  (Comp  B) 

0.0 

40.0 

3.5 

HEl-1 

1-28-91 

3250 

.66 

0.0 

.89 

3.5 

HEl-2 

1-28-91 

3256 

9.9 

0.0 

14.6 

3.5 

HE2-3 

7-16-91 

3399 

4.5 

1.0 

6.6 

14.0 

HE2-1 

7-18-91 

3401 

4.5 

0.0 

6.6 

3.5 

HEl-3 

7-20-91 

3526 

4.5 

0.0 

6.6 

3.5 

HE2-2 

7-20-91 

3404 

4.5 

1.0 

6.6 

8.7 

For  TP  #2  (HE2-1  through  HE2-3)  all  the  high  explosive  shots  had  the  same  equivalent 
weight  and  provided  an  opportunity  to  compare  reflection  factors  for  high  explosives  with  that 
for  center  split  pressure  vessel  bursts  at  a  single  pressure  level. 

Initially  pressure  vessels  were  located  in  the  arena  with  the  edge  of  the  vessel  above 
ground  zero,  creating  a  one  foot  offset  for  the  two  foot  diameter  vessels  (see  "offset"  in  Table 
2.3).  The  elevated  height  of  burst  shots,  HE2-2  and  HE2-3,  were  rigged  with  wooden  supports 
on  the  vessel  support  stand  and  were  thus  at  the  same  one  foot  offset.  Other  shots  followed 
NSWC  usual  procedures. 

After  July  1991  no  high  explosive  shots  were  detonated.  The  high  explosive  shots  needed 
to  be  detonated  before  the  pressure  vessel  bursts  because  of  ground  zero  soil  erosion  which 
frequently  accompamed  the  vessel  bursts.  Detonating  high  explosive  shots  first  prevented 
simultaneous  preparations  of  the  arena  field/transducer  recording  setup  by  instrumentation  and 
vessel  setup  by  riggers,  explosive  techmcians  and  General  Physics  Corporation  personnel.  Since 
a  quantity  of  high  explosive  data  had  been  recorded  using  the  vessel  instrument  setup  and  since 
the  data  were  of  lesser  importance,  its  collection  was  abandoned. 

2.2  Pneumatic  Vessel  Bursts  and  LSC  Detonations 

Test  events  included  LSC  detonations  and  vessel  bursts,  sometimes  separately  and 
sometimes  as  a  combined  event.  The  detonation  of  17  linear  shaped  charges  during  preliminary 
testing  are  listed  in  Table  2.2.  Field  overpressure  data  without  a  vessel  burst  were  recorded  for 
the  first  13  detonations,  through  2-2.  Number  2-3  (vessel  burst  P-1)  resulted  in  the  recording  of 
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combined  LSC  and  vessel  burst  overpressure.  No  data  were  recorded  for  the  latter  three  LSC 
detonations,  3-1  through  3-3.  All  the  LSCs  were  wrapped  around  a  two-foot  diameter  pressure 
vessel  with  a  six  inch  double  pigtail  extending  to  the  ungaged  side  of  the  arena.  These  vessels 
were  also  offset  one  foot  toward  the  ungaged  side  of  the  arena.  Events  are  discussed  under 
Section  3.1.2,  Preliminary  Pneumatic  Test. 

Figure  2.1  is  a  test  matrix  of  all  the  bursts,  hydraulic  and  pneumatic. 

Table  2.2 

Shaped  Charges  Detonated  During  Preliminary  Testing 


Prelim. 

Test 

Segment 

LSC 

Charge 

Number 

Vessel 

Burst 

Number 

.  Vessel 
Pressure, 
psig 

LSC  charge 
wt,  gr/ft. 
(lead/RDX) 

Overpressure 

i^ta 

1 

lA-1 

— 

0 

25 

.03 

yes 

1 

lA-2 

— 

0 

10 

.01 

yes 

1 

lA-3 

— 

0 

10 

:01 

yes 

1 

lA-4 

— 

0 

25 

.03 

yes 

1 

lA-5 

— 

0 

15 

.02 

yes 

1 

lA-6 

— 

0 

15 

.02 

yes 

1 

lA-7 

— 

0 

25 

.03 

yes 

1 

lB-1 

— 

3200' 

200 

.27 

yes 

1 

lB-2 

— 

3200 

50 

.07 

yes 

1 

lB-3 

— 

0 

50 

.07 

yes 

1 

lB-4 

— 

3200 

200 

.27 

yes 

2 

2-1 

— 

2850 

25 

.03 

yes 

2 

2-2 

— 

2950 

25 

.03 

yes 

2 

2-3 

P-1 

3250 

50 

.07 

with  vessel 
burst 

3 

3-1 

— 

0 

50 

.07 

no 

3 

3-2 

— 

0 

50 

.07 

no 

3 

3-3 

— 

2700^ 

200 

.27 

no 

‘Planned  p 
^Planned  p 

iressure  =  6000 
»ressure  =  6000,  2700  ,oc 

curred  on  depressuring  from  5800 

5 


g<N<N“* 
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FIGURE  2.1 

TEST  MATRIX  SHOWING  ACTUAL  BURST 


Table  2.3  lists  all  the  pneumatic  vessel  bursts  and  the  LSC  detonations  which  occurred 
as  a  combined  event.  The  TNT  equivalence  of  both  the  vessel  and  the  LSC  are  provided.  For 
vessel  TNT  equivalence  see  Section  4.1.  LSC  TNT  equivalence  was  obtained  using  Table  I  of 
Kinney  and  Graham’®. 

Table  2.3 

Summary  of  Vessel  Burst  Events 


Date 

NSWC 

Shot# 

Vessel  # 
(WB#) 

Vess 

ft-es 

psig 

1-25-92 

3573 

PC 

3975 

6-26-90 

3248 

P-1 

3250 

6-25-90 

3245 

P-2 

4700 

1-25-91 

3251 

1-1 

1475 

1-29-91 

3252 

1-2 

3450 

1-26-91 

3253 

1-3 

5425 

7-19-91 

3402 

1-4 

7125 

7-20-91 

3403 

2-1 

3450 

7-17-91 

3400 

2-2 

3450 

7-15-91 

3398 

2-3 

3475 

6-16-92 

3574 

5-1  (5) 

3600 

6-17-92 

3575 

5-2  (6) 

3600 

6-19-92 

3577 

5-3  (7) 

3600 

6-18-92 

3576 

5-4(8) 

3600 

11-2-93 

4053 

6A-1  (13) 

3280 

11-3-93 

4054 

6A-2 

4000 

11-3-93 

4055 

6A-3  (4) 

3300 

11-4-93 

4056 

6A-4  (10) 

3500 

Lead/RDX 


Lead/RDX 


Lead/RDX 


Lead/RDX 


Lead/RDX 


Lead/RDX 


Lead/RDX 


Lead/RDX 


Lead/RDX 


Lead/RDX 


AL/RDX 


AL/RDX 


AL/RDX 


AL/RDX 


AL/RDX 


AL/RDX 


AL/RDX 


AL/RDX 


TNT  Equiv.,  Offset,  HOB 

lbs.  ft  ft 


Grain 

/foot 

Vess 

LSC 

100 

2.0 

.125 

50 

31.7 

.07 

0 

47.0 

0.0 

200 

13.5 

.27 

50 

33.8 

.07 

25 

54.8 

.03 

100 

73.4 

.13 

25 

33.8 

.03 

25 

33.8 

.03 

25 

34.1 

.03 

50 

35.3 

.07 

50 

35.3 

.07 

50 

35.3 

.07 

50 

35.3 

.07 

400 

32.2 

3.41 

200 

75 

2.0 

.07 

200 

32.3 

.29 

100 

34.3 

2.18 

200 

nominal  (mfg.  guarantee  volume)  except  PC,  6A-2  (measured) 
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Table  2.4  contains  data  needed  to  compute  fragment  velocities:  pressure,  pressurant  and 
ambient  temperatures,  vessel  dimensions  and  fragment  weights. 


Table  2.4 

Data  for  Computing  Fragment  Velocities 


Vessel  #/ 

Vessel 

Temperature,  ‘T 

Inner 

Radius,'* 

Outer 

Radius 

Length  of 
ft 

(cyL  sect.) 

Vessel  or 
Fragment 
WP”,  lbs 

Fragment 

Pressure 

psig 

Vessel 

gas 

Ambient 

Air 

fl 

^  ft 

P-2 

4700 

124 

85 

.93’ 

1.0 

19’ 

5775 

P-1 

3250 

NR‘ 

NR 

.93’ 

1.0 

19’ 

5800 

1-1 

1475 

65 

31 

.929 

1.0 

19’ 

5525 

1-3 

5425 

75 

oo 

.93’ 

1.0 

19’ 

5825 

1-2 

3450 

109 

55 

.93’ 

1.0 

19’ 

5900 

2-3’ 

3475 

142 

85 

.93* 

1.0 

19* 

5400 

2-2’ 

3450 

132 

90 

.924 

1.0 

19* 

5300 

M 

7125 

133’ 

93’ 

.93’ 

1.0 

19’ 

5250 

2-1 

3450 

135 

88 

.93’ 

1.0 

19’ 

5025 

PC' 

3975’ 

125 

42 

.86 

.88 

0/0 

43.6 

5-1/Lg 

3600 

130 

80 

.93’ 

1.0 

14.0 

3500 

5-1/sm 

1.0 

5.0 

1425 

5-2/lg 

3600 

140 

84 

.93’ 

1.0 

14.2 

3600 

5-2/sni 

1.0 

4.6 

1500 

5^/lg 

3600 

123 

78 

.93’ 

1.0 

19.2 

4675 

5-4/sm 

1.0 

.1 

300 

5-3/Ig 

3600 

126 

80 

.93’ 

1.0 

18.9 

4675 

5-3/sm 

1.0 

.1 

300 

6A-1’ 

3280 

106 

50 

1352 

1.5 

7.0’ 

2x800 

6A-2= 

4000^ 

169 

44 

.86 

.88 

0/0 

43.6 

6A-3 

3300 

97 

47 

1348 

1.5 

7.0’ 

6100 

6A-4  sidewall 

3500 

108 

60 

.929 

1.0 

18.7 

12x362 

6A-4  endcap 

2x352 

*  Preliminaiy  composite  overwrapped  pressure  vessel  (COPY) 

^  Vessel  volume  =  2.7  cubic  feet,  all  others  53 

'  TP  #6  Deleted 

*  At  groove 

^  Data  not  recorded.  This  is  from  previous  day’s 
pressurization  to  7200  psig  (similar  pumpup  time  and 
ambient  temperature). 

®  Not  recorded 

^  of  cylindrical  sections,  of  vessel  6A-4  length  is  point  of 
curvature  (PC)  to  PC 
*  Nominal 

’  Height  of  burst  of  V2-2  and  V2-3  is  8.7'  and  14.0'.  All  others  is 

3.5' 

See  Table  6-2  for  individual  fragment  weights 
"  Weighed,  not  calculated 
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Of  the  18  vessels  burst  using  pneumatic  pressure  (gaseous  nitrogen)  16  were  cylindrical 
steel  vessels  and  two  were  composite  spheres.  Further  information  on  the  vessels  is  provided  in 
Figure  2.1.  All  vessels  were  53  cubic  feet  voliune  except  for  the  2.7  cubic  feet  spheres.  The 
vessel  materials  and  pressure  ratings  are  listed  in  Table  2.5. 


Table  2.5 

Test  Vessel  Ratings 


Size 

Ratings,  Materials 

24"  cylinders: 

SA-372,  2450  psi  (ASME  Section  VIII,  Div  1,  App.  22)  (14  were  burst) 

36"  cylinders: 

SA-516,  1770  psi  (ASME  Section  VIII,  Div.  1)  (2  were  burst) 

spheres: 

cryostretched  301  stainless  steel  liner  with  Kevlar-epoxy  overwrap,  4000 
psi  (MIL-STD-1522A) 

Since  the  ground  surface  conditions  in  the  center  of  the  arena  may  have  affected  blast 
ou^ut,  these  conditions  are  tabulated  in  Table  2.6. 


Ground 


Table  2.6 
Zero  Surface 


Conditions 


3.0  TEST  TECHNIQUES 


Test  techniques  include  all  the  activities  surrounding  the  tests,  including  burst  initiation, 
vessels,  pumping  and  other  hardware,  instrumentation  and  observations. 

3.1  Establishment  of  Techniques 

Two  preliminary  burst  test  series  were  used  to  establish  basic  techniques.  However 
techniques  continued  to  evolve  throughout  the  program.  The  first  burst  test  series  used  water  as 
a  pressurant  and  was  conducted  at  Cape  Canaveral  Air  Station,  FL.  For  the  second  test  series 
a  gaseous  nitrogen  pressurization  system  was  assembled  at  the  Naval  Surface  Warfare  Center  at 
Dahlgren,  VA  and  testing  was  conducted  at  that  facility. 

3.1.1  Hydroburst  Test 


Planned  burst  initiation  for  pneumatic 
burst  tests  consisted  of  detonating  a  shaped 
charge  in  a  premachined  groove  with  the 
vessel  at  the  desired  burst  pressure.  This 
permits  selection  of  the  failure  geometry, 
burst  pressure  and  time  of  burst.  The  latter  is 
important  because  of  the  limited  time 
available  with  high  speed  video  tapes  (30 
minutes),  instrumentation  tapes  (30  to  60 
minutes)  and  high  speed  photography  (as  little 
as  9  seconds). 

The  object  of  the  test  was  to  check  out  failure  stress  levels  for  burst  initiation  piu^oses. 
The  use  of  a  linear  shaped  charge  (LSC)  detonated  in  a  premachined  groove  was  considered 
essential  for  utilizing  a  small  shaped  charge.  A  small  LSC  was  necessary  to  minirmVft  the 
undesirable  explosive  blast  wave  overpressure.  Undesirable  because  the  overpressure  would  later 
be  measured  and  if  possible  attributed  to  the  vessel  burst.  Figure  3.1  is  a  drawing  of  the  vessel 
groove  geometry. 

The  desired  stress  level  was  to  be  less  than  yield  by  some  safety  margin  on  pressurization 
and  achieve  ultimate  on  shaped  charge  detonation.  It  was  felt  that  with  a  narrow  circumferential 
groove  the  longitudinal  stress  alone  would  determine  failure.  It  was  expected  that  local  yielding 
would  maintain  the  tangential  stress  to  a  value  near  the  yield  strength  at  a  pressure  condition  and 
groove  geometry  which  would  give  a  calculated  tangential  stress  at  the  groove  outer  diameter 
above  ultimate.  Additionally,  the  longitudinally  short  groove  section  was  expected  to  be 
reinforced  by  the  continuous  moment  into  the  adjacent  heavy  section  thus  preventing  large 
deflections  with  resultant  high  stresses.  A  preliminary  finite  element  stress  analysis  by  a 
consultant"  could  not  confirm  these  expectations.  It  was  suggested"  that  tradeoffs  in  cost  versus 
accurate,  positive  results  were  in  favor  of  testing  rather  than  more  detailed  modeling  and  that 
failure  could  be  expected  beyond  a  maximum  effective  stress  of  (Sy+SJ/l. 


Figure  3.1,  Typical  Vessel  Cross  Section 
Showing  Groove  and  Shaped  Charge 
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A  two  vessel  test  was  planned,  using  groove  depths  which  would  be  useful  for  a  future 
preliminary  pneumatic  burst  test.  The  vessel  wall  thickness  was  drawing  listed  as  0.817" 
minimum  with  no  roundness  or  wall  thickness  tolerance  provided.  Measurements  showed  vessel 
HI  to  be  0.190"  (indicator  reading)  out  of  round  and  vessel  H2  to  be  0.245"  out  of  rotmd.  The 
wall  thickness  (ultrasonic  thickness  gage)  varied  from  .874"  to  .926"  on  vessel  1  and  from  .847" 
to  .929"  on  vessel  2.  A  search  was  initiated  for  a  suitable  computer  numerically  controlled 
(CNC)  milling  machine  because  of  the  dimensional  variations.  However  none  was  found  and 
jury-rig  procedures  were  used  to  permit  using  a  lathe  to  machine  the  vessels  in  Central  Florida. 
General  Physics  Corporation  supported  one  of  the  procedures  with  software’^  to  facilitate 
offsetting  the  vessel  and  machining  eight  slightly  larger  radii  at  the  groove  location.  The  as- 
machined  grooves  are  shown  in  Table  3.1. 


Table  3.1 

Hydroburst  Vessel  Dimensions  at  Circumferential  Groove 


Vessel 

Average  O.D., 
in. 

Original  Wall, 
in.  (avg.  of  16) 

iiiiiiiiiiiii 

width,  in 

Groove  wall,  in. 
(avg.  of  26-31) 

HI 

24.012 

.928 

.420 

.392 

H2 

24.192 

.914 

.200 

.320 

Each  vessel  in  turn  was  filled  with  non-deaerated  water  and  air  was  bled  using  a  valve 
on  the  upper  surface.  The  maximum  pressure  was  held  to  7500  psig  due  to  limitations  of  the 
1/2"  stainless  steel  supply  tubing.  Vessel  HI  was  held  at  the  maximum  pressure  for  five  minutes 
and  depressurized  when  it  did  not  fail.  Vessel  H2  failed  at  6500  psig.  Pressurization  schedule 
after  6000  psig  was  achieved  was  to  increase  pressure  in  100  psi  increments  and  hold  briefly. 
A  fairly  loud  report  could  be  heard  950  feet  away  when  the  vessel  burst. 

Upon  vessel  H2  burst  a  high  response  transducer,  Endevco  model  8511  A,  recorded  a 
pressure  drop  to  ambient  in  2  milliseconds  (ms).  A  48  frame/second  motion  picture  camera 
recorded  the  event  and  film  records  later  showed  that  each  vessel  half  was  accelerated  to  73  feet 
per  second.  One  frame  from  the  camera  is  presented  as  Figure  3.5  in  Photo  Section,  Section  3.7. 
Accelerometers  produced  no  useful  data  due  to  a  failure  of  cement  bond  mounting.  Each  half 
vessel  stopped  in  about  57  feet  with  its  nose  partly  buried  in  sand. 

Longitudinal  stress  analyses  results  for  both  vessels  is  provided  in  Table  3.2.  Lab  tests'"* 
on  a  section  cut  from  one  of  the  fragments  provided  material  strengths  as  also  shown  in  Table 
3.2.  The  calculated  average  stress  in  vessel  H2  was  less  than  ultimate  at  failure.  Failure  is 
deemed  to  have  occurred  by  the  groove  section  stretching  to  failure  as  the  value  of  (Sy+Su)/2  = 
114  ksi  was  approached.  A  later  test,  (vessel  5-1)  Section  5.3,  with  a  13  ms  delay  to  failure 
confirms  this  belief  Note  that  tangential  stresses  in  a  cylinder  are  twice  longitudinal  and  that 
these  high  stresses,  well  beyond  failure,  never  occurred. 
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Table  3.2 

Calculated  Stresses  and  Measured  Strengths 


Vessel  Stress 
(ksi) 

Measured  Strength 
(ksi) 

H2 

Sy 

Su 

Longitudinal 

104.1 

112.0 

98 

130.3 

Tangential 

— 

— 

97 

113.7 

The  knowledge  gained  from  these  tests  confirmed  initial  expectations  and  permitted 
simplifying  the  following  preliminary  pneumatic  burst  test. 

3.1.2  Preliminary  Pneumatic  Test 

A  preliminary  pneumatic  test  was  conducted  at  the  Naval  Surface  Warfare  Center  at 
Dahlgren,  VA.  The  test  site  is  described  in  Section  3.5. 

3.1.2.1  Purpose 

The  purpose  of  the  preliminary  test  series  was  fourfold:  1)  measure  the  cutting  action  of 
various  linear  shaped  charges  (LSC)  proposed  for  use  on  test  vessels,  2)  collect  data  on  the 
overpressure  in  air  due  to  the  detonation  of  LSCs  of  various  explosive  weights  (so  that  LSC 
overpressure  effects  can  be  deleted  from  real  vessel  overpressure  data  as  one  would  delete  a  tare), 
3)  explore  the  combined  LSC/intemal  pressure  interaction  in  causing  vessel  failure  and  4) 
checkout  a  gaseous  nitrogen  (GN2)  pressurization  system  and  obtain  experience  with  the  joint 
NSWC/explosive  and  General  Physics  Corporation/pressure  testing  environment.  Some 
background  on  purpose  No.  3  is  appropriate.  If  the  LSC  cut  penetration  is  about  60%  or  more 
of  the  total  section  thickness,  the  effect  of  shock  loading  will  cause  the  section  to  break  at  the 
cut  line.  Then  for  a  section  initially  under  stress,  a  larger  section  thickness  might  break  tVian  the 
cut  thickness.  This  was  to  be  addressed  in  the  preliminary  test  by  providing  different  stress  levels 
at  3  grooves  on  each  of  two  pressurized  vessels  and  detonating  shaped  charges.  Vessels  are 
described  in  Section  3.3. 

3.1.2.2  Description  of  Test  Plan 

The  preliminary  test  was  performed  in  4  segments  (lA,  IB,  2  &  3)  as  shown  in  Table  3.3. 
With  the  first  vessel,  seven  shaped  charges  around  the  circumference  of  the  unpressurized  vessel 
were  to  be  detonated  sequentially.  An  array  of  pressure  transducers  would  record  pressure  as  a 
function  of  time  at  different  locations  in  the  adjacent  arena.  Next,  four  larger  shaped  charges 
were  to  be  detonated  sequentially  while  the  vessel  was  pressurized  to  a  planned  6000  psig.  The 
vessel  was  not  expected  to  burst  following  any  of  these  shaped  charge  detonations  as  the  shaped 
charge  cut  would  not  remove  enough  wall  metal  to  raise  the  axial  stress  above  the  steel’s  ultimate 
tensile  strength. 
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Table  3.3 

Preliminary  Pneumatic  Test  Summary 


Test 

Segment 

Vessel 

Number 

Pressure 

(psig) 

Number  of 
Shaped 
Charges 

Purpose 

lA 

— 

0 

7 

Measure  shaped  charge  overpressures 

IB 

— 

6000 

4 

Measure  shaped  charge  overpressures, 
avoid  surprise  failure  in  later  vessels 

2 

P-1 

3000 

3 

Study  shaped  charge/intemal  pressure 
interaction  in  causing  vessel  failure 

3 

P-2 

6000 

3 

Study  shaped  charge/intemal  pressure 
interaction  in  causing  vessel  failure 

The  third  segment  of  the  test  was  to  utilize  a  similar  pressure  vessel  which  had  three 
grooves  machined  around  the  circumference.  Shaped  charges  were  to  be  placed  inside  the 
grooves  and  the  vessel  would  be  pressurized  to  approximately  3000  psig.  The  shaped  charges 
were  to  be  detonated  sequentially  at  either  a  deeper  groove  or  with  a  larger  core  load  such  that 
the  axial  stress  developed  at  the  groove  would  increase  with  each  detonation. 

The  fourth  segment  of  the  test  was  to  utilize  a  similar  vessel  with  three  different  grooves 
and  pressurized  to  a  planned  6000  psig.  Larger  core  loads  of  shaped  charges  were  to  be  used  but 
the  test  would  otherwise  proceed  in  the  same  manner  as  the  second  test. 

Both  the  third  and  fourth  segments  of  the  test  were  designed  such  that  the  axial  stress  was 
high  enough  to  cause  failure  when  the  third  charge  detonated  but  failure  was  not  assured  when 
the  first  or  second  charge  detonated. 

3.1.2.3  Summary  of  Results 

Table  2.2  in  the  summary  section.  Section  2.2,  tabulates  all  the  shaped  charges  and 
accompanying  vessel  burst  for  this  test  series.  One  vessel  burst  (pressurized  imtil  burst),  vessel 
P2,  occmred  the  day  following  the  initiating  LSC  (No.  3-3)  failure  to  cause  a  burst  due  to 
premature  depressurization.  Measured  shape  charge  overpressures  are  discussed  in  Section  5.3. 
Vessel  burst  overpressures  are  included  in  Sections  5.4  through  5.15.  Fragment  flight  and 
velocities  are  included  in  Section  6. 

The  cutting  action  of  various  LSCs  on  SA-372  steel  was  measured  and  is  shown  in  Figure 
3.2.  Data  were  recorded  with  and  without  standoffs.  (Standoffs  used  were  rated  as  optimum  for 
C1018  steel.)  The  figure  shows  only  the  data  recorded  with  standoffs  which  was  typically  a  25% 
to  44%  deeper  cut  than  without  standoffs.  These  data  were  not  available  from  the  shaped  charge 
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manufacturer.  It  was  also  found  that  there 
was  no  shock  interaction  to  assist  failure  in  a 
heavy  section  of  a  pressurized  vessel.  The  cut 
penetration  must  cause  a  suitable  increase  in 
stress  due  to  pressurization  to  effect  failure. 

In  addition  to  measuring  shaped  charge 
cuts  and  overpressures,  vessel  overpressure 
and  film/video  data  for  fragment  velocities, 
experience  with  the  pressurization  system  and 
\vith  explosives  personnel  was  gained. 
Problems  with  the  pressurization  system  were 
documented  but  were  not  fully  resolved  this 
trip  to  NSWC.  Measurements  of  shaped 
charge  cross-section  dimensions  became  a 
regular  procedure  in  lieu  of  following  the 
catalog.  The  method  for  preventing  damage 
to  the  pumping  systems  by  a  fragment  jerking 
the  line  was  changed  to  using  the  safety 
anchor  block  as  a  primary  device, 
operational  problems  including  the 


Figure  3.2,  Shaped  Charge  Performance 


The  prior  method  of  explosively  cutting  the  fill  line  caused 
premature  depressurization  of  a  vessel. 


Average  dimensions  and  longitudinal  stresses  for  the  two  vessels  that  burst  are  shown  in 
Table  3.4.  Like  vessel  HI,  vessel  PI  appeared  to  have  burst  because  the  average  of  (Sy+SJ/2 
=  114  ksi  was  approached.  The  failure  mechanism  for  vessel  P2  could  not  be  explained.  The 
groove  appeared  to  be  fairly  uniformly  machined  based  on  32  measurements  more  or  less  equally 
spaced  about  the  circumference  although  the  maximum  LSC  cut  was  .035  inches  greater  than 
average  (measured  since  detonation  was  the  day  prior).  This  caused  us  to  modify  techniques  to 
use  pre-LSC  cut  stresses  generally  less  than  90  ksi  and  to  start  recorders  early  when  that  value 
was  exceeded.  Recorders  and  high  speed  video  were  also  started  early  for  some  of  the  other 
tests. 


Table  3.4 

Vessel  Dimensions  &  Longitudinal  Stresses 


■ 

Average 
O.D.  (in.) 

Average  orig. 
WT(in.) 

Groove  WT  as 
machined  (in) 

LSC  cut,  (in.) 

Average  stress 
(ksi) 

PI 

23.992 

.941 

.240 

.066’ 

112.2 

P2 

24.032 

.959 

.576' 

NA^ 

90.3 

'From  Fig.  3.2  ^Not  applicable,  cut  the  day  prior 

^Includes  measured  LSC  cut 
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3.2  Test  Plans  and  Test  Reports 

Eight  test  plans'^‘^^  were  prepared.  This  includes  a  test  plan  for  each  of  the  six  test  series 
shown  in  Figure  2.1  plus  two  additional  test  plans  which  were  not  conducted  for  lack  of  funds. 
The  deleted  test  series  are  for  Test  Plan  3  and  Test  Plan  4.  Test  Plan  3  was  to  have  been  the 
center  split  burst  of  four  cylindrical  vessels  at  3.5’  HOB.  The  vessels  were:  two  at  53  cubic  feet 
volume  with  diameters  of  24  inches  and  36  inches  and  two  of  22  cubic  feet  volume  with  diameter 
of  16  inches  and  36  inches.  One  of  these  vessels  was  burst  in  Test  Plan  6A  when  TP#3  was 
canceled.  Test  Plan  4  was  to  have  been  the  test  of  four  composite  spheres,  2  with  vertical  split 
and  2  with  horizontal  split.  One  of  these  was  also  burst  in  TP  #6A  testing. 

Each  test  plan  included  items  typical  of  the  following: 

A  Vessels  -  stresses,  inspection,  testing,  machining 

A  Schedule  -  grooves,  pressures,  shape  charges,  other  explosives. 

A  Site  preparation 

A  Instrumentation 

•  Expected  values 

•  Transducers  -  vessel  &  arena 

•  Make/break  wire  instnunentation 

•  Data  recording 

•  Calibration 

A  Support  equipment 

A  Photo  support 

A  Safety 

•  Explosives 

•  High  Pressure  Gas 

•  Vessel  Fragments:  blockhouse  penetration  and  range 

•  Blast  and  test  stand 

A  Post  Test 

A  Procedures/checklists 

A  Drawings 

Five  test  reports^^'^’  cover  testing  and  provide  much  the  same  kind  of  information  foimd 
herein  for  the  entire  test  program. 

3.3  Test  Vessels 

At  program  initiation  a  telephone  search  was  conducted  of  NASA,  Air  Force  and  Navy 
installations  for  vessels  that  could  be  contributed  for  testing.  Many  vessels  were  offered  but  were 
never  used  because  they  would  represent  one  of  a  kind  vessels  to  the  test  program.  Repeat  tests 
or  the  variation  of  some  parameters  which  were  actually  varied  would  be  impossible.  The  largest 
collection  of  similar  vessels  was  found  at  Cape  Canaveral  Air  Station.  Fifteen  of  these  vessels, 
2’  diameter  and  53  cubic  feet,  were  used  in  the  test  program.  A  typical  modification  drawing 
for  these  vessels  is  provided  in  Appendix  A. 
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In  addition  two  composite  vessels,  obtained  from  Lewis  Research  Center,  were  burst. 
These  were  stainless  steel  lined  spheres,  overwrapped  with  Kevlar/epoxy  that  were  rejected  for 
space  shuttle  use.  Two  additional  53  cubic  feet  vessels  were  fabricated  for  testing  so  that  a 
variation  of  diameter  with  the  same  internal  volume  would  be  available  for  testing.  A 
specification  drawing  for  these  vessels  is  included  in  Appendix  A. 

Data  on  all  three  vessels  are  provided  in  Table  2.5  in  the  Summary  Section,  Section  2. 
The  section  also  contains  tabulations  of  some  of  the  most  pertinent  data  relative  to  burst 
conditions. 

3.3.1  Preparation  for  Testing 

Two  penetrations  were  added  to  the  cylindrical  vessel:  pressure  transducers  at  1/4  length 
from  supply  end  (1/4  length  or  more  to  split)  and  for  a  thermocouple.  One  or  more  grooves  were 
cut  in  the  circumferential  direction  in  all  cylinders  and  four  axial  grooves  were  cut  in  vessels  6A- 
1  and  6A-4  of  Test  Plan  6A.  For  TP  #5,  machining  of  vessels  was  switched  from  a  central 
Florida  shop  where  a  lathe  was  used  for  groove  cutting  to  Wright  Patterson  Air  Force  Base  where 
a  computer  numerically  controlled  (CNC)  six  axis  milling  machine  was  recently  installed.  This 
provided  better  accuracy  and  made  it  possible  to  cut  the  longitudinal  grooves. 

Prior  to  cutting  grooves,  an  ultrasonic  volumetric  inspection  was  made  of  the  groove  area 
to  ensure  against  defects  in  the  area  to  be  machined. 


On  site  but  away  from  the  arena,  one  of  the  Preliminary  Pneumatic  Test  vessels  had 
standoffs  built  up  on  the  vessel  surface  using  tape.  These  were  for  the  LSCs  detonated  to 
measure  penetration  depth.  Heavier  pieces,  to  locate  the  LSCs  laterally,  were  also  applied. 
Measurements  of  the  LSCs  showed  a  smaller  cross-section  than  the  catalog  and  this  necessitated 
adding  thicknesses  of  tape  to  the  sides  of  grooves  for  proper  LSC  alignment. 

3.3.2  Groove  designs/dimensions/stresses 


Stress  analyses  used  the  remaining  area  available  to  resist  tensile  stress  after  reduction  in 
area  by  machining  and  detonation  of  a  shaped  charge.  These  equations  reduce  to: 


Sr  = 


(Ry 


-1 


(eq.  3.1) 


Sc=P 


2\ 


(eq.  3.2) 
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where: 


=  longitudinal  stress 

S(~  =  circumferential  stress 

P  =  vessel  pressure,  psig 

=  outer  radius  remaining  at  groove,  inches 

i?.  =  iimer  radius,  inches 

It  was  desired  to  hold  the  stress  level  below  yield  at  pressurization  conditions,  prior  to  the 
cut  or  penetration  caused  by  the  LSC  detonation  (Section  3. 1.2.3).  In  one  instance,  vessel  1-1, 
pre-detonation  stress  was  held  quite  low  to  increase  the  safety  factor  for  handling  and  machining. 
Table  3.5.  Of  course  this  required  a  larger  LSC. 

On  24  January  92  a  burst  of  vessel  5-2  was  attempted  at  3500  psig  internal  pressure  by 
detonating  a  25  gr/ft  lead/RDX  LSC.  The  vessel  failed  to  burst  at  a  calculated  average  stress  of 
1 16  ksi  using  .050"  typical  penetration.  After  depressurization,  measurements  confirmed  that  the 
typical  penetration  had  been  realized.  Yet  the  vessel  failed  to  burst  at  a  calculated  average  stress 
that  was  higher  than  that  of  every  vessel  burst  except  vessel  1-4.  This  was  also  the  first 
attempted  burst  of  a  vessel  whose  groove  had  been  machined  at  WPAFB  using  a  six  axis  CNC 
milling  machine.  It  was  concluded  that  since  the  vessel  was  machined  to  a  greater  accuracy  with 
less  local  variability  that  one  of  the  following  was  true: 

1. )  the  average  wall  thickness  using  32  measurements  more  or  less  equally  spaced  on 

the  lathe  machined  vessels  was  not  accurately  enough  known  or 

2. )  on  vessels  burst  previously  the  weak  spots  transferred  the  load  to  the  strong  spots 

to  a  greater  degree  than  was  thought,  causing  burst  at  a  lower  average  stress  level. 

The  test  of  Test  Plan  5  vessels  was  abandoned  for  January  92  and  only  vessel  PC  was 
burst.  Changes  to  burst  initiation  for  TP  #5  vessels  included: 

1. )  increasing  the  pressure  to  3600  psig 

2. )  switching  LSCs  from  25  gr/ft  lead/RDX  to  50  gr/ft  Al/RDX. 

The  aluminum  LSC’s  are  physically  smaller  than  50  gr/ft  lead/RDX  and  with  .005  inches  filed 
from  each  edge  of  the  LSCs  by  the  manufacturer.  Explosive  Technology,  they  fit  in  the  grooves. 
This  avoided  the  cost  of  sending  the  vessels  back  for  additional  machining. 

The  first  vessel  to  be  burst  in  the  continuation  of  Test  Plan  #5  in  June  1992  was  vessel 
5-1.  Since  the  aluminum/RDX  shaped  charges  require  a  greater  standoff  distance,  layered 
"gaffer’s"  tape  was  used  to  build  up  Ae  surface  and  dlow  the  LSC  to  cut  through  the  tape  and 
into  the  steel.  Tape  was  placed  intermittently  with  two  inch  strips  and  two  inch  open  areas.  On 
LSC  detonation  there  was  a  delay  of  13  ms  until  vessel  burst.  Later  inspection  of  fragments 
revealed  that  there  was  no  LSC  penetration  through  the  two  inch  lengths  of  tape.  For  succeeding 
TP5  bursts  the  standoff  tape  was  changed  to  3/4  inch  lengths  spaced  two  inches  apart  and  the 
delays  were  reduced  to  near  nothing  (Table  5.2). 
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Table  3.5 

Steel  Vessel  Groove,  Stress  Data 


Knowledge  gained  from  Test  Plan  5  provided  additional  confidence  in  using  a  pre¬ 
detonation  stress  that  was  very  close  to  yield.  However  the  practice  of  starting  tape  recorders  and 
high  speed  video  15  minutes  prior  to  estimated  burst  pressure  being  achieved  was  continued. 

Stress  analyses  for  Test  Plan  6A  required  a  knowledge  of  the  penetration  of  Al/RDX  LSCs 
into  both  SA516  and  SA372.  Samples  from  vessel  materials  were  sent  to  Explosive  Technology 
who  provided  data^*  for  Table  3.6.  Samples  of  SA516  were  also  furnished  to  the  Materials  Lab 
at  Wright  Patterson  Air  Force  Base  for  material  properties.  The  new  data^^  are  combined  with 
data  from  Table  3.2  to  form  Table  3.7. 


Table  3.6 

Penetration  of  AL/RDX  LSC  at  Optimum  Standoff 


Charge  wt.  grain/ft 

Penetration 

into  SA516,  in.  into  SA372,  in. 

50 

.098 

.069 

75 

.127 

.101 

100 

.160 

.123 

200 

.270 

.196 

325 

.365 

.315 

400 

.450 

.336 

Table  3.7 

Measured  Strengths  of  Steels  in  TP  #6A 


Measured  Strength,  ksi  j| 

SA516 

SA372  1 

_ _ 

Su 

Sy 

Su 

P’ 

RP^ 

48.6 

52.4 

76.4 

76.5 

98 

130.3 

48.7 

53.7 

76.6 

77.3 

97 

113.7 

Opiate  ^Rolled  Plate  || 
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For  the  multi-fragment  vessels  in  TP6A  no  delay  between  detonation  and  failure  could 
be  tolerated  because  the  delay  at  one  groove  would  exceed  another,  thereby  altering  breakup.  To 
avoid  a  burst  delay  and  use  minimum  LSC  core  loads  a  pre-detonation  stress  slightly  below  yield 
and  post-detonation  stress  slightly  above  ultimate  tensile  was  desired.  This  determined  LSC  sizes 
which  are  shown  in  Table  3.5  along  with  all  the  steel  vessel  stresses.  The  post-detonation  stress 
in  the  longitudinal  groove  of  vessel  6A-1  is  for  one  groove  only.  The  LSCs  failed  to  cut 
adequately  in  the  other  grooves  as  discussed  in  Section  3.7.  The  stresses  in  the  circiunferential 
grooves  of  vessel  6A-1  and  6A-4  were  not  continuous.  The  circumferential  grooves  were 
interrupted  so  that  the  axial  grooves  could  be  continuous.  A  firing  scheme  using  8  detonators 
was  used  for  these  vessels.  Two  LSCs  in  longitudinal  grooves  were  extended  6"  beyond  grooves 
for  detonators  at  angles  180  degrees  apart  at  both  the  east  and  west  end.  In  between  2  each  1/4 
circumference  length  LSCs  were  extended  beyond  their  grooves  and  one  detonator  used  for  both 
shaped  charge  segments. 

3.4  Nitrogen  Pressurization  System 

System  design  requirements  included  7500  psig  ouqjut,  capability  of  protecting  expensive 
components  and  operation  remotely  from  the  blockhouse.  The  pressurization  system^®  which 
evolved  is  designed  to  free  flow/pump  nitrogen  from  a  6000  psig,  240  cubic  feet  storage  vessel 
into  the  test  vessel.  An  Arnold  Engineering  Development  Center  (AEDC)  vessel  was  recertified^’ 
and  shipped  to  the  test  site  to  provide  gas  storage.  The  storage  vessel  is  semi-buried  and 
protected.  Commercially  available  repressurization  used  a  transporter  which  piunped  liquid 
nitrogen  through  a  vaporizer  and  into  the  systems.  Two  Haskel  gas  boosters  in  parallel  provide 
pressures  to  9000  psig  in  the  presence  of  zero  leaks.  Booster  drive  gas  was  initially  bootstrapped 
from  the  storage  volume.  Modifications^^’^®  made  to  the  system  include  provisions  for  a  high 
flow,  low  pressure  compressor  to  furnish  drive  air,  and  was  thereafter  always  used.  A  system 
diagram  is  furnished  in  Appendix  B. 

In  operation,  equalization  flow  is  used  xuitil  the  test  vessel  achieves  storage  vessel  pressure 
at  which  time  the  boosters  are  started.  After  starting,  the  boosters  draw  supply  air  from  the  high 
flow  compressor,  located  behind  the  blockhouse  from  the  arena,  which  has  been  idling.  Test 
vessel  pressure  is  read  with  a  transducer  and  on  a  gage  at  the  pumping  system.  Closure  of  a 
valve  at  the  pumping  system  shuts  off  the  flow  and  permits  both  items  to  be  exposed  to  the  same 
pressure.  The  valve  is  operated  to  get  a  2  to  3  point  calibration  within  a  few  himdred  psi  of  test 
pressure.  This  calibration  thus  includes  zero  and  sensitivity  shift  effects  due  to  solar  and 
compression  heating  and  provides  an  accurate  final  pressure. 

Pressurization  of  the  2450  psi  rated  vessel  to  an  attempted  7500  psi  was  more  difficult 
than  expected  due  to  leaks  at  the  2-3/4  nozzle  plug  weld  at  one  end.  The  leak  was  found  at  low 
pressure  and  an  imsuccessful  repair  attempted.  Test  pressure  of  7125  psig,  where  inflow  matched 
leaks,  was  accepted.  Approximate  pressurization  rates  for  this  pump  up  at  different  pressures  are 
given  in  Table  3.8.  For  equalization  flow,  the  240  ft^  storage  vessel  was  initially  at  5430  psig. 
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Table  3.8 

Pressurization  Rates  for  Vessel  1-4 


Vessel  Pressure 

Pressurization 

Flow  Type 

Storage  Vessel 

(psig) 

Rate  (psi/min) 

Pressure  (psig)  | 

250 

Equalization 

170 

4480 

IHSEH 

115 

4210 

40 

Booster 

4070 

3930 

3760 

10 

3520 

Photos  of  the  installation  are  included  in  the  photo  section,  Section  3.7.  Figure  3.6  shows 
the  storage  vessel  in  a  partial  excavation  with  output  plumbing  attached.  In  Figure  3.7  the  vessel 
has  been  covered  with  approximately  12”  of  earth  and  surplus  pressure  vessels  having  1.2  inch 
wall  thickness.  Two  more  feet  of  earth  completes  the  installation.  Figure  3.8  shows  the  free 
flow/booster  pumping  system  that  is  controlled  by  an  operator  in  the  blockhouse.  The  system 
is  skid  mounted  for  removal  when  not  in  use  for  pressure  vessel  testing.  Figure  3.9  shows  the 
storage  vessel  being  pressurized. 

3.5  Site  for  Pneumatic  Burst  Tests 

The  Preliminary  Pneumatic  Test  and  all  pneumatic  burst  tests  were  conducted  at  the  Naval 
Surface  Warfare  Center’s  (NSWC)  experimental  test  area  at  Dahlgren,  VA.  The  Center  has 
personnel  experienced  in  explosive  detonation  and  blast  data  recording  from  small  up  to  very 
large  charges  of  high  explosive.  High  speed  motion  picture  coverage  is  available  with  multiple 
cameras  and  hardened  camera  shelters.  Heavy  duty  handling  equipment  is  available  such  as 
cranes,  fork  lifts,  payloaders,  etc.  A  variety  of  transducers,  tape  recorders  and  timing  controls 
are  available  for  testing.  A  hardened  blockhouse  and  instrumentation  room  plus  the  capability 
of  tape  recorder  control  from  a  remote  site  is  available.  This  site  provides  an  already  wired  arena 
in  close  proximity  to  a  blockhouse  which  can  prevent  penetration  of  high  kinetic  energy 
fragments.  Figure  3.3  is  a  sketch  of  the  test  area.  An  isometric  drawing  of  a  pressure  vessel 
installed  in  a  blast  field  arena  at  NSWC  is  shown  in  Figure  3.4.  (Drawing  pertains  to  Test  Plan 
2  vessels.) 

3.6  Test  Hardware 

The  24"  diameter  vessels  were  initially  foimd  to  be  out  of  round  with  a  12.050"  radius 
template  and  feeler  gages. 
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INSTRUMENTATION  ROOM 


TRANSDUCER  CABLES 


22  FT.  LONG 
TEST  VESSEL 


Two  roller  assemblies  were  fabricated  for  mounting  and  turning  the  vessels  at  Dahlgren. 
These  facilitated  cleaning  preservative  from  grooves,  making  groove  measurements  and  modifying 
standoffs  in  the  field. 

Four  vessel  support  stands  were  fabricated  and  one  was  jury-rigged.  The  fabricated  stands 
were  used  for: 

1.  Center  split,  1/4  length  and  end  cap  split  vessels. 

2.  Composite  vessel  with  vertically  lavmched  fragments. 

3.  Composite  vessel  with  horizontally  laimched  fragments. 

4.  Multi-fragment  vessels. 

NSWC  supplied  height  of  burst  instrument  stands  for  3.5’  HOB.  These  were  a  standard 
type  fabrication  they  used  frequently.  A  HOB  instrument  stand  was  designed  by  General  Physics 
Corporation  and  built  by  an  NSWC  welder  because  of  the  size.  This  was  used  for  HOB 
measurements  at  8.7’  and  14’  height  and  10  ft.  to  22  ft.  distances. 

Two  make/break  wire  stands  were  fabricated  and  one  was  jury  rigged.  The  fabricated 
stands  are: 

1.  Wooden  stands  with  adjustable  height  and  center  line  wires. 

2.  Metal  stand  with  multiple  crisscross  wires. 

A  deflector  stand  was  fabricated  for  use  with  end  caps  of  vessels  5-3  and  5-4.  The 
requirements  for  the  metal  make/break  wire  stand  and  the  deflector  stand  emanated  from  the 
inability  to  predict,  a  priori,  the  angle  that  the  preliminary  composite  vessel  fragment  and  the  end 
cap  fragment  would  depart.  The  metal  wire  stand  was  made  wide  and  square  to  avoid  contact 
and  to  assure  a  wire  strike.  Computations  showed  that  the  end  cap  fragment  would  travel  too 
great  a  distance  if  it  departed  at  an  angle  greater  than  7°  up  elevation.  The  deflector  stand  was 
to  absorb  energy  and  thus  reduce  the  fragment  range  if  it  was  struck.  It  turned  out  that  both 
fragments  appeared  to  have  flown  a  straight  and  true  path  compared  to  their  original  orientation. 


3.7  Photos/Burst  Observations/Site  Rework 

Figures  3.5  through  3.9  were  discussed  earlier.  Figure  3.10  shows  vessel  1-1  ready  for 
testing.  The  legs  to  the  test  stands  are  both  aligned  with  and  normal  to  the  vessel  axis.  The  blast 
typically  blows  both  stands  off  of  their  1"  thick  round  plate  bases  and  may  invert  one  of  them. 
For  vessel  1-2  the  stands  were  inadvertently  placed  at  45°  angles  to  the  vessel  axis.  The 
aftermath  of  such  a  placement  is  shown  in  blast  Figure  3.11  where  both  stands  were  inverted  by 
the  blast.  There  was  some  concern  that  the  blast  asymmetry  may  have  been  increased  but  a  close 
screening  of  the  data^^  didn’t  seem  to  show  it. 
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Figure  3.12  shows  a  12"  diameter  tree  trunk  broken  off  by  the  1/2  vessel  1-4  fragment 
and  displaced  4’  from  the  stump.  This  is  one  of  10  trees  of  4"  to  12"  diameter  which  were 
broken  off  in  woods  approximately  850  feet  from  the  arena  center. 

Figures  3.13  through  3.16  are  all  Test  Plan  2  vessel  shots.  For  this  test  plan  the  legs  to 
the  vessel  stand  were  buried  flush  with  the  groxmd  surface  to  preserve  alignment  and  minimize 
blockage  to  the  gas  outflow.  Figure  3.13  is  a  shot  of  the  14’  HOB  vessel  2-3.  Figure  3.14  is 
a  post-test  shot  of  the  8.7’  HOB  vessel  2-2.  Figures  3.15  and  3.16  are  pre  and  post-test  shots  of 
the  3.5’  HOB  vessel  2-1.  The  latter  figure  shows  scouring  of  the  earth  in  the  center  of  the  arena. 
This  occurred  frequently  in  warm  weather  until  light  steel  plate  was  used  in  the  center  of  the 
arena. 


Figure  3.17  is  vessel  5-1  prior  to  test.  A  steel  plate  under  the  split  region  was  blown  out 
of  the  arena  on  burst.  The  contact  wire  frame  (used  earlier  in  vessel  PC  test.  Figure  3.25)  was 
also  blown  over  as  shown  in  Figure  3.18. 

The  steel  plate  length  was  increased  at  each  end  in  the  east-west  direction  and  the  ends 
bent  90°  and  buried  9"  in  the  ground  for  vessel  5-2  test.  Figure  3.19.  At  burst  the  east  edge  of 
the  plate  was  scooted  about  3"  to  the  west  with  the  center  of  the  plate  bulged.  The  makewire 
stand  was  moved  about  a  foot  by  the  blast  and  rotated  a  few  degrees  to  the  south.  Figure  3.20. 

Vessel  5-4,  shown  in  Figure  3.21,  was  tested  before  vessel  5-3.  The  makewire  stand  was 
changed  to  a  breakwire  stand  by  cutting  sections  of  steel  (bolted  in  a  modification  on  June  15). 
Reinforcement  was  added  to  the  top  of  the  makewire  stand  to  slow  down  the  end  cap  fragment 
in  case  of  penetration.  The  deflection  stand  was  added  and  welded  to  the  makewire  stand  braces 
added  for  vessel  5-2.  At  burst  the  makewire  stand  was  destroyed  as  the  blast  compressed  the 
eastern  part  of  the  entire  stand,  straight  areas  being  compressed  more  than  comer  areas.  Figure 
3.22.  Total  compression  was  approximately  1-1/2  feet.  The  stand  was  also  pushed  by  the  blast 
about  40  feet  which  required  overcoming  8  trailer  stakes  and  3/16"  cables.  Some  were  yanked 
out  and  some  cables  pulled  through  the  clamps.  The  end-cap  fragment  appeared  (by  its  location 
in  the  woods)  to  depart  on  a  line  which  would  be  an  extension  of  the  vessel  axis. 

For  vessel  5-3,  Figure  3.23,  the  makewire  (now  breakwire)  stand  was  cut  loose  and  the 
barely  damaged  deflection  stand  was  located  to  its  original  position.  Locally  available  steel  pipe 
was  added  to  form  a  holding  device  for  a  single  row  of  four  breakwires  using  the  music  wire. 
This  was  placed  6"  below  center  so  that  another  device  could  be  used  on  centerline.  A  pipe, 
supported  in  two  places,  was  added  to  the  centerline  of  the  vessel  to  break  graphite  pencil  leads 
at  a  location  (greater  than  6  feet  from  ground  zero)  where  the  pipe  would  arrive  before  the  blast 
wave.  The  graphite  leads  were  supported  in  a  wood  framework  by  modifying  the  wood 
makewire  stand.  Pre-test  calculations  had  shown  that  a  pole  attached  to  this  fragment  could  be 
a  difficult  design  problem  due  to  potential  bending  due  to  acceleration.  However,  the  effort  was 
minimal  and  NSWC  favored  pencil  leads  where  useable  due  to  the  sharp  break  obtained.  The 
devices  can  be  seen  in  the  photo. 


24 


Upon  burst  the  structure  added  to  the  front  of  the  stand  was  destroyed,  Figure  3.24,  by 
being  bent  outward  from  fragment  flight  path.  The  blast  also  scooted  the  structure  about  10’  to 
the  west  and  rotated  the  eastern  edge  about  20“  to  the  north.  The  wooden  legs  of  the  support  for 
the  pencil  leads  were  broken  and  the  east  end  of  the  base  was  rotated  about  50“  to  the  north. 
Sharp  break  data  were  obtained  from  the  graphite  leads.  However  the  graphite  lead  data  showed 
a  considerably  lower  velocity  than  the  wire  data.  The  lower  velocity  can  only  be  explained  by 
a  bending  of  the  pole  as  the  endcap  accelerated.  The  graphite  lead  data  were  not  included  in 
Section  6. 

Figure  3.25  is  a  pre-burst  shot  of  the  preliminary  composite  sphere,  vessel  PC.  Fragments 
were  laimched  vertically.  The  bottom  portion  fragmented  into  very  small  pieces  upon  striking 
the  steel  plate.  The  upper  fragment  was  launched  through  a  wire  grid  with  crossing  contact  wires 
spaced  a  foot  apart.  Figure  3.26  is  the  second  composite  sphere,  vessel  6A-2.  Fragments  were 
launched  horizontally.  The  left  hemisphere  carried  the  support  with  it,  which  broke  off  at  the 
bottom  steel  plate.  Figmes  3.27  and  3.28  are  both  fragments  from  vessel  6A-2.  The  latter  is  the 
largest  liner  fragment  found. 

Figure  3.29  is  vessel  6A-1  ready  for  test.  No  shaped  charge  was  used  in  the  center 
circumferential  groove.  The  axial  grooves  were  supplied  with  angle  iron  for  lateral  positioning 
of  the  shaped  charges.  Shaped  charge  standoff  was  multiple  strips  of  "gaffers"  tape  to  build  up 
to  .221  inches  (to  nearest  .010"  tape  thickness).  The  tape  was  used  in  1"  widths  with  the  angle 
iron  on  top  (screwed  through  tape  gaps).  Vessel  6A-1  was  supposed  to  be  a  multi-fragment  test, 
however  at  burst  only  the  two  end  caps  separated.  Figure  3.30.  The  center  section  remained  in 
place  and  shifted  about  a  foot.  The  blast  broke  the  shear  pins  at  the  frame  to  support  and  the 
frame  was  driven  to  the  ground.  The  frame  was  distorted  beyond  on-site  repair. 

The  axial  grooves  were  examined  closely  in  vessel  6A-1  after  the  edges  failed  to  part. 
Some  of  the  grooves  were  not  adequately  cut  by  the  LSC  with  the  jet  impingence  on  the  side  of 
some  grooves.  Apparently  the  LSCs,  which  were  very  stiff,  could  not  be  adequately  held  in  place 
by  the  duct  seal  which  was  pressed  on  top  of  it  as  a  holding  device.  At  least  one  groove  was 
adequately  cut  with  the  LSC  cut  located  along  the  groove  center.  The  reason  this  groove  did  not 
split  is  unclear. 

Figure  3.31  is  an  end  cap  fragment  from  vessel  6A-1. 

Figure  3.32  is  the  center-split  3’  diameter  vessel  6A-3  ready  for  test.  Vessel  support  was 
jury-rigged  because  1.)  economics  and  time  did  not  permit  building  3  frames,  2.)  vessel  6A-1 
seemed  unsteady  on  the  inadequately  strengthened  frame  and  3.)  it  was  important  to  reserve  a 
good  frame  for  vessel  6A-4.  The  support  blocked  the  gas  outflow  to  the  extent  that  one  of  the 
lower  12  x  12s  was  foimd  half  way  to  the  woods,  about  400’  away.  This  could  have  affected 
data  asymmetry  but  no  affect  is  apparent.  Section  5.9.  The  center  (cut  edge)  of  each  of  the  two 
fragments  was  seen  to  tip  up  on  the  video  monitor.  This  never  occurred  on  other  center-split 
tests. 
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After  the  unsuccessful  multi-fragment  burst  of  vessel  6A-1  the  second  multi-fragment 
vessel,  6A-4  was  painstakingly  prepared  for  test,  Figure  3.33.  Steel  bars,  under  angle  iron 
positioning  devices,  had  been  prepared  for  precise  stand-offs.  However  the  duct  seal  would  not 
hold  the  stiff  LSCs  in  the  bottom  of  the  space  between  the  angles.  The  angles  sloped  away  from 
each  other  due  to  the  vessel  curvature.  A  combination  of  duct  seal  and  clamping  the  angles 
against  the  LSCs  as  each  mounting  screw  was  tightened  was  finally  successful.  Four 
circumferential  grooves  were  cut  but  no  LSCs  were  placed  in  the  two  center  grooves  because  the 
fragment  velocity  program^^  did  not  anticipate  a  split  circumferentially  except  at  end  caps  when 
side  fragment  are  considered.  However  14  fragments  were  found  at  separate  locations.  Figure 
3.34  is  taken  from  high  speed  film  of  the  vessel  burst.  The  film  appears  to  show  two  axial 
fragments  that  landed  separately  flying  together  as  an  integral  unit. 

The  support  for  vessel  6A-4  incorporated  shear  pins  at  the  vertical  supports  to  permit  the 
frame  to  be  driven  to  the  ground,  Figure  3.35.  This  design  minimized  damage  to  the  vertical 
supports  and  permitted  acceleration  of  the  lower  fragments  with  a  possible  result  of  a  stronger 
blast  wave.  The  side  fragments  in  Figure  3.33  were  intended  to  depart  horizontally  and  clear  the 
instrumented  arena. 

Figure  3.36  is  a  lineup  of  vessel  6A-4  fragments  (rotated  90°). 

3.8  Instrumentation 

Most  of  the  follovving  is  repeated  from  the  Test  Plan  5  test  reporF^  and  is  typical  for  all 
tests.  The  following  parameters  were  measured  and  recorded  dining  the  two  test  series: 

•  vessel  pre-burst  pressure 

•  vessel  pre-burst  gas  (and  wall)  temperatures 

•  vessel  fragment  location  versus  time 

•  nitrogen  pressurization  line  pressure 

•  time  of  shaped  charge  detonation 

•  time  of  each  contact  wire  signal 

•  field  pressure  at  each  transducer  location 

Instrumentation  used  to  measure  the  above  parameters  is  listed  in  Table  3.9. 


Table  3.9 
Instrumentation 


pressure  transducer  (vessel) 

Endevco  8511 A 

1 

pressure  transducer  (field) 

PCB  137A11/21  (HOB) 

PCB  102A02  (ground) 

PCB  102A  (ground) 

thermocouple 

Omega  TJ36-CPSS-18U-12 

Omega  TJ36-CPSS-116U-18 

pressure  gage  (line) 

Ashcroft  45-1377SS-04B 

pressure  transducer  (line) 

Omega  PX303-10KG-5V 
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Figure  3.11,  Stands  Inverted  During  Test  Figure  3.12,  Broken  Trunk  Displaced  From  Stump 


Ground  Scoured  by  Vessel  2-1 


Figure  5.28,  Vessel  6A^2  Fragment 


Figure  3.35,  Vessel  6A-4  Support  After  Test  Figure  3.36,  Vessel  6A-4  Fragments 


The  field  pressure  transducers  were  calibrated  by  NSWC  prior  to  the  test  program.  Many 
of  the  transducers  were  to  be  used  to  measure  expected  pressures  well  below  their  full  scale 
ranges,  as  piezoelectric  transducers  are  not  normally  available  at  pressures  below  50  to  100  psig. 
As  a  result,  a  special  calibration  stratagem  was  devised  to  ensure  the  accuracy  of  transducer 
measurements.  Rather  than  selecting  one  sensitivity  for  the  full  scale  range  (industry  practice  for 
transducers  used  to  measure  pressures  near  full  scale  rating),  the  transducers  were  calibrated  at 
a  range  of  pressures  with  a  different  sensitivity  determined  at  each  pressure.  This  stratagem 
permitted  the  selection  of  a  sensitivity  for  the  actual  pressure  anticipated  and  for  which  the 
recorder  was  ranged.  Calibrating  to  full  scale  rather  than  the  expected  value  could  produce  many 
results  with  poor  signal  to  noise  ratio  on  the  FM  tape  recorder. 

A  line  pressure  gage  was  used  for  calibrating  each  test  vessel  pressure  transducer  and  for 
accuracy  checks  during  pressurization.  The  gage  connection  was  450  feet  away  from  the  test 
vessel  in  a  1/2"  O.D.  pressurization  line.  At  low  flow  rates  and  high  pressmes  (booster 
operation)  there  is  a  negligible  pressure  drop  in  the  line.  At  low  test  vessel  pressures  and 
appreciable  mass  flow  rates  (equalization  flow)  there  was  a  measurable  pressure  drop  in  the  line. 
A  spot  check  on  the  transducer  calibration  (needed  due  to  potential  thermal  drifts)  was  generally 
obtained  when  flow  was  stopped  as  maximum  pressure  was  reached.  This  gage  was  calibrated 
by  the  Johns  Hopkins  Applied  Physics  Laboratory  after  testing  in  July  1991.  It  was  recalibrated 
by  High  Purity  Systems,  Titusville,  Florida  in  October  1993. 

Vessel  work  instruments  were  also  used.  General  Physics  Corporation  used  a  combination 
of  ultrasonic  wall  thickness  (UT)  measurements  and  depth  micrometer  to  determine  remaining 
groove  thickness  of  the  SA372  vessels.  The  UT  device  used  was  a  Panametrics  model  22  with 
a  Microscan  number  Ml  12-RM  transducer  having  a  range  of  .035"  to  3".  WPAFB  used  a  similar 
device.  Two  gage  blocks  of  SA372  from  vessel  material  were  furnished  by  WPAFB  to  General 
Physics  Corporation.  For  measuring  LSC  detonation  cuts  in  vessels  General  Physics  Corporation 
used  a  Starrett  model  449AZ-3R  depth  micrometer  with  the  non-rotating  narrow  blade  ground 
thinner  as  required  to  fit  in  slots  cut  by  the  smaller  LSCs. 

3.8.1  Instrumentation  Channels 

Approximately  45  instrumentation  channels  were  used  in  each  test  for  field  pressure 
transducers  plus  a  vessel  pressure  transducer,  a  pressurization  line  pressure  transducer,  vessel  and 
fragment  make/break  wires.  Some  strain  gage  chaimels  were  recorded  during  multi-fragment 
testing,  but  since  the  vessels  did  not  burst  prematurely  the  data  were  not  analyzed.  Field  and 
vessel  pressures  and  timing  data  were  all  recorded  on  Fairchild  Model  9  FM  tape  recorders  at  a 
speed  of  60  inches/s  for  a  frequency  response  of  40  khz. 

Figure  3.37  is  a  block  diagram  illustrating  the  sequence  by  which  data  are  measured, 
recorded,  filtered  and  prepared  for  output  display  and  analysis.  The  block  diagram  is  from  an 
NSWC  paper^^  on  the  pressure  vessel  burst  testing. 
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Figure  3.37,  NSWC  Data  Conditioning  Block  Diagram  [Courtesy  NSWCJ 


Photography 

Typical  photographic  support  requirements  requested  of  NSWC  are  as  follows: 

2  to  5  high  speed  film  cameras.  Typically  2  were  used  for  dual  fragment  tests  and  5  for 
the  multi-fragment  tests.  Camera  placement  for  dual  fragment  tests  varied  and  was 
unsatisfactory  for  TP  #5  endcap  vessels.  Desired  camera  placement  for  TP  #6 A  is  shown 
in  Figure  3.38.  (However  no  film  was  furnished  from  die  two  cameras  at  position  A). 
Frame  rates  of  1 000  to  2000  frames  per  second  (FPS)  were  typically  used. 

400  and  200  FPS  video  cameras  were  used  except  in  TP  #6A.  The  400  FPS  camera  used 
a  field  of  view  of  about  60  feet  and  the  200  FPS  camera  approximately  300  feet  at  the 
vessel  distance.  These  cameras  were  moimted  in  a  tower  above  the  blockhouse  and  were 
not  used  for  TP  #6A  because  fragments  would  be  launched  in  the  blockhouse  direction. 

2-4  standard  video  cameras  for  an  overall  view 

Still  photos  before  and  after  vessel  bursts 

Photo  data  reduction  for  fragment  velocity 

Composite  video  tapes,  selected  still  photos  and  viewgraphs 


3.10  Computer  Software 

Over  90  different  computer  programs  were  used  or  written  as  part  of  the  Burst  Test  Study. 
This  includes  the  following  programs  written  by  others:  Burst,  Fortran  compiler,  Lotus  123©, 
Pvsafe-d,  Sigma-plot©,  Traj,  WordPerfect©,  WP  Presentations©,  Microsoft  Excel©,  Wordscan©. 
It  also  includes  over  80  Fortran  programs  and  more  than  35  variations  of  the  80+  Fortran 
programs  that  were  written  for  use  and  a  scientific  subroutine  package,  written  in  Fortran,  that 
was  purchased  for  use.  Several  programs  were  often  used  in  preparing  files  for  input  to  the 
plotting  program  as  shown  in  the  block  diagram  of  Figure  3.39  prepared  for  Test  Plan  5  report. 
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TPS  BURST  SOFTWARE/FiLES 
FOR  BLAST  OVERPRESSURE 


REGRESSION 


PLOTTING 


CONTOUR 


Y 

REGCK7 

OR 

REGCK8 
(If  offset)* 


3577.DAT 
(  =shot#) 


3577.RG7 

OR 

35ZZ.RG8 


P35ZZ.CON 


PC3SZZ.SPG 


P3877.F3 
col  33-47 

F1  -  Col  1-22  to  plot  measured  test  (symbol)  data 
F2  -  Col 23-30  to  plot  measured  test  (symbol)  data 
F2  -  Co!  31,  32  not  plotted,  for  reference  only 
F3  -  Col  33-45  to  plot  calculated  data  on  angles  from  strongest  array 
F3  -  Co! 46,  47  to  plot  calculated  For  I  for  high  explosive 

*  Some  vessels  offset  1  foot  from  0"-180* 


Figure  3.39 


4.0  BLAST  OVERPRESSURE  MODELING 


Eqviations  and  techniques  used  for  estimating  burst  pressures  or  for  comparing  to  TNT  are 
included  in  this  section.  The  equations  and  techniques  have  also  been  incorporated  into  a 
computer  program,  BLAST  (discussed  in  more  detail  in  a  workbook''^  written  for  the  Air  Force 
and  Navy). 

4.1  Vessel  Stored  Energy 

The  following  equation  gives  the  isentropic  energy  released  by  the  failure  of  a  vessel 
containing  a  volume  of  ideal  gas,  Vi,  at  a  pressure  of  P,.  Pj  is  the  surrounding  atmospheric 
pressure,  y  is  the  specific  heat  ratio: 

(eq.  4.1) 


Multiplying  the  above  energy  in  foot  pounds  by  a  conversion  factor  yields  the  TNT 
equivalence  in  pounds  TNT.  All  the  TNT  equivalence  calculations  in  this  report  use  the  ideal 
gas  energy  equation  and  a  conversion  factor  of  1.545  x  10*  ft  lbs  per  lb  TNT  as  used  by  Kinney 
and  Graham'*.  A  literature  survey'  as  part  of  the  Burst  Test  Study  found  values  in  use  from  1.4 
X  10*  (by  Baker^)  to  1.6  x  10*  (by  the  U.S.  Air  Force^''). 


Ideal  gas  behavior  is  considered 
adequate  for  most  low  pressure  applications 
(1500  psi  or  less).  The  ideal  gas  assumption 
for  high  pressure  ruptures  gives  expansion 
energies  that  can  be  unrealistically  high  due  to 
compressibility  effects.  For  example,  using 
Figure  4.1,  the  real  gas  energy  for  the  highest 
pressure  burst  at  7125  psi  and  133°F  is  less 
than  the  ideal  gas  energy  by  15%  to  20%. 

Figure  4.1  was  published  as  Figure  4  by 
Wiedermann*  who  used  van  der  Waal’s 
equation  of  state  for  evaluating  the  energy  of 
the  integral  of  PdV.  (Wiedermann  uses  1.426 
X  10*  ft.  lbs  per  lb.  TNT  as  a  conversion 
factor.) 

The  kinetic  energy  imparted  to  the  fragments  reduces  the  energy  available  to  create  a  blast 
wave  and  should  be  subtracted.  This  energy  is  not  subtracted  in  this  report.  The  kinetic  energy 
is  computed  and  is  presented  in  Table  6.3  in  Section  6.3  as  a  percentage  of  the  isentropic 
expansion  energy. 


Figure  4.1,  TNT  Equivalency  for 
Pressurized  Vessels  Filled  with  Nitrogen 


Y-1 


1- 
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4.2  Blast  Overpressure  Equations 


Some  authors  use  the  one-dimensional  shock  tube  equation  for  computing  the  blast  wave 
initial  shock  overpressure  at  the  vessel  surface.  This  equation  is  provided  as  equation  4.2  and 
may  be  obtained  from  Liepman  and  Roshko^^  by  renaming  variables.  Efforts  to  measure 
pressures  close  to  the  vessel  to  approach  the  initial  shock  overpressure  were  unsuccessful  due  to 
side  loading  as  the  venting  fragment  departed.  No  measurements  taken  in  the  test  program  show 
the  shock  tube  equation  to  be  inappropriate. 


Pso=Pl 


(Y^-1) 


1- 


f  p  ^ 

_ ££  _  j 

p. 


fp  I 

2y^+(Yo+1) 

SO  J 
<  ^0  } 

2y« 


(eq.  4.2) 


where  Pj 


a 

a 


O 

1 


Po 

y 


initial  vessel  pressure,  psia 
initial  shock  pressure,  psia 
ambient  soxmd  velocity 
sound  velocity  in  vessel  gas 
ambient  pressure 

specific  heat  ratio,  either  ambient  (yj  or  vessel  gas  (y,J 


The  shock  velocity  may  also  be  found  in  Liepman  &  Roshko  and,  by  renaming  variables, 
becomes: 


2/Y, 

P  o 


(eq.  4.3) 


One  data  point  was  obtained  for  checking  the  actual  versus  computed  shock  velocity.  On 
1/26/91  a  2450  psi  rated  pressure  vessel  was  burst  at  5425  psig.  One  gage  was  mounted  on  a 
stand  8"  from  the  vessel  surface  and  axially  at  the  groove.  Another  gage  was  located  on  the 
ground  at  6.1'  from  the  groove  and  edge  of  the  vessel.  The  8"  gage  was  destroyed  when  the 
stand  was  carried  away,  however  shock  arrival  time  was  recorded.  The  difference  in  arrival  times 
provides  an  average  shock  velocity  between  the  two  points  of  1900  fit/sec  compared  to  a 
calculated  velocity  of  2230  ft/sec  at  the  vessel  which  then  slows  as  the  shock  wave  advances. 

A  graph  of  initial  overpressure  and  velocity  versus  vessel  pressure  using  the  equations 
cited  is  provided  as  Figure  4.2.  The  values  were  computed  for  air  in  the  vessel  vrith  both  vessel 
gas  and  ambient  air  at  60°?. 
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Kinney  and  Graham'°  provide  Bode 
type  equations  (their  6-2  and  6-12)  for  the 
incident  shock  overpressure  and  positive 
impulse.  These  equations  are  used  to 
compute  the  TNT  equivalence  overpressures 
in  the  plots  which  are  furnished.  Section  5.1. 
The  equations,  after  conversion  to  English 
units  as  used  in  the  Burst  Test  Data  Reports 
furnished  by  NSWC,  are  provided  as 
equations  4.4  and  4.5.  Both  equations  are 
based  on  chemical  explosions. 


Figure  4.2,  Initial  Shock  Overpressure 
and  Velocity  V5.  Vessel  Pressure 


Where: 


(eq.  4.4) 


A 


4.756 


1  + 


\4 


.580 


/  7 

1+  — ^ 
1,3.907 


(TNT)^  (eq.  4.5) 


P°/Pa  =  ratio  of 
shock 
overpressure  to  ambient  pressure 

I/A  =  impulse  per  unit  area,  psi-ms  (Ref.  10  considers  impulse  to 
have  units  of  force  x  time  rather  than  pressure  x  times) 


and: 


Z  =  scaled  distance 

D/(TNT)*^ 

D  =  distance,  ft 

TNT  =  TNT  equivalency,  lbs. 


The  term  (TNT)’^  is  factored  into  the  Kinney  &  Graham  equation  to  produce  equation 
4.5  so  that  the  result  will  not  be  a  scaled  value.  Hopkinson’s  scaling  law^®  permits  the 
comparison  of  overpressures  of  large  and  small  blasts  (by  comparison  to  one  pound  TNT)  by 
dividing  distance  by  the  cube  root  of  TNT  equivalence.  Impulse  may  also  be  compared  if  the 
impulse  is  also  scaled.  Scaled  values  are  used  in  Sections  5.2  and  5.13. 
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4.3  Methods  from  references 


Baker^  uses  P^^,,  Section  4.2,  as  a  starting  point  in  calculating  the  overpressure  due  to  a 
bursting  pressure  vessel.  Baker’s  curves  (based  on  one-dimensional  hydrocode  calculations), 
assume  sudden  vessel  wall  disappearance,  hence  the  theoretical  predictions  are  typically  higher 
than  can  be  achieved  in  a  real  vessel  burst  since  a  finite  time  is  required  for  the  gas  to  exhaust. 
Baker’s  curves  were  used  for  estimating  some  of  the  overpressures  for  test  planning  and  are  used 
to  compare  test  results  for  two  type  bursts  with  calculated  values.  Section  5.15. 

Held^  presents  a  method  for  estimating  shock  overpressure  using  a  TNT  explosion  located 
at  a  greater  distant  than  the  bursting  pressure  vessel.  Held’s  method  is  based  on  the  fact  that  the 
impulse  from  a  bursting  vessel  is  the  same  as  that  from  a  TNT  detonation  and  is  based  on 
equations  reduced  from  Baker’s  curves.  Like  Baker’s  curves,  a  spherical  vessel  with  walls  that 
suddenly  disappear  is  assumed.  This  approach  was  not  used  for  burst  calculations. 

4.4  Reflected  Shocks 

It  is  characteristic  of  an  explosive  shock  front  that  amplification  of  the  overpressure  occurs 
when  a  nearby  reflecting  surface  exists.  The  amplification  is  typically  termed  a  reflection  factor 
or  reflection  coefficient.  Kinney  and  Graham’®  develop  equations  for  reflection  coefficients  and 
provide  limiting  values  of  2.0  for  very  weak  shocks  and  8.0  for  very  strong,  normal  (head  on) 
shocks.  They  also  provide  data  for  oblique  shocks  and  reflection  factors  due  to  Mach  stem 
formation  (grazing  incident  shock  and  reflected  shock  joined  at  the  reflective  surface).  Figure 
5.8  of  Kinney  and  Graham  is  reproduced  here  as  Figure  4.3  (shock  angle  sketch  was  added  to 
the  figure).  The  Mach  number  (ratio  of  shock  velocity  to  sound  velocity)  is  related  to  blast 
overpressure  by  equation  4.6,  thus  strong  blast  waves  are  propagated  at  high  velocity. 

Swisdak®  provides  curves  for  reflection  factor  based  on  overpressure.  Baker’s  workbook^ 
(paragraph  2.3.4)  provides  some  approximate  reflection  factors  for  pressure  vessel  bursts.  Reed^’ 
presents  a  curve  (his  Figure  2)  of  HOB  yield  factor  versus  scaled  burst  height  and  is  reproduced 
here  as  Figure  4.4  (English  units  added).  Reed’s  software  also  computes  an  overpressure-distance 
function. 

Shock  wave  overpressures  close  to  an  explosion  can  be  very  high  as  shown  by  Swisdak’s® 
Figure  3A  which  is  reproduced  here  as  Figure  4.5.  Shock  waves  emanating  from  pressure  vessels 
are  typically  much  less  at  the  source  as  shown  in  Figure  4.2  and  thus  can  be  expected  to  generate 
low  reflection  factors. 
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REFLECTION  COEFFICIENT  (A) 


4.5  Pretest  Calculation  of  Burst  Overpressures 

For  the  initial  vessel  burst,  vessel  P-2,  TNT  equivalence  was  used  to  estimate 
overpressures.  The  data  results  were  a  fraction  of  the  recorder  scale  and  rescaling  was 
accomplished  for  the  following  burst  which  was  to  occur  at  a  lower  burst  pressure.  To  prepare 
for  Test  Plan  #1,  the  overpressures  from  vessels  P-2  and  P-1,  4700  psig  and  3250  psig  burst 
pressure,  were  plotted  versus  distance  and  extended  toward  the  vessel  to  the  nominal  vessel 
pressure.  Starting  at  new  vessel  burst  pressures,  new  overpressure  versus  distance  lines  were 
drawn  parallel  to  the  original  lines.  This  resulted  in  good  recorder  scaling  for  the  following  test 
series. 


New  estimates  for  3.5  ft  HOB  burst  for  the  following  test  series  #3  were  based  on 
normalizing  overpressure  data  by  converting  measurements  from  gage  to  absolute  and  dividing 
by  the  calculated  initial  shock  overpressure,  P^^  from  equation  4.2.  An  average  was  obtained  for 
each  distance  in  the  arena  for  normalized  overpressures  from  5  vessel  bursts.  A  log-log  plot  of 
these  averages  was  used  to  establish  a  base-line  for  converting  to  the  expected  results  from  3500 
psi  and  7500  psi  testing.  With  some  scatter  the  normalized  overpressures  ranged  from  0  24  at 
10'  to  0.153  at  50'. 


Elevated  height  of  burst  tests  were  also 
to  be  conducted  as  part  of  test  series  3  (Test 
Plan  #5  and  completion  of  Test  Plan  #2).  To 
determine  the  necessary  variation  in  height  of 
burst,  the  optimum  heights  of  burst  for  the  10 
to  50  foot  range  of  the  arena  were  found 
using  four  techniques.  The  results  are  plotted 
in  Figure  4.6.  The  methods  used  were:  (1) 
Swisdak’s®  Figures  5B  and  5C,  (2)  Figure  8-8 
of  Kinney  and  Graham’®,  and  (4)  a  HOB 
curve  by  Reed”.  Input  data  were  the  TNT 
equivalent  as  computed  using  the  overpressure 
at  each  distance  from  the  3250  psi  burst  of  the 
preliminary  test^'’.  It  was  decided  to  use  an 
optimmn  height  for  the  mid-range  of  the  arena 
of  22  feet.  This  required  a  height  of  14  feet 
to  cover  the  requirement  of  3  of  the  4 
techniques  and  19  feet  to  cover  technique 
number  3.  For  practical  reasons  it  was 
decided  to  limit  the  height  to  14  feet. 


Figure  4,6,  Peak  Overpressure 
Versus  Distance  Comparison 
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Two  methods  were  used  to  estimate 
overpressure  variation  with  HOB  in  order  to 
determine  if  the  overpressure  variation  would 
be  large  enough  to  measure.  One  used  the 
shock  overpressure  ratio  and  the  other  used 
the  shock  Mach  number.  These  two  methods 
are  not  necessarily  independent  because  both 
use  a  shock  incident  angle  and  the  shock 
Mach  number  and  overpressure  ratio  are 
related  by  an  equation®  given  below  (equation 
4.6).  Figure  4.7  is  a  plot  of  the  equation  and 
pressure  vessel  burst  overpressure  data  versus 
average  shock  velocity  between  two  stations 
prior  to  the  pressure  measuring  point.  It 
shows  reasonable  correlation,  with  the  better 
scaled  data  for  vessel  P-1  showing  better 
correlation  than  P-2  data. 
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(eq.  4.6) 


Shock  Velocity  (ft/s) 


Figure  4. 7,  Peak  Overpressure  Versus  Shock  Velocity 
for  3250  psi  Vessel  P-1  and  4700  psi  Vessel  P-2 


where  p. 


So 


Y 

P 

P 


0 


=  shock  velocity  (ft/s) 

=  sonic  velocity  ahead  of  shock  front  for  air  at  ambient  temperature 
=  ratio  of  specific  heats 
=  peak  overpressure 
=  ambient  pressure  (atmospheric) 


The  shock  Mach  number  was  estimated  from  Preliminary  Pneumatic  Burst  Test  data^''. 
These  data  were  obtained  by  time  of  arrival  differences  between  two  distances  in  the  array.  An 
estimate  of  Mach  1.2  was  used  at  5'  to  10'  from  the  vessel  and  Mach  1.1  at  50'. 


Incident  overpressure  data  were  estimated  using  incident  angles,  based  on  3.5'  HOB. 
Figure  4.3  (Figure  5-8  of  Kiimey  &  Graham'®)  and  the  3500  psi  burst  data  found  earlier.  The 
resultant  overpressures  ranged  from  10.5  psi  at  10  feet  to  3.9  psi  at  50  feet  on  a  log-log  plot  with 
slight  scatter. 

Overpressure  data  were  then  computed  by  the  reverse  of  the  procedure  in  the  preceding 
paragraph  with  incident  angles  for  3.5',  8.7'  and  14'  HOBs. 
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Figure  4.8 

Estimated  Overpressures  for  Height  of  Burst  Testing 

Overpressure  data  were  then  computed  using  the  same  incident  overpressures,  the 
overpressure  ratios  (to  ambient)  and  the  incident  angle  using  Figure  13a  of  Swisdak’s’  report. 
The  results  are  plotted  in  Figure  4.8.  The  two  methods  yielded  similar  results.  The  expected 
difference  in  overpressures  from  3.5'  to  14'  were  deemed  to  be  adequately  measurable  at  22' 
distance  and  less  so  for  lesser  or  greater  distances. 

Overpressures  for  the  preliminary  composite  sphere,  vessel  PC,  and  the  multi-fragment 
bursts,  vessels  6A-1  and  6A-4,  were  estimated  using  the  methods  of  Baker,  Section  4.3  for 
spheres  with  a  volume  factor  of  2.  The  closest  distance  to  the  cylindrical  vessels  was  used. 
Baker’s  cylinder  correction  factors  were  not  used  but  a  safety  factor  of  2.5  was  used  to  allow  for 
ground  reflection  to  avoid  saturation  recorder  conditions. 
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5.0  BLAST  OVERPRESSURE  TEST  RESULTS 


Graphs  are  provided  to  show  blast  strength  and  patterns.  Most  of  the  graphs  will  be 
provided  in  Appendix  E  except  where  specific  references  are  made.  NSWC  data  sheets  which 
led  to  the  graphs  are  included  as  Appendix  C. 

5.1  Data  Reduction  and  Presentation 

The  following  equation  is  used  wholely  or  in  part  for  all  data  reduction  of  blast 
overpressures: 


log^P  =  +  B^og^D  +  B^A  +  B^Oog^Df 


(eq.  5.1) 


where: 


P  =  pressure,  psig  (or  1=  impulse,  psi-ms) 

D  =  distance  (slant  height  except  for  height  of  burst  measurements)  from 
vessel  center  (and  burst  point) 

A  =  absolute  value  of  angle  from  reference  blast  an^e  (reference  angle  found 
by  selecting  one  of  36  overpressure  curve  fits,  0®  to  180®  using  5® 
increments,  having  least  sum  square  errors) 

Bj  =  coefficients,  Bj  is  zero  unless  2nd  order  provides  an  accuracy  improvement 
(vessels  2-2  and  2-3  only) 


The  terms  involving  B„  and  B,  are  always  used.  This  provides  a  regression  of 
overpressure  or  impulse  (the  area  under  the  positive  portion  of  the  overpressure  curve)  versus 
distance  without  consideration  of  arena  angle.  The  addition  of  the  B2A  term  permits  a 
consideration  of  asymmetry.  Note  that  the  angle  is  taken  from  a  reference  angle  as  defined 
above.  The  second  order  term  in  logj)  was  used  only  for  8.7'  and  14'  height  of  burst  data 
reductions  where  a  non-linear  relationship  between  log^P  and  logJD  was  anticipated  from  the  test 
plan.  All  data  were  reduced  twice,  with  and  without  the  asymmetry  term. 

A  typical  curve  fit,  using  burst  of  vessel  P-1,  is  shown  in  Figure  5.1.  The  asymmetric 
curve  fit  which  calculates  an  angle  term  BjA  is  provided.  This  vessel  burst  had  a  calculated  peak 
overpressure  at  about  80®,  or  almost  normal  to  the  vessel.  The  same  best  fit  array  angle  is 
assigned  to  the  impulse  curve  fit  and  in  this  case  it  fits  very  well  as  can  be  seen  by  the  impulse 
error  deviation  being  nearly  the  same  as  that  for  pressure.  None  of  the  height  of  burst  data  were 
used  for  the  regression.  Pressures  were  measured  beyond  50'  but  these  data  are  also  not  used. 
Errors  are  provided  in  percent  of  reading  instead  of  full  scale  in  order  to  weigh  all  errors  equally. 
Data  shown  were  collected  at  ground  distances  of  10',  15',  22',  34'  and  50'.  The  curve  fit 
computes  slant  height  after  correcting  for  the  vessel  edge  rather  than  centerline  lying  on  the 
origin. 
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DATA  PROGRAM  VERSION  =  REGCKIO.FOR,  JAN  94 


CENTER  SPLIT 

VESSEL  P- 

-1  §  3250  PSIG 

(FILE=3248 

.RIO) 

3248 

TP: PRELIM,  TEST  DATE: 

6-26-90,  ORD=l 

,  OFFSET 

=1.,  NSWC  SHOT#: 

PHI= 

80.=BEST  FIT  CENTER  ARRAY  &  ANGLE  REFERENCE, 

HOB=  3. 

5  FT 

REGRESSION  MODEL:  LN (V)=B(N)  +B(1)  *LN(SL)+B(2)  *A+B  (3)  *  (LN  (SL)  )  **2+. 

AKHiJMA 

A 

D 

V(1) 

“PRESSURE 

Vf21=IMPULSE 

ANGLE 

ANGLE 

SLANT  MS 

MEAS 

CALC 

ERROR 

MEAS 

CALC 

ERROR 

DEG 

DEG 

FT  HT 

PSIG 

PSIG 

%RDG 

PSI-MS 

PSI-MS 

%RDG 

30. 

55. 

11.10  G 

9.90 

10.52 

6.3 

31.0 

36.7 

18.3 

60 . 

23. 

11.43  G 

9.00 

12.78 

42.0 

40.0 

58.8 

47.0 

75. 

6 . 

11.51  G 

14.30 

14.21 

-.6 

60.0 

75.1 

25.2 

90. 

10. 

11.54  G 

13.10 

13.82 

5.5 

60.0 

70,9 

18.1 

105 . 

26. 

11.51  G 

10.00 

12.36 

23.6 

50.0 

55.2 

10.4 

120. 

43. 

11.43  G 

13.60 

11.11 

-18.3 

38.0 

43.2 

13.6 

150. 

75. 

11.10  G 

7.50 

9.15 

22.0 

28.0 

26.9 

-3.9 

0. 

84. 

15.44  G 

6.70 

6.09 

-9.1 

18.0 

18.0 

-.  1 

30. 

53. 

15.91  G 

7.50 

7.31 

-2.6 

27.0 

28.2 

4.3 

60. 

22. 

16.26  G 

9.00 

8.90 

-1.1 

53.0 

45.0 

-15.2 

90. 

10. 

16.38  G 

10.10 

9.59 

-5.0 

64.0 

53.6 

-16.2 

0. 

83. 

22.30  G 

4.90 

4.18 

-14.6 

17.0 

13.7 

-19.6 

15 . 

67. 

22.55  G 

5.20 

4.60 

-11.6 

19.0 

17.1 

-10.0 

30. 

52. 

22.79  G 

6.50 

5.06 

-22.2 

20.0 

21.5 

7.4 

60  • 

21. 

23.14  G 

5.60 

6.18 

10.4 

38.0 

34.2 

-10.0 

98 . 

18. 

23.26  G 

6.80 

6.28 

-7.7 

58.0 

35.6 

-38.6 

105 . 

26. 

23.23  G 

5.70 

5.97 

4.8 

44.0 

31.9 

-27.6 

120 . 

41. 

23.14  G 

5.60 

5.38 

-4.0 

28.0 

25.1 

-10.3 

150 . 

72. 

22.79  G 

4.00 

4.40 

10.0 

19.0 

15.8 

-17.0 

0. 

82. 

34.19  G 

3.70 

2.70 

-27.1 

12.0 

9,9 

-17.9 

15 . 

67. 

34.45  G 

3.10 

2.97 

-4.1 

10.0 

12.4 

23.6 

30. 

51. 

34.69  G 

4.40 

3.28 

-25.4 

16.0 

15.5 

-2.9 

45. 

36. 

34.89  G 

3.60 

3.63 

.8 

19.0 

19.6 

3.1 

60 . 

21. 

35.04  G 

4.90 

4.02 

-17.9 

27.0 

24,7 

-8.4 

75 . 

5. 

35.14  G 

5.50 

4.47 

-18.8 

38.0 

31.3 

-17.6 

90. 

10. 

35.17  G 

5.00 

4.32 

-13.5 

37.0 

29.1 

-21.2 

0. 

81. 

50.13  G 

1.80 

1.82 

.9 

8.0 

7.3 

-8.5 

15. 

66. 

50.39  G 

1.50 

2.01 

33.8 

8.0 

9.2 

15.0 

30. 

51. 

50.63  G 

1.60 

2.22 

38.7 

9.0 

11.6 

28.6 

45 . 

36. 

50.83  G 

2.20 

2.46 

11.7 

11.0 

14.6 

32.6 

60  • 

21. 

50.99  G 

2.60 

2.73 

4.8 

16.0 

18.4 

15.0 

75. 

5. 

51.09  G 

2.80 

3.03 

8.1 

18.0 

23.3 

29.3 

98. 

18. 

51.11  G 

3.10 

2.77 

-10.8 

22.0 

19.1 

-13.3 

105 . 

25. 

51.09  G 

2.90 

2.63 

-9.2 

16.0 

17.1 

6.8 

120. 

41. 

50.99  G 

1.90 

2.37 

24.7 

13.0 

13.5 

4.0 

150. 

ERROR 

71. 

STD  DEV 

50.63  G  1.50 

IN  %  READING  OF 

1.93 
MEAS. ; 

28.6 

17.9 

7.0 

8.5 

21.4 

19 . 5 

PRESSURE  COEFS:  B(N)= 
B(l)  = 
B(2)  = 

5.24688 

-1.04286 

-.00700 

IMPULSE 

COEFS : 

C(N)=  6. 

C(l)= 
C(2)=  -. 

36686 

79778 

01544 

Figure  5.1,  Typical  Blast  Curve  Fit 
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Copies  of  the  regression  analyses  are  typically  not  furnished.  Plots  are  furnished  that 
together  with  the  furnished  error  deviation  provide  a  clear  picture  of  the  accuracy  and  distribution 
of  errors.  Table  of  coefficients  from  the  regression  analyses  may  be  found  in  Appendix  D. 
Table  D-1  provides  overpressure  curve  fit  coefficients  for  all  the  LSC  detonations.  Table  D-2 
and  D-4  provide  overpressure  coefficients  for  all  the  vessel  burst  data  and  Tables  D-3  and  D-5 
provide  impulse  coefficients  for  all  the  vessel  burst  data. 

Three  types  of  individual  burst  plots  will  be  used  for  overpressure  herein  and  will  be 
demonstrated  using  vessel  burst  P-1.  These  are  pressure  (or  impulse)  versus  angle,  pressure  (or 
impulse)  versus  distance  and  a  contour  plot. 

Pressure  versus  angle  and  impulse  versus  angle  plots  for  vessel  P-1  are  shown  in  Figure 
5.2  and  5.3  respectively.  The  data  points  are  connected  by  a  cubic  spline  for  clarity.  The  lines 
are  labeled  10'  through  50'  for  convenience.  Only  points  where  measurements  were  made  on  the 
ground  are  shown.  No  measurements  were  made  beyond  50'  for  this  test.  The  graphs  of  impulse 
versus  angle  are  only  furnished  for  a  limited  number  of  cases  such  as  vessel  bursts  5-3  and  5-4, 
Section  5.7,  where  overpressure  and  impulse  asymmetry  differ. 

The  dashed  line  shown  on  the  plots  were  computed  with  the  regression  coefficients.  The 
lines  are  truncated  to  the  maximum  angle  at  each  horizontal  distance  from  ground  zero  where 
measurements  were  made.  The  distances  shown  are  nominal,  ground  based  arena  measurements. 
The  computations  corrected  the  distances  shown  for  vessel  offset  and  slant  height.  This  would 
provide  a  one  to  one  graphical  correspondence  of  measured  to  computed  pressures  if  there  was 
a  perfect  fit  of  coefficients  to  all  data. 

Pressure  versus  distance  and  impulse  versus  distance  plots  for  vessel  P-1  are  shown  in 
Figures  5.4  and  5.5  respectively.  Two  regression  lines  are  shown.  The  solid  line  is  an  "average", 
obtained  without  an  angle  coefficient.  This  line  is  shown  solid  from  10'  to  50'  where  data  were 
typically  recorded,  and  dotted  at  closer  or  farther  distances.  A  dashed  line  above  (for  negative 
coefficient)  the  solid  line  and  labeled  "asymmetric"  in  the  legend  is  for  values  computed  with 
angle  "A"  (equation  5.1)  set  equal  to  zero,  i.e.  values  computed  for  the  best  fit  array  angle  which 
is  80°  for  vessel  P-1.  Two  dotted  lines  are  shown.  The  upper  line  represents  the  vessel  TNT 
equivalent  as  computed  using  ideal  gas  isentropic  expansion,  equation  4.1.  The  lower  TNT 
equivalence  line  represents  the  TNT  equivalence  of  the  LSC(s)  detonated  to  initiate  vessel  burst 
(see  Section  5.12). 

Measurements  made  at  the  vessel  height  of  burst  are  not  used  for  coefficient  computation. 
These  values  are  shown  in  Figures  5.4  and  5.5  with  open  circles.  Values  beyond  50'  distance  are 
not  shown  and  are  discussed  in  Section  5.8.  Measurements  made  on  the  ground  closer  than  10' 
are  shown  with  filled  circles.  These  were  included  in  coefficient  computation  if  the  data  were 
a  reasonable  extension  of  the  regression  line  with  scatter  typical  of  the  main  body  of  data.  Some 
measurements  were  a  considerable  departure  from  the  regression  line  as  shown  in  Section  5.5. 
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Overpressure,  psig  Overpressure,  psig 


Figure  5.2 

Overpressure  versus  angle  for  Vessel  P-1 


Figure  5.3 

Impulse  versus  Angle  for  Vessel  P-1 


Figure  5.4 

Overpressure  versus  Distance  for  Vessel  P-1 


Figure  5.5 

Impulse  versus  Distance  for  Vessel  P-1 
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Scoled  Dist.  (ft/lb“1/3) 


Arena  ongle,  deg. 


Figure  5,6,  Explosive  Blast  Overpressures  Figure  5.7,  Overpressures  front 

Detonation  of  High  Explosive  PHE-1 

The  overpressure  (impulse)  versus  angle  plots  graphically  present  the  asymmetry, 
irregularities  and  data  scatter.  The  overpressure  (impulse)  versus  distance  plots  permit 
comparison  of  blast  strength  versus  distance  and  permit  comparing  vessel  bursts  versus  each  other 
and  with  TNT  equivalence. 

A  special  regression  analysis  for  vessel  6A-4  is  discussed  in  Section  5.10.1. 

5.2  High  Explosive  Detonations 

High  explosives  were  detonated  as  discussed  in  Section  2.1.  Since  it  was  not  the  purpose 
of  the  test  program  to  evaluate  high  explosives,  these  are  only  discussed  in  Section  5.6  where 
vessel  HOB  tests  are  compared  to  high  explosive  tests.  However  data  plots  are  included  in 
Appendix  E  and  a  combined  plot  of  several  explosive  blasts  is  provided  as  Figure  5.6.  The  TNT 
equivalence  line  was  computed  using  equation  4.4  with  Z  =  D  (i.e.,  (TNT)*''^  =  1)  to  compare  all 
explosions  to  one  pound  of  TNT.  Table  5.1  compares  measured  data  in  Figure  5.6  to  equation 
4.4.  The  ratio  of  measured  overpressure  to  TNT  equivalence  should  be  the  reflection  factor  (see 
Section  4.4).  An  overpressure  angle  plot  is  provided  here  as  Figure  5.7.  The  figure  shows  the 
overpressure  resulting  from  detonating  33.3  pounds  of  composition  B  (TNT  equivalency  40.0  lbs). 
As  expected  the  asymmetry  in  the  high  explosive  data  is  minimal  compared  to  vessel  burst  data. 
This  is  evidenced  by  the  nearly  horizontal  regression  lines  in  Figure  5.7.  Considerable  data 
scatter  exist  in  the  results  from  this  detonation  as  it  does  in  all  the  data,  including  the  pressure 
vessel  burst  tests. 
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Table  5.1 

High  Explosive  Overpressures  versus  Equation  4.4 


Charge 

iJm, 

Equiv. 

Measured/Calculated  -  Eq.  4,4 

Ground  zero 
surface  condition 

Near  (10’) 

Far  (50’) 

65  lbs.  C-4 

94.0 

1.50 

1.82 

earth 

33.3  lbs.  Comp  B 

40.0 

1.48 

1.44 

earth 

9.9  lbs.  Pentolite 

14.6 

1.79 

2.23 

steel  plate 

4.5  lbs.  Pentolite 

6.6 

1.28 

1.50 

earth 

.66  lbs.  Pentolite 

.89 

2.58 

1.77 

steel  plate 

5.3  Linear  Shaped  Charge  Overpressures 

Linear  shaped  charge  overpressures  were  measured  without  attendant  vessel  bursts  during 
preliminary  testing,  Section  3.1.2.  LSCs  from  10  grains/foot  to  200  grains/foot  were  detonated. 
Coefficients  from  the  regression  analyses  are  provided  in  Table  D-1  in  Appendix  D. 
Overpressures  from  the  regression  analyses  are  plotted  in  Figure  5.8. 

The  purpose  in  making  the  measurements  was  to  quantify  the  LSC  overpressure  output 
to  get  an  idea  of  how  much  of  the  vessel  burst  output  might  have  to  be  attributed  to  the  initiating 
LSC.  Note  that  1/2  of  the  LSC  is  on  the  far  side  of  an  unburst  vessel.  It  was  also  desired  to 
subtract  the  LSC  output  from  the  vessel  burst  output  as  a  tare  if  that  was  practical.  (It  seems  not 
to  be.) 


The  LSC  overpressure  for  many  of  the  vessel  bursts  can  be  completely  ignored.  Delays 
were  observed  on  the  overpressure  versus  time  graphs  furnished  by  NSWC.  These  are  delays 
between  LSC  overpressure  peak  and  the  sudden  rise  to  the  vessel  peak  overpressure  at  burst.  The 
delay  is  caused  by  a  lower  than  ultimate  stress  level  after  LSC  detonation.  A  sample  delay  of 
9.8  ms  from  LSC  peak  to  vessel  burst  peak  is  shown  in  Figure  5.9.  The  delay  includes  a  gradual 
rise  in  pressure  after  initial  vessel  overpressure  rise.  The  delays  at  sample  angles  and  given 
distances  are  listed  in  Table  5.2.  For  those  delays  that  show  a  superscript  of  5  or  6  the  LSC 
overpressure  is  not  a  consideration  because  either  an  LSC  was  not  used  or  the  LSC  overpressure 
has  returned  to  zero  prior  to  the  vessel  burst  response.  Consideration  may  also  be  ruled  out  for 
superscript  of  1,  2  and  3  because  the  peak  vessel  pressure  did  not  occur  at  initial  vessel  pressure 
rise. 


Comparison  of  the  LSCs  in  Figure  5.8  that  were  actually  used  for  vessel  bursts  (25,  50 
and  200  gr/ft)  shows  that  measurements  were  35%  to  40%  more  than  the  LSC  TNT  equivalent 
at  10'  range  and  30%  to  60%  more  at  50'  range.  The  TNT  equivalence  is  computed  using  an 
equation  (equation  4.4)  from  Kinney  and  Graham’®  and  is  for  side-on  overpressrires. 


52 


Table  5.2 

Peak  Overpressure  Delays 
LSC  Peak  to  Vessel  Burst  Peak 
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Distance  (ft) 


Time  (ms) 


Figure  5.8,  Overpressures  from 
LSC  Detonation  on  2*  Diameter  Steel  vessel 


Figure  5.9,  Overpressure  vs  Time  for 
3250  psig  Vessel  P-1  at  9(T  and  15  ft 


5.4  Composite  Spheres 

Burst  overpressures  for  the  composite  sphere  vessels  are  shown  in  Figure  5.10  for  vessel 
PC  and  Figure  5.11  for  vessel  6A-2.  The  angle  plots  for  both  bursts  show  a  high  degree  of 
symmetry  about  the  origin  (regression  lines  almost  horizontal).  Vessel  PC  was  split  in  the 
horizontal  plane  with  1/2  of  the  vessel  flying  up  and  1/2  driven  to  the  ground.  The  reason  for 
the  dip  in  data  at  10'  distance  for  vessel  PC  are  unknown.  The  slanted  regression  lines  for  the 
vessel  6A-2  angle  plot  can  be  explained.  This  vessel  was  cantilever  supported  from  the  90°-180° 
side  with  the  split  in  the  vertical  plane  along  the  90“  axis.  The  90“-180“  side  fragment  accelerated 
more  slowly  than  the  0“-90“  side  because  it  carried  the  support  structure  with  it,  after  breaking 
it  off  at  the  base.  This  shielded  the  90“-180“  side  slightly  and  caused  the  highest  burst  pressures 
on  the  0“-90“  side. 

5.5  Burst  Pressure  Variation 

Vessel  burst  pressure  was  intentionally  varied  in  four  vessel  bursts  of  TP  #1 :  1-1  through 
1-4.  These  were  all  center-split  2'  diameter  steel  cylinders  burst  at  a  center  line  height  of  3.5 
feet.  In  addition  vessel  bursts  P-1,  P-2  and  2-1  were  conducted  at  3.5'  height  of  burst  (HOB)  and 
may  be  compared  to  TP  #1  vessels.  The  comparison  is  made  in  Figure  5.12  where  pressures  at 
10'  and  50'  nominal  ranges  are  computed  from  the  peak  value  with  symmetric  regression  analyses 
using  all  the  appropriate  bursts  and  then  are  plotted  in  the  figure.  The  dashed  lines  are  a  further 
regression  of  the  10'  and  50'  computed  points. 
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Figure  5.13  and  Figure  5.14  each  show  a  data  point  at  a  slant  height  of  3.5  feet  (the  vessel 
centerline  height).  These  bursts,  conducted  at  1475  psig  and  5425  psig,  respectively,  show  an 
overpressure  directly  under  the  vessel  that  is  quite  high,  higher  than  would  be  computed  with  the 
regression  coefficients  or  by  the  methods  of  Baker^  or  with  weak  shock  reflection  coefficient. 
The  pressure  rise  was  gradual  instead  of  a  step  rise  as  shown  in  Section  5.14.  (A  similar 
measurement  was  also  attempted  with  the  3450  psig  vessel  1-2  but  no  data  were  collected).  This 
shows  that  the  limited  flow  escape  area  rmder  the  vessel  causes  a  pressure  rise  which  would  not 
be  predicted. 

One  of  the  purposes  of  the  pressure  variation  Test  Plan,  TP  #1,  was  to  evaluate  real  gas 
effects.  However  data  scatter,  see  Figure  5.12,  was  too  great  for  an  evaluation.  The  overpressure 
caused  by  the  vessel  burst  is  comparable  to  its  TNT  equivalency  overpressure  at  50'  distance. 
Using  equation  4.4  with  a  50'  distance  and  a  17%  energy  decrease  (from  the  value  computed  with 
7125  psig  and  equation  4.1),  the  overpressure  decrease  is  1 1%.  The  overpressure  at  50'  distance 
from  the  three  vessels  burst  at  3.5'  HOB  and  3250  to  3475  psig  varied  more  than  1 1%  and  in  a 
random  fashion.  The  7125  psig  vessel  burst  overpressure  decrease  from  the  regression  lines  in 
Figure  5.12  must  include  real  gas  effects  but  of  unknown  magnitude. 

5.6  Height  of  Burst  Variation 

Three  steel  vessels  of  2  ft.  diameter,  center  split,  were  burst  at  centerline  heights  of 
3.5',  8.7'  and  14'.  (The  3.5'  HOB  was  included  in  Section  5.5.)  Preparations  for  these  tests 
included  special  pressure  probe  moimts  to  obtain  side  on  peak  overpressure,  prior  to  ground 
reflection  if  possible.  These  were  used  for  the  8.7'  and  14'  HOB.  The  vessels  were  burst  at  the 
same  nominal  pressure,  3450  to  3475  psig.  Data  from  both  8.7'  and  14'  HOB  were  combined 
and  are  shown  in  Table  5.3.  Three  probes  were  provided  imder  the  vessel  (90"  in  the  table)  and 
three  at  HOB  (75"  in  the  table).  The  closest  probes  at  locations  both  vmder  the  vessel  and  at 
height  of  burst  broke  off  due  to  side  loading  and  only  one  is  included.  A  regression  analysis  for 
the  pressure  values  yielded  the  coefficients  B^  =  4.36426  and  B,  =  -.95341. 

Two  regression  analyses  were  performed  for  all  the  groimd  based  data  at  all  three  HOBs: 
1.)  1st  order  in  logg  (distance)  with  no  angle  (or  asymmetry)  coefficient  and  2.)  1st  or  2nd  order 
with  an  angle  coefficient.  Both  are  used  for  pressure  versus  distance  plots  but  the  asymmetric 
coefficients  are  used  in  this  discussion.  First  order  was  used  for  3.5'  HOB  but  2nd  order  for 
8.7'  and  14'  HOB  because  a  non-linearity  with  logg  of  distance  was  anticipated  from  pre-test 
calculations.  The  data  in  Figure  5.15  resulted  from  division  of  the  pressures  versus  distance  at 
each  height  of  burst  by  the  probe  pressure  data.  Coefficients  were  used  to  compute  the  pressure 
values.  The  results  are  plotted  in  the  figure. 
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Overpressure,  psig  Peak  Overpressure,  psig 


Vessel  Internal  Pressure,  psig 


Figure  5.12,  Overpressure  Comparison 
from  Mid-Length  Split 


Figure  5.13,  Vessel  Burst  1-1  at  1475  psig 


Figure  5.14,  Vessel  Burst  1-3  at  5425  psig 


Figure  5.15,  Reflection  Factors  for  Vessel  Burst 
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Table  5.3 

Pressure  Probe  Measurements  for  Pressure  Vessel  Bursts 


Vessel 

Burst 

Angle, 

deg. 

Distance,  ft. 
Vessel  Center 
to  transducer 

Pressure, 

psig 

Imptilse, 

psi-ms 

2-2 

90 

1.83 

56.3 

Bad* 

2-2 

90 

4.25 

26.00 

35.0 

2-2 

90 

6.67 

6.50 

19.0 

2-2 

75 

15.97 

5.80 

38.0 

2-2 

75 

22.97 

3.70 

25.0 

2-3 

90 

7.50 

9.00 

17.0 

2-3 

90 

13.17 

8.00 

25.0 

2-3 

75 

15.97 

4.70 

32.0 

2-3 

75 

22.97 

6.20 

42.0 

*Probe  broke  after  initial  response  complete 


Three  high  explosive  (HE)  shots  of  4.5  lbs  pentolite  were  also  detonated  as  part  of  height  of 
burst  testing.  The  HE  shots  were  sized  to  yield  approximately  the  same  overpressure  at  the  10'  range 
as  the  pressure  vessel  burst.  The  HE  shots  for  8.7'  and  14'  were  mounted  from  the  vessel  support  and 
used  the  same  V  offset.  No  probes  were  used  under  the  shots  because  gain  changes  were  not  desired 
and  the  pressure  values  would  have  been  off  scale.  The  probe  data  that  were  collected  is  shown  in 
Table  5.4.  The  resulting  coefficients  from  regression  analyses  are  Bn  =  6.72653  and  B,  =  -1.81578. 

Reflection  factors  were  computed  in  a  manner  similar  to  that  for  the  pressure  vessels.  Dividing 
by  the  overpressure  for  TNT  equivalent  of  the  pentolite  charge  should  produce  reflection  factors  that 
compare  to  the  probe  based  factors.  The  comparison  is  poor.  The  factors  and  pressures  used  for 
computation  are  shown  in  Table  5.5.  Coefficients  for  14'  HOB  are  valid  for  17'  distance  (slant  height 
to  the  groimd)  and  beyond. 

The  poor  comparison  between  probe-based  and  TNT  equivalence  based  reflection  factors  for 
high  explosive  may  be  partly  caused  by  the  small  niunber  of  pressure  points  which  also  appear  to  be 
low  in  pressure.  This  also  makes  the  reflection  factor  for  pressure  vessels  suspect.  It  can  be  shown 
however  that  quantitatively  the  reflection  factors  for  pressure  vessels  do  change  with  distance  and 
pressure  in  a  manner  similar  to  that  for  high  explosives  but  not  necessarily  the  same  magnitude.  In 
Table  5.6  the  computed  pressures  at  8.7'  and  14'  HOB  for  both  pressure  vessels  and  high  explosives 
is  divided  by  the  3.5'  HOB  values  to  show  trends. 
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This  test  would  have  been  much  improved  if  many  more  height  of  burst  pressure  measurements 
were  made.  However  these  measurements  are  expensive  to  make  for  a  one  time  application  to  the 
support  structure  required. 


Table  5.4 

Pressure  Probe  Measurements  for  High  Explosive  Detonations 


HE 

angle 

Distance^  to 

Pressure, 

Impulse, 

shot 

:  transducers,  ft. 

psig 

psi-ms 

HE  2-2 

75” 

10.97 

10.60 

13.0 

HE  2-2 

750 

15.97 

5.30 

13.0 

HE  2-2 

75” 

22.97 

2.60 

9.0 

HE  2-3 

75" 

15.97 

5.80 

13.0 

HE  2-3 

75” 

22.97 

3.00 

6.0 

Table  5.5 

High  Explosive  Overpressures  and  Reflection  Factors 


Dist., 

ft. 

side-on  overpressure 

HOB,  ft 

pressure  @  : 
ground 

Reflection  factor 
based  on 

probe, 
(coefs.)  psig 

equiv.* 

(equa.),psig 

I»t)be. 

IKT 
equiv.  . 

17 

4.865 

7.30 

3.5 

13.109 

2.69 

1.79 

8.7 

16.256 

3.34 

2.23 

14 

20.572 

4.28 

2.82 

50 

.686 

1.26 

3.5 

2.021 

2.95 

1.60 

8.7 

2.684 

3.91 

2.13 

14 

3.754 

5.47 

2.98 

*  5.7  lbs.  used  based  on  average  probe  pressure  and  Figure  2F  of  Swisdak®. 
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Table  5.6 

Reflection  Factor  Trends  for  High 
Explosive  (HE)  and  Pressure  Vessels  (PV) 


Distance, 

ft 

Ratio  of  pressure  to  3.5'  HOB  test 

8,7' 

14' HOB  1 

HE 

PV 

HE 

PV 

17' 

1.24 

1.34 

1.57 

1.51 

50' 

1.33 

1.23 

1.86 

1.53 

5.7  Split  Location  Variation  with  2'  Diameter  Cylinders 

Prior  to  Test  Plan  5  several  center  split  bursts  were  conducted.  The  purpose  of  TP#5  was 
to  compare  results  using  a  different  split  location  with  the  same  2  ft.  diameter,  53  cubic  feet  steel 
cylinders.  Locations  used  for  the  split  were  3/4  -  1/4  length  and  full-endcap  length.  Two  bursts 
were  conducted  for  each  location.  For  each  burst  the  split  location  was  placed  at  the  arena  origin 
with  the  vessel  along  the  0“  -  180°  axis  and  the  short  fragment  end  toward  0°. 

The  3/4  -  1/4  bursts,  vessels  5-1  and  5-2,  provided  a  variation  in  fragment  velocity  but 
very  little  change  in  blast  results  from  center  split  vessels.  Blast  strength  was  reduced  slightly 
but  symmetry  was  approximately  the  same  as  center  split.  Table  D-4,  Appendix  D  shows  that 
all  the  center  split  vessels  plus  vessels  5-1  and  5-2  had  an  angle  of  blast  symmetry  within  15°  of 
normal  to  the  vessel  axis.  There  was  also  a  close  correlation  between  bmsts  of  vessels  5-1  and 
5-2.  Vessel  5-1  had  a  line  of  symmetry  of  75“  (using  overpressure)  with  a  deviation  of  16.9% 
of  reading  whereas  5-2  symmetry  angle  was  80°  with  a  deviation  of  16.6%. 

The  endcap  failures  had  a  greater  degree  of  asymmetry.  Again  both  bursts  correlated  well 
with  each  other,  vessel  5-3  having  an  angle  of  symmetry  of  55°  (using  overpressure)  with  a 
deviation  of  15.9%  versus  50°  and  15.8%  for  vessel  5-4.  The  impulse  data  also  correlated  well 
from  burst  of  5-3  to  burst  of  5-4.  The  impulse  regression  errors  are  high  due  to  apparent  data 
scatter,  however  the  appearance  of  impulse  versus  arena  angle  graphs  of  bursts  5-3  and  5-4  are 
similar  as  shown  in  Figure  5.16  and  5.17. 

A  strong  flow  along  the  0°  array  is  expected  from  the  endcap  vessels  and  high  impulse 
measurements  and  high  pressure  along  0°  might  also  be  expected.  However  the  distribution  of 
overpressure  versus  angle  does  not  correlate  with  impulse  versus  angle  for  either  endcap  vessel 
5-3  or  5-4.  This  can  be  seen  by  comparing  the  impulse  versus  angle  for  vessel  5-3,  Figure  5.16, 
with  the  overpressure  versus  angle  for  vessel  5-3,  Figure  5.18.  The  best  fit  angle  using  pressure 
is  55°;  but  the  best  fit  angle  using  impulse  appears  to  be  0°.  This  was  confirmed  with  a 
temporary  program  change  which  picked  the  best  fit  angle  due  to  impulse.  The  selected  angle 
was  0°  with  coefficients  C„,  Cj,  Cj  of  6.41162,  -.75905  and  -.01516  respectively.  However  the 
deviation  was  only  slightly  improved,  from  44.3%  to  37.4%. 
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Figure  5.16,  Impuke  for  Vessel  Burst  5-3 


Figure  5.17,  Impulse  for  Vessel  Burst  5-4 


Figure  5.18,  Overpressure  for  Vessel  Burst  5-3 


Figure  5.19  shows  the  overpressure  and  impulse  versus  time  at  15°  and  10  ft  distance. 
The  measurement  dipped  below  atmospheric  pressure  and  then  experienced  a  second  positive 
pressme  phase.  This  occurred  to  several  transducers  in  the  vicinity.  Most  of  the  waveform 
included  in  the  impulse  occurred  in  about  7  milliseconds.  But  at  15°  and  34  ft,  Figure  5.20,  the 
pressure  gradually  dropped  to  ambient  over  a  period  of  about  21  milliseconds.  This  seems  to 
explain  high  impulses  occurring  at  not  so  high  pressures. 

A  contour  plot  of  all  four  of  Test  Plan  5  vessels  is  included  as  Figure  5.21.  The  contours 
are  computed  from  the  coefficients  and  indicate  where  in  the  field  the  maximiim  overpressure  at 
50'  and  the  strongest  blast  angle  would  appear. 

5.8  Far  Field  Overpressure  Measurements 


Far  field  measurements  were  made  for  all  but  the  first  (preliminary)  and  last  (TP#6A)  test 
series  to  see  how  far  away  window  breakage  (approximately  1/2  psig  but  is  a  function  of 
thickness  and  size)  would  occur.  These  measurements  were  made  at  the  typical  HOB  of  3.5  ft 
due  to  uneven  terrain  beyond  the  50'  arena.  Measurements  used  side-on  pressure  probes  typically 
used  for  measurements  without  a  reflection  factor.  No  HOB  measurements  were  made  for  the 
last  test  series  due  to  expected  fragment  flight  toward  the  arena  transducer  field. 

Far  field  can  be  loosely  defined  as  the  distance  at  which  amplification  due  to  ground 
reflection  is  reduced  to  a  factor  of  approximately  one.  Baker^  provides  a  reflection  factor  of  1.1 

for  i?  =  1  and  1.0  for  R  greater  than  1 .  With  a  reflection  factor  of  1,  then  measurements 
made  on  the  ground  or  at  HOB  should  be  the  same.  The  distance  for  i?  =  1  can  be  found 
using  equation  5.2  (Baker’s  1  =  1  equation  solved  for  r). 


r  =  R 


(^-1)(2F^) 

_ 

Yl-l 


(eq.  5.2) 


where: 


Pi  =  initial  vessel  pressme 
=  atmospheric  pressmre 
2  is  for  burst  at  ground  level 
Fj.  =  initial  vessel  volume 

Yi  =  ratio  of  specific  heats  for  the  pressurant 

R  =  a  normalized  distance 
r  =  distance 
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Figure  5.19,  Pressure  and  Impulse  versus  Time  for  Vessel  Burst  5-3  at  15°  and  10  ft 
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Figure  5.20,  Pressure  and  Impulse  versus  Time  for  Vessel  Burst  5-3  at  15°  and  34  ft 
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Table  5.7 


Distance  to  R  bar  =  1 


Vessel 

Distance,  ft 

volume, 
cubic  feet 

pressure, 

psig 

53 

1475 

27.8 

53 

7125 

50.4 

2.7 

3975 

15.4 

Table  5.7  provides  distances  for  /?  =  1  for 
vessels  tested. 


Dist  (ft.)  along  0-180  deg.  axis 

Figure  5,21  y  Lines  of  equal  overpressure 
for  1/4  length  and  end  cap  vessels. 


The  trend  for  HOB  overpressure  and  impulse  is  to  be  less  than  the  average  regression  line 
at  10  ft.  distance  and  increase  as  the  far  field  is  approached.  At  distances  of  85  ft.  to  200  ft.  the 
TNT  equivalent  line  for  the  vessel  energy  sometimes  appears  to  be  a  better  predictor  than  the 
average  regression  line  as  shown  in  Figure  5.22  for  overpressure  of  vessel  5-1  and  in  Figure  5.23 
for  impulse  from  vessel  1-3.  However,  the  overpressure  for  the  composite  spheres  is  above  the 
TOT  equivalence  line  for  much  of  the  distance.  See  Figure  5.10  for  vessel  PC  and  Appendix  E, 
Figure  E24  for  vessel  6A-2.  Reflection  does  not  account  for  the  higher  pressures  since  the 
reflection  factor  should  be  approximately  one  for  15  to  50  feet.  (See  also  Section  5.15). 

5.9  Diameter  Variation  with  Center  Split 

Three  2'  diameter  center  split  vessels  were  burst  at  a  pressure  range  of  3250  to  3450  psig 
at  3.5'  HOB:  P-1, 1-2  and  2-1,  P-1  will  be  used  to  compare  to  the  center  split  3'  diameter  vessel, 
6A-3,  because  its  burst  pressure  of  3250  most  closely  matches  the  6A-3  burst  pressure  of  3300 
psig.  Based  on  the  inner  diameters,  the  3'  vessel  has  2.1  times  the  flow  area  of  the  2'  vessel  and 
was  expected  to  develop  higher  blast  overpressures.  Since  impulse  should  be  a  function  of  total 
energy  rather  than  release  rate,  a  higher  impulse  for  the  3'  diameter  vessel  was  not  expected.  Ehie 
to  scatter  it  is  difficult  to  compare  overpressures  in  Figures  5.2  and  5.24  or  impulses  in  Figures 
5.3  and  5.25.  Using  symmetric  coefficients  the  overpressure  and  impulse  for  both  vessels  is 
compared  at  distances  of  10’  and  50’  in  Table  5.8.  Based  on  the  table  it  is  unclear  if 
expectations  were  realized. 


Table  5.8 

Comparison  of  Vessels  P-1  and  6A-3 


Distance, 

ft. 

Pressure,  psig 

Impulse,  psi-ms 

P-1 

6A-3 

P-1 

6A-3 

10’ 

13.079 

14.082 

50.6 

37.2 

50’ 

3.22 

2.39 

19.45 

19.52 
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5.10  Dual  versus  Multi-fragment  in  2'  Cylinders 

The  multi-fragment  vessel  was  expected  to  produce  a  high  blast  overpressure  due  to  the 
sudden  release  of  energy,  however  the  distribution  of  overpressure  with  respect  to  geometry  was 
uncertain.  For  center  split  vessels  the  distance  to  the  vessel  from  the  semi-circular  field  of 
overpressure  transducers  may  be  taken  as  the  distance  to  the  split.  An  overpressure  plot  with 
symmetry  about  the  origin  might  be  expected.  However  the  exhausts  from  each  section  impinge 
at  the  center  and  a  strong  blast  field  normal  (or  nearly  normal)  to  the  vessel  axis  is  produced. 
Typical  plots  are  shown  in  Figure  5.26.  The  closest  distance  to  choose  for  the  multi-fragment 
burst,  with  the  sudden  disintegration  of  the  whole  vessel,  is  not  so  clear.  Figure  5.26  indicates 
that  the  end  points  of  the  10  ft  and  15  ft  semicircles  experienced  high  pressures  due  to  proximity 
to  the  ends  of  the  vessel.  This  trend  is  even  displayed  to  a  lesser  extent  at  the  34  ft  distance. 
There  is  also  a  tendency  for  a  strong  blast  field  normal  to  the  vessel  centerline,  showing  that  the 
confined  jet  impingement  as  in  the  center-split  is  not  necessary  to  cause  that  type  of  asymmetry, 
A  special  regression  analysis  of  test  data  for  a  hazard  assessment  workbook''®  is  repeated  as 
Section  5.10.1. 

Overpressures  and  impulses  were  also  considerable  higher  for  the  3500  psi  multi-fragment 
vessel  than  any  other  due  to  the  more  sudden  release  of  energy.  This  is  even  true  when 
considering  the  7125  psi  center  split  vessel. 

5.10.1  Special  Regression  Analysis,  Vessel  6A-4 

Figure  5.26  contains  an  angle  plot  of  overpressure  from  the  multi-fragment  vessel,  6A-4. 
The  distances  are  measured  radially  from  the  center  of  the  arena,  and  center  of  the  vessel.  (The 
distances  are  all  nominal  and  refer  to  distances  along  the  ground.  Data  reduction  used  slant 
height  distances.)  Data  at  10’,  15°;  15’,  0°  and  15’,  180°  (points  labeled  1,  2  and  3)  showed  a 
trend  in  opposition  to  the  balance  of  the  arena  which  showed  high  overpressures  near  90°  and  low 
near  0°  and  180°.  For  this  reason  the  distance  (slant  height)  for  all  the  points  was  computed  at 
90°  to  the  vessel  axis.  The  data  are  replotted  in  Figure  5.27. 

In  Figure  5.27  the  three  points  are  again  labeled  and  it  is  seen  that  these  points  are  not 
high  but  are  either  average  (point  1)  or  quite  low  (points  2  and  3)  compared  to  the  regression 
line,  "average  90°  distance".  Thus  it  is  seen  that  the  overpressure  is  high  along  a  line  normal  to 
the  center  of  the  vessel  and  is  low  at  more  remote  distances  from  the  normal  line.  (See  also 
Appendix  F.) 

Four  regression  lines  are  shown  in  Figure  5.27.  The  two  dashed  lines,  "  maximum  radial 
distance"  and  "average  radial  distance"  are  the  average  and  maximum  pressure  from  the  type  of 
data  reduction  described  in  Section  5.1,  with  radial  distances.  The  line  "average  90°  distance" 
results  from  a  regression  of  the  Figure  5.27  data  with  the  solid  symbols  excluded.  These  were 
beyond  the  vessel  length  and  appeared  to  be  sufficiently  low  for  both  pressure  and  impulse  to 
discard.  The  "maximum  90°  distance"  line  results  from  a  regression  using  points  at  or  above 
"average  90°  distance"  line  and  replaces  "maximum  radial  distance".  The  effect  of  the  new 
regression  lines  is  to  lower  the  "measured"  pressure  for  vessel  6A-4  for  close  distances. 
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Figure  5.27,  Overpressure  of  Vessel  6A-4 


Figure  5.28,  Impulse  of  Vessel  6A-4 


The  overall  appearance  of  the  impulse  data  is  very  similar  to  the  pressure  data  and  is 
shown  in  Figure  5.28.  Coefficients  based  on  the  special  regression  analyses  are  shown  on  the  last 
line  of  Tables  D-2  through  D-5  of  Appendix  D. 

5.11  Dual  versus  Multi-fragment  in  3'  Cylinders 

Vessel  6A-1  was  expected  to  burst  along  horizontal  grooves  as  well  as  at  the  two  end  cap 
grooves.  This  did  not  happen  and  a  double  end  cap  failure  resulted  with  the  center  section 
remaining  in  place  (shifted  a  foot). 

The  double  end  cap  failure  results  in  twice  the  pressurant  exhaust  area  than  the  center  split 
and  can  be  expected  to  result  in  higher  overpressures.  This  did  occur  as  shown  in  Figure  5.24 
versus  Figure  5.29.  In  Figure  5.29  three  peaks  can  be  seen  in  the  overpressure  pattern  at  10  ft 
distance.  The  side  peaks  represent  the  outrush  of  gas  at  each  end  and  the  center  peak  is  evidently 
due  to  the  intersection  of  the  flow  from  each  end.  The  high  pressure  tendency  at  the  center  of 
the  field  remains  out  to  the  50'  distance.  This  burst  also  exhibited  a  greater  impulse  field  than 
the  center  split  3'  diameter  or  the  2'  diameter  center  split  vessels  at  any  pressure  range. 

5.12  Bursts  compared  to  TNT  equivalence 

TNT  Equivalence  lines  in  the  pressure  (or  impulse)  versus  distance  plots  were  identified 
in  Section  5.1.  The  TNT  equivalence  for  the  upper  line,  vessel  TNT  equivalence,  was  computed 
using  the  isentropic  expansion  energy  of  the  vessel,  with  a  factor  of  1.545  x  10®  applied  to 
convert  from  foot  pounds  of  energy  to  lbs  TNT.  The  conversion  factor  is  within  a  range  of  1.4 
to  1.6  presented  by  various  authors*.  The  factor  selected  is  from  Kinney  and  Graham*®.  No 
corrections  are  made  for  energy  expended  in  propelling  fragments  rather  than  creating  a  blast 
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wave.  Once  the  TNT  equivalence  is  known,  the 
overpressure  (or  impulse)  is  computed  using 
equation  4.4  (or  4.5  for  impulse)  from  Kinney 
and  Graham’®.  The  overpressure  (or  impulse) 
equations  are  for  side-on  overpressures  and  do 
not  consider  reflection  (amplification)  factors 
other  than  1.0.  No  corrections  are  applied  for 
ground  reflection  factors.  Some  experimentation 
was  conducted,  Section  5.6,  but  insufficient  data 
is  available  for  quantifying  vessel  burst  reflection 
factors. 

A  constant  TNT  equivalence  is  used  in 
the  high  explosive  plots  and  was  taken  from  the 
mid  range  (34  ft)  distance  overpressure  using 
Swisdak’s’  report.  The  report  shows  that  the 
conversion  factor  for  pounds  of  an  arbitrary 
explosive  to  pounds  of  TNT  is  not  a  constant. 

Most  vessel  bursts  used  a  linear  shaped  charge,  typically  in  a  machined  groove,  to  initiate 
a  vessel  burst.  Overpressures  due  to  LSCs  detonated  on  a  steel  vessel  surface  are  presented  in 
Section  5.3.  These  are  not  used  for  the  LSC  overpressure  (or  impulse)  versus  distance  lines  in 
the  plots.  Rather  the  conversion  factor  for  RDX  to  TNT  was  obtained  from  Kinney  and 
Graiham’®  and  the  values  for  the  lines  determined  as  discussed  above.  Neither  the  TNT  equivalent 
overpressure  for  the  initiating  LSC  from  reference  10  nor  the  overpressure  from  Figure  5.8  is 
added  to  the  vessel  TNT  equivalence  to  compare  to  actual  data  in  the  data  overpressure  versus 
distance  plots.  This  might  be  convenient  to  have,  but  how  and  whether  to  add  it  is  unclear. 
Timing  eliminates  some  LSCs  from  consideration.  For  others  the  rated  conversion  to  TNT  from 
reference  10  may  not  be  a  constant  but  may  vary  with  overpressure  per  reference  9.  Additionally 
reflection  factor  is  not  a  constant  but  varies  with  distance. 

From  Table  5.2,  the  contribution  of  LSC  overpressure  for  vessel  burst  overpressure  can 
probably  be  ignored  for  the  following  bursts:  P-1,  P-2,  1-1,  1-2,  1-3,  2-1,  2-1,  2-3,  5-1,  6A-3 
and  possibly  5-3. 

5.13  All  3.5  ft  HOB  bursts  compared 


Areno  ongle,  deg. 


Figure  5.29,  Overpressure from  Vessel  6A-1 


It  is  convenient  to  borrow  the  Hopkinson  or  "cube  root  scaling"  law  for  high  explosives^® 
and  apply  it  to  vessel  bursts.  This  requires  dividing  the  actual  distance  by  the  cube  root  of  the 
vessel  TNT  equivalence.  In  this  way  large  and  small  vessel  bursts  (or  explosions)  may  be 
compared  to  each  other  and  to  die  overpressure  (which  is  not  scaled)  for  one  pound  TNT. 
Similarly  if  the  impulse  is  divided  by  the  cube  root  of  the  vessel  TNT  equivalence,  then  the 
impulse  for  large  and  small  vessel  bursts  may  be  compared. 
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Figure  5.30,  Overpressure  Comparison  Figure  5.31,  Impulse  Comparison 

at  10’,  50’ for  all  3.5’  HOB  Bursts  at  10’,  50’ for  all  3.5’  HOB  bursts 

Figure  5.30  is  a  plot  of  peak  overpressure  versus  scaled  distance  for  all  vessels  burst  at 
a  centerline  height  of  3.5  ft.  All  the  points  are  evaluations  from  coefficients  (no  angle  term)  at 
either  10  ft  or  50  ft  distance  and  are  connected  by  a  dashed  line.  The  distances  represent  distance 
extremes  for  the  regression  analysis  (actually  10.6  or  11.5  to  50.1  or  51.1  when  corrections  for 
offset  and  slant  height  are  applied).  The  open  circles  are  all  replotted  from  Figure  5.12,  of 
overpressure  versus  vessel  pressure  (i.e.,  coefficients  used  to  compute  overpressure  at  10'  and  50' 
distances).  A  regression  line  is  plotted  for  the  open  circled  points  (cylinders  burst  at  mid-length). 
The  two  points  considerably  below  the  line  are  for  vessel  1-1.  No  explanation  for  the  low 
pressures  can  be  offered.  The  filled  points,  all  at  or  above  the  regression  line,  have  a  vessel 
and/or  burst  geometry  which  permits  a  faster  venting  of  the  pressurized  gas.  The  open  squares 
and  triangles  have  a  burst  geometry  which  increases  the  exhaust  time. 

Also  plotted  in  Figure  5.30  is  an  equation  showing  the  overpressure  versus  scaled  distance 
for  the  TNT  equivalence  of  the  vessel  (or  for  overpressme  versus  distance  for  a  one  pound  TNT 
explosion).  It  is  apparent  that  increasing  the  exhaust  rate  increases  the  overpressure  towards  the 
TNT  equivalence  line.  Some  of  the  points  appear  to  be  greater  than  the  TNT  equivalence 
pressure.  This  can  be  explained  (except  for  the  spheres  at  15’  and  50’)  by  the  fact  that  the  TNT 
pressures  are  incident  pressures,  not  amplified  by  ground  reflections,  whereas  the  vessel 
overpressures  are  reflected  pressures.  See  height  of  burst  variation.  Section  5.6. 

Regression  coefficients  were  also  obtained  for  log^  (impulse)  versus  log^  (distance)  with 
no  angle  term.  As  before  the  coefficients  were  used  to  compute  the  impulse  at  the  10'  and  50' 
range.  The  distances  were  then  converted  to  scaled  distance  by  dividing  by  the  TNT  equivalence 
of  the  vessel  energy.  The  data  is  plotted  in  Figure  5.3 1  where  again  the  points  are  connected  by 
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a  dashed  line  and  a  regression  line  is  found  for  the  mid-length  split  vessels  with  2'  diameter.  The 
TNT  equivalence  impulse  equation  is  also  plotted.  The  ordinates  are  now  scaled  values  as  are 
the  abscissas. 

The  measured  impulse  values  tend  to  be  close  to  the  TNT  equivalence  for  any  distance. 
The  fact  that  some  points  appear  to  be  greater  than  TNT  equivalence  can  again  be  explained  by 
the  fact  that  the  TNT  equivalence  uses  incident  overpressure  and  the  measured  values  are 
reflected. 


5.14  Pressure  versus  Time  Waveform 

Ideal  high  explosive  waveform  characteristics  include  a  sharp  rise  followed  by  an 
exponential  decay.  Two  of  the  traces  in  Figure  5.32  are  similar  to  high  explosive  response.  All 
recordings  shown  were  made  from  transducers  located  in  the  same  general  location:  at  a  range 
of  10  feet  to  15  feet  and  within  30“  of  normal  to  center  of  the  long  vessel  axis. 


Figure  5.32  presents  data  from  vessel 
bursts  P-1,  6A-2  and  6A-4.  P-1  is  a  cylinder 
burst  about  its  mid-length.  The  square  wave 
(also  shown  in  Figure  5.9  in  Section  5.3)  was 
seen  at  distances  of  22  feet  and  closer.  At 
greater  distances  the  waveform  transitioned  to 
an  exponential  decay.  The  pre-burst  pressure 
rise  is  the  LSC  detonation.  Vessel  6A-2  is 
the  latter  composite  sphere.  The  equatorial 
split  presented  a  large  ratio  of  exhaust  area  to 
vessel  volume  compared  to  the  cylinder. 

Additionally  the  lightweight  fragments 
accelerate  rapidly  and  minimally  reduce  the 
vent  area  at  the  beginning  of  laimch.  The 
sphere  waveform  has  the  sharp  peak,  semi¬ 
exponential  decay  and  second  shock  similar  to 
a  high  explosive  blast.  Vessel  6A-4  was  burst 

into  14  fragments,  providing  a  large  vent  area  for  a  cylindrical  burst.  The  waveform  has  a  sharp 
peak,  a  generally,  although  somewhat  ragged,  exponential  decay  and  second  shock.  The 
overpressme  measured  was  greater  than  with  any  other  cylindrical  bturst. 

Figure  5.33  shows  the  13  millisecond  delay  between  LSC  detonation  and  burst  of  vessel 
5-1  as  mentioned  in  Section  3.1.1. 


Figure  5.32,  Vessel  Burst  Overpressure  Time  Response 


Figures  5.34  and  5.35  show  the  pressure  on  the  ground  directly  under  the  burst  of  vessels 
1-1  and  1-3  as  mentioned  in  Section  5.5.  In  each  case  the  transducer  experienced  a  gradual 
pressure  rise  rather  than  a  step  rise  as  expected  for  a  shock  passing.  This  characteristic  was  not 
observed  anywhere  else  in  the  arena. 
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Figure  5.55,  Burst  Delayed  13  ms  after  LSC  Detonation,  Vessel  5^1 
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Figure  5.34,  Pressure  Measured  on  the  Ground  Under  Vessel  1-1 


Figure  5.35,  Pressure  Measured  on  the  Ground  Under  Vessel  1-3 


2 


Distance,  ft. 


Figure  536 

Comparison  of  Composite  Sphere  Overpressures 


Figure  537 

Comparison  of  Multi  fragment  Cylinder  Overpressures 


5.15  Selected  Data  Compared  to  Literature 

Curves  published  by  Baker^  assxime  sudden  vessel  wall  disappearance.  Two  types  of 
bursts  conducted  under  the  Burst  Study  approach  that  assumption  and  can  be  compared.  One  is 
the  composite  spheres  which  have  a  large  vent  area  to  volume  ratio  and  fragments  which  rapidly 
accelerate  to  high  velocity.  The  other  is  the  multi-fragment  vessel  that,  with  14  fragments,  also 
has  a  very  large  vent  area  to  volume  ratio.  Figure  5.36  presents  data  for  the  initial  composite 
sphere  burst,  vessel  PC.  The  measured  data  was  recorded  at  ground  level  and  hence  is  reflected 
data.  The  calculated  line  using  "volume  correction"  uses  a  volume  of  twice  actual  as  a  partial 
reflection  correction.  This  decreases  the  distance  by  20%  and  increases  the  pressure  by  about 
40%  in  the  near  field  and  90%  in  the  far  field.  The  curve  having  additional  ground  corrections 
includes  a  factor  of  2  times  the  "volume"  curve  at  near  field  and  1 .1  at  far  field.  (Calcvilated  data 
procedure  includes  picking  points  from  curves  of  normalized  values,  hence  the  lack  of  smooth 
curves.) 

Figure  5.37  presents  data  for  cylinder  burst  6A-4.  This  cylinder  broke  into  14  fragments 
of  approximately  equal  weight,  which  also  approximates  the  disappearing  wall  assumption.  The 
figure  presents  curves  for  measured,  calculated  with  volume  correction  and  a  curve  which 
contains  both  cylinder  correction  factors  and  ground  reflection  factors.  The  cylinder  correction 
factors  (which  Baker  et  al  stated  are  "very  crude")  varied  from  4.5  in  the  near  field  to  1.4  in  the 
far  field.  The  reflection  correction  factors  varied  from  2.0  in  the  near  field  to  1 . 1  in  the  far  field. 

It  is  doubtful  that  much  greater  overpressures  from  these  two  vessels  would  ever  occur 
accidentally  since  this  would  require  a  greater  number  of  fragments  with  attendant  faster  venting. 
The  corrections  to  the  cylinder  calculations  which  are  required  by  the  reference  appear  to  be 
excessive  except  at  the  50  ft.  distance. 
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6.0  FRAGMENT  MODELING  AND  TEST  RESULTS 


The  explosive  failure  of  a  pressure  vessel  not  only  may  generate  a  blast  wave  but  may 
also  propel  fragments.  The  fragments  may  be  massive  and  the  velocities  quite  high,  presenting 
a  significant  hazard.  The  lethality  and  damage  potential  of  a  fragment  propelled  by  a  bursting 
pressure  vessel  depends  on  the  impact  velocity  and  mass  of  the  fragment.  Both  the  impact 
velocity  and  distance  a  fragment  may  travel  depend  on  the  initial  velocity.  The  Burst  Test  Study 
used  a  two  pronged  approach  to  the  study  of  fragments.  The  initial  velocity  of  fragments  was 
measured  and  a  computer  program  for  the  prediction  of  initial  velocities  was  rewritten  and 
improved.  Maximum  distance  travelled  by  fragments  was  also  measured  and  for  initial  tests  the 
ricochet  points  were  also  obtained. 

6.1  Measurement  of  Initial  Velocities 

This  section  is  concerned  with  fragment  initial  velocity.  Initial  velocity  will  mean  the 
maximum  velocity  the  fragment  achieves  after  vessel  breakup  prior  to  a  decay  in  velocity  due  to 
aerodynamic  drag.  During  the  acceleration  phase  the  fragment  is  acted  upon  by  the  expanding 
containment  gas  and  by  drag.  Initial  velocity  is  the  resulting  velocity  at  the  end  of  that  phase. 
Drag  can  almost  be  ignored  in  the  calculation  of  initial  velocity  of  fragments  from  high  pressure 
vessels,  but  greater  accuracy  is  obtained  by  considering  drag. 

6.1.1  Methods 

Two  methods  were  used  to  capture  fragment  displacement  versus  time.  The  primary 
method  was  high  speed  photography  or  in  some  instances  high  speed  video  or  both.  For  instance 
only  high  speed  video  was  available  for  the  burst  of  vessel  P-2  since  the  time  of  burst  was 
unknown:  the  shaped  charge  was  detonated  the  day  before  in  a  failed  test.  The  running  time  of 
the  high  speed  cameras  is  dependent  on  the  frame  speed  in  use  but  is  typically  less  than  15 
seconds.  The  second  method  was  the  use  of  contact  wires.  Sometimes  the  wires  were  contact 
making  and  sometimes  contact  breaking. 

High  speed  photography  was  the  preferred  method  and  was  useable  for  all  of  the 
fragments  for  most  of  the  tests.  It  has  the  disadvantage  that  the  fragment  may  be  obscmed  from 
the  camera  for  a  period  of  time  by  the  cloud  of  dust  and  condensation  caused  by  the  vessel  biorst. 
For  most  fragments  this  was  not  a  problem,  the  fragments  would  be  visible  for  a  sufficient  period 
of  time  to  track  them  for  several  distance  measurements.  However  the  endcaps  from  vessels  5-3 
and  5-4  were  never  seen  with  the  high  speed  cameras.  The  field  of  view  was  limited  to  get 
adequate  resolution  for  a  high  velocity  fragment.  They  were  seen  with  standard  video  cameras 
having  a  wider  field  of  view,  however  the  camera  angles  were  not  good  and  the  frame  rate  was 
slow. 


Contact  wires  were  used  for  backup  and  as  an  attempt  to  obtain  acceleration  data 
immediately  after  rupture  when  the  vessel  may  be  obscured  by  the  dust/condensation  cloud. 
Contact  break  wires  were  initially  used.  These  are  preferred  by  the  NSWC  instrument  personnel 
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when  brittle  wires,  such  as  pencil  leads  that  break  cleanly,  can  be  used.  However  since  the  exact 
fragment  launch  direction  was  unknown,  long  wires  were  used  that  would  be  broken  by  the 
leading  surface  of  the  fragment,  hopefully  with  the  wire  framework  imdamaged.  This  caused 
some  concern  about  accuracy  because  the  copper  wires  in  use  (Test  Plans  1  and  2)  could  stretch 
before  breaking.  Additionally  the  possibility  existed  for  shock  loading  to  deflect  the  wires  and 
effect  the  data. 

For  vessel  PC  a  new  contact  wire  arrangement  was  used  that  utilized  pretensioned  piano 
wire  in  a  more  rigid  frame  and  a  change  of  the  signal  sense  from  contact  breaking  to  contact 
making.  Crossing  grids  of  wires  one  inch  apart  was  used  to  ensure  continuity  in  the  face  of  high 
acceleration  loads  (a  common  wire  was  previously  attached  to  the  vessel)  and  to  allow  for  a 
slightly  angled  fragment  motion.  Six  sets  of  grids  were  used  at  one  foot  intervals.  For  the  next 
vessel  test,  5-3,  the  makewire  contact  sensing  was  changed  to  breakwire  due  to  readability 
problems  and  the  frame  was  turned  on  its  side  .  Damage  during  testing  and  modifications  to  the 
stand  are  discussed  in  Section  3.7. 

Distances  to  contact  wires  for  three  vessels  are  shown  in  Table  6.1.  For  these  vessels;  PC, 
5-3  and  5-4;  the  fragment  velocities  are  listed  in  Table  6.2  using  the  wire  data  rather  than 
photography  data. 


Table  6.1 

Distances  to  wires,  ft. 


Vessel  I 

Wire# 

PC’ 

5-3^ 

5-4^ 

1 

.68 

.69 

.70 

2 

1.68 

2.07 

1.65 

3 

2.68 

3.40 

2.66 

4 

3.68 

5.10 

3.64 

5 

4.68 

— 

4.63 

6 

5.68 

— 

5.64 

’Makewire 

^Breakwire 
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Figure  6.1,  Typical  Type  of  Figure  6.2,  Regression  Line  Omitting 

Distance/Velocity  Plots  from  NSWC  Acceleration  Distance  Points 

6.1.2  Data  Reduction 

Typically  each  camera  was  placed  normal  to  the  expected  fragment  track  and  aimed  at  a 
location  where  the  fragment  would  be  after  motion  commenced.  Frequently  the  unburst  vessel 
was  at  the  edge  of  the  field  of  view.  A  series  of  2"  by  4"  lumber  placed  slightly  off  the  fragment 
track  provided  data  for  a  length  calibration  and  for  angle  corrections  when  the  fragment  was  not 
exactly  normal  to  the  camera  view.  High  speed  motion  pictures  were  digitized  to  distance  versus 
time  data  by  NSWC  using  a  special  projection/digitizing  table  and  computer  system.  Light  bursts 
on  the  edge  of  the  film  provided  the  actual  frame  rate  which  was  used  in  place  of  the  nominal 
value.  Spread-sheet  software  was  then  used  for  calculations. 

Data  furnished  by  NSWC  included  point  to  point  velocity  versus  time  as  well  as  distance 
versus  time.  Figure  6.1,  plotted  by  GP,  shows  the  type  of  typical  velocity  data  and  plots 
furnished  by  NSWC. 

GP  then  used  the  distance  versus  time  data  and  a  computer  program  to  obtain  the  velocity. 
The  program  provided  a  slope-intercept  regression  analysis,  a  slope  printout  from  the  intercept 
to  each  point  and  the  option  to  exclude  points  from  the  regression  analysis.  The  data  printout 
made  it  obvious  which  points  were  part  of  the  initial  acceleration  and  which  points  were  part  of 
the  deceleration  due  to  drag.  By  discarding  only  these  points  an  accurate  and  maximum  velocity 
was  obtained.  Usually  two  to  four  trips  through  the  routine  yielded  a  final  value.  Figure  6.2  is 
a  plot  of  the  data  and  the  regression  line  for  the  east  fragment  from  vessel  5-4. 
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6.1.3  Accuracy  of  photo-based  data 


Accuracy  in  obtaining  velocities  from  high  speed  photography  probably  varied  with  each 
vessel  bmst.  The  accuracy  achieved  with  vessel  5-4  was  examined  and  is  probably  typical  for 
fragments  that  travelled  in  a  predictable  direction,  normal  to  the  camera  field  of  view  as  noted 
in  Table  6.2.  Figure  6.1  shows  that  some  data  scatter  exists,  resulting  from  inaccuracies  in 
reading  the  vessel  position.  However  65  points  were  used  in  the  regression  analysis  for  slope 
(velocity)  and  intercept,  providing  an  error  variance  for  computed  distance  of  a  small  fraction  of 
an  inch.  There  is  still  a  potential  calibration  error  of  scale  distance  versus  actual  distance.  This 
calibration  scale  was  established  using  2"  X  4"  lumber  at  known  distances  that  could  be 
accmately  measured,  and  this  error  is  also  expected  to  be  small,  less  than  0.5% 

In  spite  of  the  inherent  high  accuracy  to  be  expected  from  photo  data  recorded  in  the 
manner  described,  certain  things  can  affect  the  data  and  render  a  result  outside  of  expectations. 
These  cases  are: 

•  Due  to  a  high  pressure  area  between  supports,  fragments  of  vessel  6A-3  may  have 
been  pivoting  during  acceleration,  resulting  in  a  lower  actual  velocity  than  might 
otherwise  have  been  recorded. 

•  One  sidewall  fragment  from  vessel  burst  6A-4  was  tracked  at  500  ft/second. 
Others  that  were  tracked  (over  half  of  the  12  could  not  be  seen)  at  considerably 
lower  velocities  used  cameras  with  less  than  optimum  angles. 

•  The  lightweight  spherical  vessels  PC  and  6A-2  were  composite  overwrapped 
pressure  vessels  having  stainless  steel  liners  and  Kevlar  overwrap.  It  is  believed 
that  the  overwrap  was  tracked  for  vessel  fragments  PC  top  and  6A-2  west, 
resulting  in  a  lower  velocity  than  would  otherwise  have  been  measured.  These 
two  fragments  will  be  discussed  in  Section  6.1.4. 

6.1.4  Lightweight  sphere  fragment  velocities 

Two  lightweight  spherical  vessels  were  burst,  resulting  in  very  high  fragment  velocities. 
Considering  both  vessel  bursts  together,  the  correct  velocity  for  fragments  from  both  vessels  may 
be  between  978  Q)s  and  1252  ^s.  The  uncertainty  lies  in  the  wire  firework  used  for  vessel  PC 
and  in  the  condition  of  overwrap  shedding  for  both  vessels.  The  Kevlar/epoxy  overwrap  from 
these  vessels  appeared  to  be  shedding  at  the  start  of  photo  tracking.  The  overwrap  had  been 
painted  white  and  being  light  in  color  and  having  the  tendency  to  fluff  out,  it  was  not  only  easier 
than  the  stainless  steel  liner  to  track  but  it  also  decelerated  faster.  This  can  be  seen  in  the  high 
speed  film  for  vessel  PC.  It  is  believed  to  have  also  occurred  for  fragment  6A-2  west  since  the 
liner  was  foimd  at  a  greater  distance  than  the  overwrap. 
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The  velocity  of  the  top  fragment  of 
vessel  PC  was  measured  using  both  contact 
times  through  a  wire  grid  and  photo  tracking 
at  the  exit  to  the  wire  framework.  The  data 
from  both  types  of  measurements  is  plotted  in 
Figure  6.3.  The  photo  data  shows  a  lower 
velocity  than  the  wire  data  and  shows  a  much 
more  rapid  deceleration  due  to  drag  than 
would  be  computed  for  the  intact  vessel  half. 

The  wire  data  shows  no  deceleration  within 
the  wire  grid  framework,  however  there  might 
have  been  an  effect  on  the  maximum  velocity 
attained,  during  the  acceleration  phase. 

The  contact  wire  data  for  vessel  PC 
appeared  to  be  very  erratic.  Negative  voltages 
were  measured  where  none  were  applied  (a 
possible  capacitance  effect  as  the  shock  passed 
through  the  framework)  and  these  were  seen  prior  to  fragment  contact.  However,  by  taking  the 
initial  sharp  positive  break,  as  shown  in  Figures  6.4a  and  6.4b,  consistent  data  was  achieved. 

The  west  fragment  from  vessel  6A-2  showed  an  initial  velocity,  for  three  points,  of  1252 
fps.  The  velocity  then  dropped  suddenly  to  889  fps.  Figure  6.5,  and  decayed  rapidly  to  705  fys 
at  .090  seconds.  The  initial  three  points  tracked  may  be  the  correct  value.  The  velocity  is  at 
least  889  fys  and  may  be  as  high  as  1252  ^s.  The  east  fragment  carried  additional  weight  of 
part  of  the  vessel  support  and  had  a  lower  velocity. 

6.1.5  Endcap  velocities  of  vessels  5-3  and  5-4 

The  velocities  for  the  west  (end  cap)  fragments  of  vessel  bursts  5-3  and  5-4  were 
measured  with  contact  wire  instrumentation  that  was  itself  being  pushed  by  the  blast  and  must 
be  considered  minimum  velocities.  The  velocities  may  actually  be  fairly  accurate  based  on 
repeatability,  differences  in  acceleration  capability  of  wire  frame  versus  end  cap  and  direction  of 
initial  exhaust  during  fragment  acceleration.  Planned  three  dimensional  hydrocode  modeling  will 
include  this  case  as  a  check.  The  displacement  versus  time  data  for  both  of  these  fragments  is 
plotted  in  Figme  6.6. 


Figure  6.3,  Contact  Wire  and  Photo- 
Based  Velocity  Data  for  Vessel  PC 
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Figure  6.4a  Contact  Times  for  Vessel  PC,  Wires  1  through  3 
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Figure  6.4b  Contact  Times  for  Vessel  PC,  Wires  4  through  6 
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Figure  6.5,  Initial  Distance  versus  Figure  6.6,  Contact  Wire  Distances  Versus  Time 

Time  for  Vessel  6A-2  West  Fragment  for  West  Fragment  of  Vessel  5-3  and  5-4 

6.2  Calculation  of  initial  velocities/program  BURST 

As  fragment  velocities  were  measured  during  the  Bmst  Test  Study,  improvements  were 
made  to  a  computer  program  for  estimating  the  velocities. 

6.2.1  Background 

Taylor  and  Price^  in  1971  presented  a  technique  for  calcxilating  the  velocity  of  fragments 
resulting  from  a  bursting  gas  reservoir.  The  technique,  based  on  an  original  work  by 
Grodzovskii,  et  al^*,  removed  two  important  restrictions:  acceleration  with  small  displacements 
and  velocity  not  too  large.  The  technique  analyzed  the  displacement  of  two  fragments  from  a 
bursting  sphere  and  derived  a  differential  equation  for  the  motion.  Specific  cases,  to  compare 
to  empirical  data,  were  solved  using  numerical  integration  on  a  digital  computer. 

In  1972,  R.  A.  Lorenz,  in  an  appendix  to  Pittman’s’  report,  presented  computer  program 
FRAG  for  calculating  the  velocities  of  Wsting  spheres  and  cylinders.  The  work  was  part  of  test 
planning  for  those  pneumatic  bursts.  In  1975  Baker^  included  in  his  workbook  SPHER  and 
CYLIN,  computer  codes  for  calculating  the  initial  velocity  of  two  fragments  from  either  spherical 
or  cylindrical  vessels  respectively. 

Program  BURST  is  based  on  a  Pascal  version  of  the  Taylor  and  Price  code  as  written  and 
modified  by  Baicher^’.  Ambrosio'*®  has  rewritten  the  program,  made  a  number  of  improvements 
to  the  program  from  1990  through  1994,  and  documented  the  velocity  model  in  1993.  These 
improvements  removed  the  original  assumptions  of  perfect  gas  law  and  fragments  propelled  into 
a  vacuum.  The  improvements  were  made  as  a  cooperative  effort  on  the  Burst  Test  Study  to 
better  predict  the  fragment  velocities  that  were  being  measured.  The  latter  program  is  called 
BURST  and  is  the  program  that  will  be  discussed  herein. 
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For  additional  background  material  the  reader  is  referred  to  other  authors  whose  work  will 
not  be  covered  in  this  report.  Moore'"  presents  equations,  based  on  experiments  with  explosives 
by  Gumey^^  and  Steme^^  for  estimating  initial  fragment  velocity  by  assuming  that  60%  of  the 
initial  energy  is  converted  into  kinetic  energy.  Baum'*^  presents  data  from  his  own  research  and 
that  of  others  for  examining  the  accuracy  of  the  Moore  equations.  Bessey  and  Kulesz'*^  discuss 
an  analytical  model  and  based  on  modeling  results  show  that  fragment  velocities  increase  with 
the  number  of  fragments  until  ten  is  reached.  They  also  provide  curves  to  show  the  effect  of 
pressurant  to  fragment  mass  ratio  for  spheres  and  cylinders.  Baker^  presents  two  computer 
programs,  SPHER  for  spheres  and  CYLIN  for  cylinders,  for  estimating  initial  fragment  velocities. 
Wiedermann*  shows  through  analysis  that  fragment  velocities  are  significantly  reduced  by  real 
gas  effects.  Held  and  Jager^  present  equations  for  estimating  the  initial  velocity  of  fragments 
based  on  gas  discharge  time  constants.  Brown**®  and  Baker**’  have  both  written  summaries  of  the 
then  current  (1984)  state  of  the  art  in  fragmentation, 

6.2.2  Program  BURST  Characteristics*** 

The  Burst  program  assumes  a  variable  length  cylinder  capped  at  each  end  with 
hemispherical  domes.  The  latest  version  of  the  program,  9.0,  contains  four  rupture  geometry 
models: 

•  Multi-fragment  burst:  The  cylindrical  portion  ruptures  along  longitudinal  lines  and 
fragments  are  propelled  radially.  The  domes  are  propelled  longitudinally. 

•  Fractional  split  burst:  The  cylinder  ruptures  across  the  diameter  at  a  specified 
location.  The  cylindrical  portion  remains  intact  and  attached  to  its  respective 
dome. 

•  Dual  end  cap  split:  Each  dome  separates  from  the  cylinder  and  is  propelled 
longitudinally.  The  cylindrical  section  remains  intact. 

•  The  diial  end  cap  split  model  will  calculate  the  initial  velocity  of  a  sphere  by 
inserting  a  zero  len^  cylindrical  section. 

The  following  gas  expansion  and  fragment  motion  modeling  considerations/assumptions 
are  incorporated  in  the  program: 

•  Real  gas  properties  are  used  (air,  N2,  H2,  and  C02  are  available) 

•  Expansion  of  gas  is  assumed  to  be  adiabatic  (isentropic) 

•  Gas  temperature,  speed  of  sound  and  gas  compressibility  computed  from  gas 
expansion 

•  Gas  expansion  is  computed  independently  for  each  fragment 

•  Amount  of  gas  escaping  into  gaps  between  fragments  proportioned  according  to 
the  square  of  the  critical  velocities 
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•  Flow  coefficients  (which  detennine  transition  from  choked  to  unrestricted  flow) 
are  fimctions  of  ratio  of  gap  to  outer  radius 

•  Atmospheric  drag  is  accounted  for  in  acceleration  computation 

•  Fragment  tumbling  is  not  considered 

•  Velocity  of  fragments  is  taken  into  account  in  computing  pressure  exerted  on  the 
fragment 

•  Computations  for  a  section  are  stopped  when  critical  pressure  drops  below 
atmospheric  pressure 

6.3  Measured  velocities  compared  to  BURST 

A  comparison  of  measured  initial  velocities  and  initial  velocities  as  computed  with 
program  BURST  are  provided  in  Table  6.2.  Also  included  are  measurements  of  displacement 
and  time  to  maximum  (initial)  velocity.  Most  of  the  initial  displacement  and  time  data  were  not 
available  because  reference  locations  in  use  on  the  fragments  were  not  visible  at  the  start  of 
acceleration  due  to  the  dust  and  condensation  cloud.  Displacement  and  time  of  maximum 
velocity  firom  the  BURST  model,  however,  were  also  calculated  and  are  included  in  the  table. 

Various  investigators  have  estimated  the  percentage  of  total  energy  that  appears  as  kinetic 
energy  in  the  velocity  and  mass  of  the  fi'agments.  This  value  has  been  computed  for  each  burst 
and  is  provided  as  Table  6.3.  The  table  includes  an  "effective  fragment  velocity".  This  velocity, 
when  used  with  the  total  mass  of  the  fi-agments,  produces  the  kinetic  energy,  i.  e.  w^,v^jjf  = 
ZWjV;^.  A  value  of  500  ft/sec  is  used  for  the  sidewall  fragments  from  vessel  6A-4.  This  is  not 
an  average  velocity  of  the  sidewall  firagments  but  is  rather  the  maximrun  velocity.  It  is  used  for 
an  average  because  it  is  the  only  fragment  that  was  measured  on  a  track  almost  normal  to  the 
camera  field  of  view.  The  actual  average  may  be  less  but  is  considered  to  be  unknown  because 
of  camera  angles.  The  fragment  tracked  at  500  fit/sec  is  a  south  side  panel  that  was  launched 
horizontally  and  struck  a  dirt  moxmd  just  outside  the  arena.  Four  other  fi’agments  were  tracked. 
These  originated  fiom  the  top  or  north  side  of  the  vessel.  Their  velocities  and  angles  from  the 
horizon  are  228  @  30,  314  @  33,  314  @  81  and  338  fit/sec  @  69®.  For  reasons  yet  unknown 
these  were  not  measured  with  the  camera  having  the  best  angle.  In  fact  two  cameras  are  probably 
needed  to  determine  the  velocity  if  the  track  was  not  normal  to  the  field  of  view.  A  total  of  7 
cameras  were  in  use  as  shown  in  Figure  3.38. 
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Table  6.2 

Burst  Comparisons  with  Fragment  Test  Data 


'Computed  by  Chrostowski 

^ei^t  includes  part  of  support  structure 

^Program  revised  to  compute  dual  endcap  but  data  not  available  at  time  of  writing 
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Table  6.3 

Kinetic  Energy  of  Fragments  as  a  Percent 
of  Isentropic  Expansion  Energy 


Vessel 

Weight,  lbs 

Effective 
Velocity,  fps 

KE/IE,  % 

P-1 

5800 

239 

10.5 

P-2 

5800 

306 

11.6 

1-1 

5525 

148 

9.0 

1-2 

2900 

247 

10.7 

1-3 

5825 

315 

10.7 

1-4 

5250 

379 

10.4 

2-1 

5025 

251 

9.4 

2-2 

5300 

256 

10.4 

2-3 

5400 

265 

11.2 

PC 

44 

978 

21.0 

5-1 

4925 

261 

9.6 

5-2 

5100 

253 

9.3 

5-3 

4975 

234 

7.8 

5-4 

4975 

233 

7.7 

6A-1 

1600* 

445 

10.0 

6A-2 

112 

770 

33.2 

6A-3 

6100 

238 

10.8 

6A-4 

5048 

474 

33.3 

"Excluding  center  section  that  remained  in  place. 


6.4  Further  Velocity  Correlation 

Although  improvements  were  made  to  the  BURST  program  for  computing  fragment 
velocities,  there  are  still  some  limitations  to  the  code,  namely  there  are  failure  cases  not 
addressed  by  the  program.  In  his  1988  paper,  Baum'*  enumerated  eleven  failure  cases  for  either 
spheres  or  cylinders.  Seven  failure  cases  for  pressurized  gas  vessels  were  treated  in  the  paper. 
Five  equations  covered  the  seven  cases.  These  equations  determine  the  failure  type  in  a 
computer  program  UMITV,  written  for  a  hazard  assessment  workbook*®.  Three  of  the  cases 
are  appropriate  for  correlation  of  Burst  Test  Study  data  and  are  included  herein. 
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Baum’s  equations  provide  guidelines  on  fragment  velocity  that  are  realistic  upper  limit 
values  derived  from  experimental  data.  The  data  were  either  correlated  by  parameters  indicated 
by  dimensional  analysis  or  for  a  limited  number  of  cases  as  a  function  of  the  total  energy 
available.  In  some  cases  a  simple  experimental  model  was  used.  Pertinent  graphs  from  Baum’s 
paper  are  included  in  this  section.  New  graphs,  showing  correlation  with  experimental  data 
gathered  for  this  report  are  also  included.  His  equations  used  in  LIMITV  are  provided  in  Table 
6.4. 


The  term  F,  equation  6.4,  is  a  dimensionless  initial  acceleration.  The  term  E,  equation 
6.5  is  the  maximum  expansion  work  available  from  the  pressurized  contents.  This  is  an 
additional  equation,  beyond  equation  4.1  to  describe  the  energy.  Note  that  equation  6.5  using 
only  the  first  term  from  6.6  is  the  same  as  equation  4.1.  The  second  term  of  equation  6.6  is  a 
negative  correction  term,  making  6.5  less  than  4.1. 

All  of  the  models  used  herein  assume  an  ideal  gas.  Section  6.4.2  will  show  that,  for 
certain  burst  geometries,  the  results  are  also  applicable  for  a  real  gas  but  are  more  conservative. 

Some  of  the  burst  tests  conducted  under  this  study  are  appropriate  to  compute  with  either 
the  BURST  model  or  the  LIMITV  model.  These  are  shown  in  Table  6.5  along  with  the  test 
data.  Typically  BURST  calculated  velocities  that  were  slightly  lower  than  test  data  although 
some  are  high.  Since  the  LIMITV  velocities  are  limiting  velocities  they  can  be  expected  to  be 
high. 


6.4.1  End  cap  fragment 


Figure  6.7  is  Baum’s**  Figure  1,  for  endcap  data,  with  an  average  velocity  for  the  endcaps 
of  Test  5.3  and  5.4  from  Table  6.5  added.  Also  included  is  the  velocity  computed  by  the 
BURST  model,  (5%  less  than  measured).  (See  Section  6.1.5:  the  test  velocities  must  be 
considered  to  be  minimum  velocities  although  the  results  appear  to  accurate.) 

f 


Figure  6.7,  Test  5.3  and  5.4  endc(g>s  added  to  Baum’s  endcap  data 
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Table  6.4 

Equations  used  in  LIMITV 


*  See  the  workbook^®  for  types  3  and  4 
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Table  6.5 

Comparison  of  Test  and  Model  Velocities 
(Rocket  and  end  cap  fragments) 


Vessel 

:  Pressure 
psig 

Frag. 

Wei^t, 

lbs 

1-1 

1475 

West 

2700 

East 

2825 

1-2 

3450 

West 

2850 

East 

3050 

P-2 

4700 

West 

2715 

East 

3085 

1-3 

5425 

West 

2825 

East 

3000 

1-4 

7125 

West 

2550 

East 

2700 

5-1 

3600 

West 

1425 

East 

3500 

5-2 

3600 

West 

1500 

East 

3600 

5-3 

3600 

West 

300 

East 

4675 

5-4 

3600 

West 

300 

East 

4675 

6A-3 

3300 

West 

3050 

East 

3050 

6A-4 

3500 

End 

Caps 

2@352 

Side 

Frags 

12@362 

Cy!.  Gas 

length,  teiiq> 

ft.  i  “F 


9.5  .93  65 


9.5  .93  109 


9.5  .93  124 


9.5  .93  75 


‘Rocket  model  unless  shown  otherwise 

^End  cap  model 

^Cylinder  disintegration  model 


^Average  of  both  cq)s 

*Max.  - 1  frag,  others  may  be  less 

*A  minimum,  to  be  checked  with  hydrocode 


6.4.2  Rocketing  Fragment 


Most  of  the  Burst  Test  Study  data  falls  into  this  category  and  data  comparisons  are  shown 
in  two  graphs. 

Baum’s  plot  for  rocketing  fragments  (his  Figure  2)  is  presented  as  Figure  6.8.  Baum  used 
the  vessel  length  in  the  F(L/R)^'^  expression.  In  Figure  6.9  the  Held  and  Jager^  data  is  plotted 
using  L  both  as  a  vessel  length  and  as  a  fragment  length  (some  unavailable  ancillary  data  was 
estimated).  This  has  no  change  on  the  ordinate  for  each  point  but  the  increase  in  abscissa  with 
a  longer  L  places  the  points  lower  with  respect  to  the  line.  For  this  plot  the  use  of  the  fragment 
length  was  a  slightly  under-conservative  estimate.  The  BURST  computed  velocity  for  each  test 
is  also  plotted  at  the  fragment  length  abscissas. 

For  the  Burst  Test  Study  data  shown  in  Figure  6.10,  L  was  taken  as  the  fragment  length. 
The  recommended  limit  appears  to  be  a  conservative  estimator  of  the  test  data  (filled  symbols), 
typically  more  accurate  than  the  BURST  computations  (open  symbols)  except  at  higher 
pressures.  The  sonic  velocity,  Uo,  in  Figure  6.10  was  computed  for  air  as  an  ideal  gas.  Sonic 
velocity  for  the  actual  test  gas,  nitrogen,  as  a  real  gas,  was  used  in  Figure  6.11  to  replot  the  data 
from  Figure  6.10.  The  sound  velocity  values  are  available  from  BURST  output.  In  the  new 
plot  the  recommended  limit  is  no  longer  a  good  conservative  estimator  but  the  BURST  computed 
values  are  still  typically  low.  The  percentage  difference  is  not  large  as  can  be  seen  in  Tables 
6.2  and  6.5. 


SYMBOL 

ai 

17^ 

msm 

SOURCE 1 

• 

teg 

narH 

HI 

• 

10.4 

jj 

U 

15.4 

-27.0 

1  o 

sy 

ni 

Ba 

H 

HERZOG  et  ol  (1981) 

in 

,13.5 

□□ 

100-294 

HELD8.JAGER{1982) 

Figure  6.8,  Rocketing  Fragments  (from  Ref.  4) 
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Figure  6,9,  Held  and  Jaget^  rocketing  fragments  using 
L  for  both  vessel  length  and  fragment  length 


Figure  6,10,  Burst  Test  Study  Rocketing 
Fragments  using  a^  for  ideal  gas 


Figure  6,11,  Burst  Test  Study  Rocketing 
Fragments  using  a^  for  real  gas 
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6.4.3  Disintegration  of  sphere  or  cylinder 

Disintegration  of  a  cylinder  is  modeled  by  BURST  providing  the  endcaps  remain  intact  and 
sidewall  fragments  are  full  length  strips  of  selected  quantity.  One  such  test  case  was  conducted, 
vessel  6A-4.  Disintegration  of  a  sphere  is  not  modeled  by  BURST.  If  BURST  is  needed  to 
predict  fragment  velocities  then  velocities  from  the  centerline  split  of  a  sphere  could  be  used. 
The  predicted  velocities  would  be  lower  than  expected  since  there  would  be  only  two  fragments 
and  the  fragment  velocity  increases  as  the  number  of  fragments  increases  until  approximately 
ten  is  reached.**^ 

Baum  presented  a  separate  graph  for  cylinders,  reproduced  here  as  Figure  6. 12,  and  for 
spheres,  reproduced  here  as  Figure  6.13.  Two  new  plots  are  presented  that  include  both 
cylinders  and  spheres.  (The  recommended  limit  velocity  equation  does  not  differentiate  between 
these  two  configurations.)  Figure  6.14  plots  the  data  for  the  function  F  versus  V/aj,  where  a,, 
is  the  sonic  velocity  using  the  ideal  gas  assumption. 
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Y 

Pc/Pe 
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o 

1.4 
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X 
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2 
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♦ 
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6 
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L 

1.4 

4.3 

2 

PITTMAN  (1972) 

Figure  6.12,  Disintegrating  cylinders  (from  ref. 4) 
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5 
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V 
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Figure  6.13,  Disintegrating  Spheres  (from  Ref.  4) 
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Figure  6. 14  includes  comparisons  to  test 
cases,  to  Baum’s  recommended  limit  and  to 
BURST  calculated  velocity.  BURST  can 
model  only  one  of  the  cases,  vessel  6A-4  -  a 
multi-fragment  burst.  However,  BURST 
computations  were  made  for  comparison 
purposes  with  other  test  configurations  in 
spite  of  the  fact  that  the  technique  does  not 
strictly  apply.  Baum’s  limit  velocity 
technique  does  not  apply  to  vessel  burst  6A- 
4,  since  the  end  caps  remained  intact,  and 
burst  6A-2  that  resulted  in  only  2  fragments. 
BURST  computed  velocities  higher  than 
measured  for  the  Esparza  and  Baker*®  data 
shown  in  Table  6.6.  BURST  computed 
velocities  less  than  measured  for  the  Pittman’ 
bursts,  some  of  which  occurred  at  pressures 
around  8000  psig.  One  point,  the  625  psig 
Pittman  cylinder  burst  was  slightly  above 
Baum’s  recommended  limit  for  the  V/a„  used. 


Figure  6.14,  Disintegrating  spheres  and  cylinders 
from  Burst  Test  Study  and  the  literature 

Burst  gas  temperatures  of  60°F  were  assumed 


for  both  the  Esparza  and  Baker  bursts  and  the  Pittman  bursts.  The  actual  sidewall  fragment 
velocity  was  also  above  the  recommended  limit  while  the  end  cap  velocities  were  considerably 
below.  The  greatest  error,  measured  versus  BURST,  occurs  at  the  highest  pressure  Pittman 
burst  for  which  BURST  is  36%  low  and  Baum’s  recommended  limit  is  26%  high. 


Table  6.6 

Esparza  &  Baker  Fragment  Data 


Test 

Yol 

illBIII 

Uncorrected 
Burst  Pres, 
psig 

Radius,  Ft‘ 

Vessel 
Weight,  lbs 

Measured 

Velocity, 

FPS 

4 

.0023 

300 

.082 

.044 

233.9 

5 

.0200 

200 

.168 

.258 

205.4 

8 

.0028 

750 

.088 

.180 

249.0 

17 

.0027 

401 

.086 

.066 

216.2 

‘Calculated  from  volume,  added  .001  ft  for  outer  radius 
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Baum  did  not  present  a  case  for  two  sphere  halves  propelled  by  a  pressurized  gas  burst  (only 
a  liquid,  not  covered  herein).  It  can  be  seen  by  Figure  6.16  that  BURST  did  a  good  job  of 
estimating  fragment  velocity  for  both  the  heavy  side  (includes  support  hardware)  and  the  light 
side  of  6A-2.  The  velocities  of  fragments  from  vessel  6A-4  were  also  closely  estimated. 

6.5  Fragment  Recovery  Distances 

A  fragment  from  a  pressure  vessel  that  bursts  at  high  pressure  can  be  massive  and  yet  travel 
a  considerable  distance,  presenting  a  significant  hazard.  The  distance  a  fragment  may  travel  is 
of  safety  interest  both  for  a  pressure  vessel  installation  and  for  test  planning  for  the  Burst  Test 
Study  (Section  3.2).  With  a  known  breakup  geometry  and  an  estimated  initial  fragment  velocity 
using  the  method  of  Section  6.2,  the  fragment  range  can  be  estimated  by  assigning  a  launch 
angle  and  using  a  ballistic  trajectory  computer  program,  a  number  of  which  are  available.  Baker 
published  computer  programs  FRISB  for  lifting  type  fragments  and  TRAJE  for  drag  type 
fragments  in  his  workbook^.  TRAJE  was  keyed  from  the  literature  and  modified  (it  was  written 
to  produce  a  family  of  solutions).  It  appeared  to  work  reasonably  well  using  a  test  case  of  a 
vertically  launched,  light  weight  fragment,  vessel  PC.  However  in  preparation  for  end  cap 
bursts  in  the  Burst  Test  Smdy  a  better  known  program  was  desired.  Program  TRAJ  was  offered 
for  use  by  the  Naval  Surface  Warfare  Center.  TRAJ  is  a  two  dimensional  drag  trajectory 
computer  program  that  will  run  on  any  IBM  compatible  machine  using  DOS  3.1  or  higher 
operating  system.  This  program  is  easier  to  use,  also  corrects  for  drag  and  will  calculate 
ricochets  using  an  empirical  ricochet  model.  The  ricochet  model  is  described  by  Chrostowski 
and  Collins^*.  The  program  was  documented  by  Montanaro^^. 

The  only  known  fragment  launch  angles  during  Burst  Test  Study  were  horizontal.  Thus 
only  the  ricochet  model  is  relevant  in  comparing  calculated  distances  with  actual  distances 
travelled  by  fragments.  Table  6.7  shows  distances  travelled  by  fragments  from  vessels  ruptured 
during  the  Burst  Test  Study  and  includes  comments  on  damage  inflicted  to  trees  by  fragments. 
All  of  these  fragments  from  the  study,  except  for  two  vessels,  were  laimched  horizontally  and 
reached  their  final  location  by  a  variety  of  means  such  as  cartwheeling,  bouncing,  skidding  or 
rolling.  The  type  of  ricochet  depended  on  such  things  as  ground  condition,  fragment  length  and 
mass.  A  half  vessel  fragment  that  buried  its  nose  in  soft  clay  such  that  the  aft  end  was  flipped 
over,  began  to  cartwheel.  A  horizontal  launch  was  confirmed  by  measuring  the  location  of  the 
first  ground  strike  position  where  it  could  be  identified  and  comparing  that  distance  to  a 
computed  distance  using  the  initial  velocity  and  the  acceleration  due  to  gravity.  The  final 
distance  where  the  fragment  was  found  was  therefore  achieved  by  bouncing  for  most  fragments. 
For  some  tests  a  video  camera  recorded  the  action  of  the  fragments.  The  west  fragments  of 
vessel  bursts  P-1  and  P-2  were  seen  to  cartwheel  many  times  and  the  ricochets  are  recorded  in 
Tables  6.8  and  6.9.  The  west  fragment  of  vessel  1-4  was  observed  to  bounce  only  three  times 
in  reaching  the  woods  approximately  800  feet  away.  The  east  fragment  of  vessel  1-4  spun  in 
the  horizontal  plane  as  it  bounced.  This  was  typical  of  the  east  fragments.  The  pressurant 
tubing,  attached  to  that  end  was  pulled  in  two  by  &e  motion  of  the  fragment  and  the  restraining 
action  of  an  anchor  block  approximately  80  feet  from  the  arena  center.  Since  there  was  no  plan 
to  improve  a  ricochet  model  the  ground  strike  depressions  were  not  filled  so  that  each  fragment 
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would  have  a  unique  and  identifiable  pattern.  The  east  half  fragment  typically  went  a  lesser 
distance  because  of  the  restraining  force  necessary  to  break  the  supply  tubing.  A  half  vessel 
striking  frozen  ground  would  skid  and/or  bounce.  It  is  likely  that  some  of  the  fragments  would 
have  travelled  a  greater  distance  had  they  been  launched  at  an  upwards  angle,  especially  a  45° 
angle  for  maximum  range. 

TRAJ  provided  a  reasonable  and  conservative  range  estimate  through  the  ricochet  model  for 
the  horizontally  launched  fragments,  even  though  the  correlation  with  number  of  ricochets  is 
poor.  The  ricochet  depends  in  part  on  the  "soil  constant".  A  soil  constant  of  0.5  is  reasonable 
for  the  Dahlgren  test  area.  Table  6. 10  provides  ricochets  and  distances  for  the  west  fragment 
of  vessel  P-2  for  various  values  of  the  soil  constant  as  calculated  by  the  program. 

In  vessel  test  PC  one  fragment  was  launched  vertically  up  and  the  other  down.  After  vessel 
6A-4  was  prepared  for  testing  the  vessel  grooves  were  measured  with  respect  to  the  vessel 
support  to  determine  if  the  top  fragments  would  be  launched  straight  up.  A  small  angle  of  1.7° 
toward  the  north  was  discovered.  This  angle  and  an  initial  velocity  of  500  ft/sec  would  not 
produce  the  distance  travelled  by  some  of  the  top  fragments,  892’  to  1640’,  using  a  ballistic 
trajectory.  An  initial  velocity  of  8500  ft/sec  would  be  required  (program  TRAJ).  Therefore  the 
trajectory  was  not  ballistic,  the  shape  of  each  quarter-circular  section  affected  its  travel  by 
aerodynamics;  lift  and  variation  of  angle  of  attack.  Recovery  distances  for  vessel  6A-4  are 
shown  in  Table  6.7. 
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Table  6.7 

Fr^ment  Distances 
from  Burst  Test  Study 
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7.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  test  program  fulfilled  its  purpose  of  augmenting  the  available  pneumatically  pressurized 
vessel  burst  test  data.  Both  overpressure  and  fragment  data  were  collected  on  a  variety  of 
vessels  and  burst  scenarios  and  both  types  of  data  were  compared  to  existing  models.  An  existing 
model  for  fragment  velocity  was  rewritten  and  improved  as  a  cooperative  effort  with  Acta  Inc. 
Damage/injury  data  developed  during  explosives  testing  can  now  find  greater  application  to 
pressure  vessel  bursts. 

For  vessels  that  exhaust  rapidly  on  bursting  it  was  found  that  the  overpressure  and  impulse 
could  be  approximated  using  Baker’s  workbook  for  spheres  but  not  for  cylinders.  With  regard 
to  his  workbook,  slow  exhaust  vessels  produced  less  overpressure  and  cylinder  corrections  are 
much  less  than  shown.  The  Pressure  Vessel  Hazard  Assessment  Workbook  provides  models  for 
computing  overpressure  and  impulse  for  Burst  Test  Study  tested  vessels.  (For  lack  of  available 
information,  reflected  impulses  were  assumed  to  have  the  same  reflection  factors  as  their 
attendant  overpressures.)  If  there  is  a  very  close  reflecting  surface  then  reflected  overpressures 
and  impulses  may  be  produced  that  exceed  computations  by  either  workbook. 

It  is  felt  that  height  of  burst  (HOB)  testing  might  have  produced  more  meaningful  results 
if  rapid  exhaust  vessel  bursts  were  conducted  and  at  more  than  three  heights.  The  use  of 
reflection  factors  based  on  incidence  angle  and  incident  overpressure  from  testing  with  high 
explosives  is  recommended  until  data  are  available  from  additional  testing  or  hydrocode 
simulation.  (See  Section  5.6.) 

Pressure  vessel  burst  overpressures  cannot  be  computed  using  TNT  equivalence  except 
in  the  far  field  (Section  5.8).  In  the  near  field  TNT  equivalence  result  is  higher  than  actual. 
In  the  far  field  measured  overpressures  approached  TNT  results  for  most  vessels.  Even  in  the 
far  field  the  two  spheres  produced  overpressures  that  exceeded  TNT  equivalence  results. 

The  initial  velocities  of  fragments  were  computed  with  BURST  software  that  was  improved 
during  the  course  of  the  effort.  Close  results  were  obtained  with  computations  being  typically 
slightly  lower  than  actual  velocity.  Baum’s  limit  models  were  found  to  provide  useful  limiting 
velocities  for  the  cases  not  appropriate  for  BURST  (Section  6.4). 

Additional  testing  and/or  hydrocode  solutions  to  bursting  vessels  would  be  very  useful. 
Efforts  should  be  directed  toward  obtaining  additional  data  in  the  areas  of  lower  pressures,  larger 
vessels,  reflection  factors  for  both  overpressures  and  impulses,  cylinder  corrections  to 
overpressure  and  impulse  and  bursts  with  other  fluids  such  as  argon,  helium  or  a  propellant. 
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APPENDIX  A 


Drawings  and  Sketches  of  Test  Vessels 


A-2 


APPENDIX  B 


Pressurization  System  Diagram 


APPENDIX  C 


NSWC  Data  Sheets 


High  Explosives 

C-l  . PHE-2,  65#  C-4,  SHOT  3201  &  PHE-1,  33.3#  COMP  B,  SHOT  3249 

C.2  . HEl-1  .66#,  SHOT  3250  &  HEl-2,  9.9#,  SHOT  3256 

C-3  . HE2-3,  4.5#,  SHOT  3399  &  HE2-1,  4.5#,  SHOT  3401 

C-4  . HEl-3,  4.5#,  SHOT  3526  &  HE2-2,  4.5#,  SHOT  3404 

Vessels 

C-5  . PC,  3975  psi  @  3.5’  HOB,  SHOT  3573  &  1-1,  1475  psi  @  3.5’  HOB,  SHOT  3251 

C-6  .  P-1,  3250  psi  @  3.5’  HOB,  SHOT  3248  &  P-2,  4700  psi  @  3.5’  HOB,  SHOT  3245 

C-7  .  1-1,  1475  psi  @  3.5’  HOB,  SHOT  3251 

C-8  .  1-2,  3450  psi  @  3.5’  HOB,  SHOT  3252  &  1-3,  5425  psi  @  3.5’  HOB,  SHOT  3253 

C-9  .  1-4,  7125  psi  @  3.5’  HOB,  SHOT  3402  &  2-1,  3450  psi  @  3.5’  HOB,  SHOT  3403 

C-10  .  2-2,  3450  psi  @  8.7’  HOB,  SHOT  3400  &  2-3,  3475  psi  @  14.0’  HOB,  SHOT  3398 

C-11  .  5-1  (5),  3600  psi  @  3.5’  HOB,  SHOT  3574  &  5-2  (6)  ,  3600  psi  @  3.5’  HOB,  SHOT  3575 

C-12  .  5-3  (7),  3600  psi  @  3.5’  HOB,  SHOT  3577  &  5-4  (8),  3600  psi  @  3.5’  HOB,  SHOT  3576 

C-13 .  6A-1  (13),  3280  psi  @  3.5’  HOB,  SHOT  4053  &  6A-2 , 4000  psi  @  3.5’  HOB,  SHOT  4054 

C-14 .  6A-3  (4),  3300  psi  @  3.5’  HOB,  SHOT  4055  &  6A-4  (10),  3500  psi  @  3.5’  HOB,  SHOT  4056 


HOB  GAGES 


HEl-3,  4.5#,  SHOT  3526  HE2-2,  4.5#,  SHOT  3404 


1-1,  1475  psi  @  3.5’  HOB,  SHOT  3251 


Pmax  (psi)  I  MAX  IMPULSE 
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B41 

105® 

B47 

60 

u 

B49 

30 

®  ( 

B68 

0 

®  i 

B71 

150® 

B72’‘ 

135® 

B73 

120® 

B75 

105® 

B82 

45® 

B83 

30® 

B84 

15® 

75*  @  85' 


Pmax  (Dsi 


122.8 


NO  DATA 


10.8 


NO  DATA 


5.2 

2.8 


MAXIMPULSE(psi-sec 


1.569 


NO  DATA 


0.033 


NO  DATA 


_ 0.032 _ 

0.015 


0.030 

0.041 

0.054 

0.035 


0.032 


0.014 


0.025 


0.018 


0.010 

0.005 


22' 

1.0 

0.006 

22' 

BAD 

BAD 

50' 

1.2 

50' 

1.0 

50' 

1.2 

50' 

1.8 

50' 

1.2 

50' 

1.1 

0.008 


0.007 


07 

06 


>  130’  2" 

BAD 

BAD 

>  130'  4" 

1.5 

0.013 

B130* 

75® 

192'  5" 

0.4 

0.004 

HOB  GAGES 


C-8 


HOB,  SHOT  3400  2-3,  3475  psi  @  14.0’  HOB,  SHOT  3398 


5-1  (5),  3600  psi  @  3.5’  HOB,  SHOT  3574  5-2  (6)  ,  3600  psi  @  3.5’  HOB,  SHOT  3575 
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MAX  IMPULSE 


6A-1  (13),  3280  pst  @  3.5’  HOB,  SHOT  4053  6A-2 , 4000  psi  @  3.5’  HOB,  SHOT  4054 
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APPENDIX  D 

Overpressure/Impulse  Coefficients 


Index 

Pressure  Coefficients  for  LSC  Detonations . D-1 

Symmetric  Pressure  Coefficients  for  Vessel  Bmrst . D-2 

Symmetric  Impulse  Coefficients  for  Vessel  Burst . D-3 

Asymmetric  Pressure  Coefficients  for  Vessel  Burst . D-4 

Asymmetric  Impulse  Coefficients  for  Vessel  Burst  . D-5 


Table  D-1 

Pressure  CoefGcients  for  LSC  Detonations 


lilC  Charge 
N+ 

^ains 
per  ft 

•  Coefficient 

Error,  std 
dev.,  %  of 
reading 

: 

lA-2,  lA-3 
combined 

10 

2.57403 

-1.15727 

16.9 

lA-5,  lA-6 
combined 

15 

3.22303 

-1.29661 

19.3 

lA-1,  lA-4,  lA- 
7  combined 

25 

3.38636 

-1.24933 

19.9 

lB-2,  lB-3 
combined 

50 

3.11826 

-1.01572 

20.5 

lB-1,  lB-4 
combined 

200 

4.29318 

-1.22225 

25.9 

2-1 

25 

3.31845 

-1.26838 

30.0 

2-2 

25 

not  reduced,  similar  to  2-1 
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Table  D-2 

Symmetric  Pressure  Coefficients  for  Vessel  Burst 


Vessel 

CoefScient 

Error,  Std. 

Burst 

dev,,  %  reading 

P-2 

5.05345 

-.98558 

29.6 

P-1 

5.00237 

-1.05581 

26.5 

1-1 

3.71934 

-.90354 

46.8 

1-2 

4.44026 

-.83577 

24.0 

1-3 

5.14076 

-.98428 

26.1 

1-4 

4.83746 

-.86672 

30.0 

2-1 

4.39180 

-.88109 

23.9 

2-2 

4.92157 

-.96390 

21.1 

2-3 

4.75478 

-.86476 

19.2 

PC 

5.96755 

-1.45037 

16.1 

5-1 

4.20666 

-.86358 

25.4 

5-2 

4.43625 

-.93645 

29.7 

5-3 

4.32214 

-.93868 

47.4 

5-4 

4.48667 

-1.01354 

44.9 

6A-1 

6.41996 

-1.32488 

30.3 

6A-2 

5.61240 

-1.33255 

23.7 

6A-3 

4.51066 

-.81016 

38.0 

6A-4' 

7.31726 

-1.45505 

33.2 

6A-42 

6.47962 

-1.18289 

38.2 

'For  radial  distance 

^For  distances  normal  to  vessel  axis 
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Table  D-3 

Symmetric  Impulse  Coefficients  for  Vessel  Burst 


Vessel 

Burst 

Coefficient 

error,  std. 
dev.,  % 
of  reading 

B, 

P-2 

5.44016 

-.66163 

42.0 

P-1 

5.82760 

-.82635 

45.9 

1-1 

5.07290 

-.79185 

55.2 

1-2 

5.32406 

-.64697 

49.8 

1-3 

5.42485 

-.62706 

42.0 

1-4 

4.95254 

-.45454 

33.9 

2-1 

5.41929 

-.72406 

37.9 

2-2 

6.84018 

-1.12561 

37.4 

2-3 

7.72843 

-1.33472 

49.0 

PC 

4.96020 

-1.00060 

13.3 

5-1 

4.86168 

-.58916 

53.3 

5-2 

5.40752 

-.75853 

56.4 

5-3 

5.04190 

-.74506 

144.7 

5-4 

5.24812 

-.82165 

135.4 

6A-1 

6.33977 

-.82148 

24.8 

6A-2 

4.60186 

-.87193 

27.6 

6A-3 

4.54099 

-.40119 

42.5 

6A-4‘ 

7.30874 

-1.13451 

48.8 

6A-4^ 

6.79787 

-.93861 

39.7 

*For  radial  distances 

^For  distances  normal  to  vessel  axis 
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Table  D-4 

Asymmetric  Pressure  Coefficients  for  Vessel  Burst 


Vessel 

Coefficient 

error,  std. 

jtJurst 

B, 

Bs 

■ 

dev.,  % 
of  reading 

P-2 

5.80369 

-1.08903 

-.00932 

— 

80 

15.1 

P-1 

5.24688 

-1.04286 

-.00700 

— 

80 

17.9 

1-1 

3.90297 

-.91138 

-.00390 

— 

80 

43.4 

1-2 

4.67180 

-.81933 

-.00680 

— 

90 

11.9 

1-3 

5.45171 

-.98185 

-.00764 

— 

80 

14.3 

1-4 

5.22369 

-.90703 

-.00662 

— 

95 

22.4 

2-1 

4.63716 

-.87346 

-.00629 

— 

95 

13.0 

2-2 

.12293 

2.12668 

-.00425 

-.46848 

85 

15.2 

2-3 

6.98762 

-2.06836 

-.00424 

.17331 

80 

15.9 

PC 

5.90204 

-1.44416 

.00092 

— 

100 

15.7 

5-1 

4.54765 

-.88896 

-.00542 

— 

75 

16.9 

5-2 

4.73660 

-.91404 

-.00798 

— 

80 

16.6 

5-3 

5.00056 

-.99360 

-.00960 

— 

55 

15.9 

5-4 

5.13864 

-1.06984 

-.00862 

— 

50 

15.8 

6A-1 

6.60403 

-1.29551 

-.00640 

— 

95 

22.6 

6A-2 

5.80860 

-1.33704 

-.00216 

— 

5 

20.8 

6A-3 

4.82324 

-.76885 

-.01018 

— 

95 

20.6 

6A-42 

7.49824 

-1.42274 

-.00660 

— 

80 

27.4 

6A-4^ 

6.73789 

-1.20199 

0 

0 

1 

1 

1 

19.8 

’  Best  fit  center  array  angle,  and  strongest  blast  angle  if  less  than  zero 
^  For  radial  distances 
^  For  distances  normal  to  vessel  axis 
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Table  D-5 

Asymmetric  Impulse  Coefficients  for  Vessel  Burst 


Vessel 

Coefficient 

error,  shi. 
dev.,  % 
of  reading 

c„ 

IBIIIII 

Ca 

C3 

P-2 

6.62655 

-.82523 

-.01474 

— 

80 

16.3 

P-1 

6.36686 

-.79778 

-.01544 

— 

80 

19.5 

1-1 

5.65255 

-.81658 

-.01230 

— 

80 

39.0 

1-2 

5.80677 

-.61270 

-.01418 

— 

90 

19.6 

1-3 

5.94240 

-.62301 

-.01271 

— 

80 

25.5 

1-4 

5.44440 

-.50587 

-.00843 

— 

95 

25.6 

2-1 

5.84232 

-.71089 

-.01085 

— 

95 

18.8 

2-2 

8.56167 

-1.88809 

-.01127 

.11325 

85 

23.5 

2-3 

20.38713 

-8.41012 

-.01482 

1.02299 

80 

19.4 

PC 

4.98115 

-1.00259 

-.00030 

— 

100 

13.1 

5-1 

5.52768 

-.63873 

-.01059 

— 

75 

27.8 

5-2 

5.91690 

-.72052 

-.01354 

— 

80 

26.5 

5-3 

6.51361 

-.86414 

-.02082 

— 

55 

44.3 

mm 

6.77268 

-.95328 

-.02015 

— 

50 

48.7 

6A-1 

6.45241 

-.80351 

-.00392 

— 

95 

22.2 

6A-2 

4.90008 

-.87874 

-.00329 

— 

5 

20.9 

6A-3 

4.91664 

-.35155 

-.01223 

— 

95 

29.6 

6A-42 

7.58337 

-1.08548 

-.01002 

— 

80 

34.8 

6A-4^ 

6.98007 

-.95196 

0 

0 

— 

36.8 

’  Best  fit  center  array  angle,  and  strongest  blast  angle  if  B2,  Table  D-4,  less  than  zero 
(based  on  pressure  measurements) 

^  For  radial  distances 
^  For  distances  normal  to  vessel  axis 
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APPENDIX  E 


OVERPRESSURE/IMPULSE  PLOTS 


Overpressure,  psig 


(file  =  P3201DTR.SPG) 


(file  =  P3201ATR.SPG) 


(file  =  B201DTR.SPG) 


Figure  E.1,  65  lbs  C-4  detonated  at  3.5'  HOB  above 
center  of  arena,  6-21-90  (PHE-2)  (94.0  lbs  TNT  equiv.) 


Overpressure,  psig 


(file  =  P3249DTR.SPG) 


a 

jT 


30  60  90  120  150 

Arena  angle,  deg. 


(file  =  P3249ATR.SPG) 


(file  =  B249DTR.SPG) 


Figure  E.2, 333  lbs  of  composition  B  detonated  at  3.5'  HOB 
above  center  of  arena,  6-26-90  (PHE-1)  (40.0  lbs  lOT  equiv.) 


180 


O  ,•  HOB,  dose  ground  pomls 
.  .  Coefs.,  osymmetric,  mox. 

_ ^ _  Coefs,  overage 

_  TNT  equiv,  vessel  or 

high  explosive 


10  30 

Distance,  ft 


1 00  200 


30  60  90  120  15 

Arena  angle,  deg. 


(file  =  P3250DTR.SPG) 


(file  =  P3250ATR.SPG) 


O  •  #  HOB.  dose  ground  points 

- Coefs.,  osymmetric,  mox. 

- ^ -  Coefs,  overoge 

_  TNT  equiv,  vessel  or 

high  explosive 


. 


10  30 

Distance,  ft 


1 00  200 


0  30  60  90  120  150 

Areno  angle,  deg. 

Distonces  per  oreno  layout,  on  the  ground 
not  corrected  for  offset 


(file  =  B250DTR.SPG) 


(file  =  I3250ATR.SPG) 


Figure  E.3,  0.66  lbs  of  pentolite  detonated  at  3.5’  HOB 
above  center  of  arena,  1-28-91  (HE  1-1)  (.89  lbs  TNT  equiv.) 


Overpressure,  psig 


(file  =  P3256DTR.SPG) 


(file  =  P3256ATR.SPG) 


(file  =  I3256DTR.SPG) 


Figure  E.4,  9.9  lbs  of  pentolite  detonated  at  3.5'  HOB 
above  center  of  arena,  1-28-91  (HE  1-2)  14.6  lbs  TNT  equiv.) 


Overpressure* 


(file  =  P3526DTR.SPG) 


(file  =  P3526ATR.SPG) 


(file  =  B526DTR.SPG) 
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Figure  E.5,  4^  lbs  of  pentolite  detonated  at  3.5'  HOB, 
above  center  of  arena,  7-20-91  (HE  1-3)  (6.6  lbs  TNT  equiv.) 


Overpressure^  psig 


(file  =  P3401DTR.SPG) 


(file  =  P3401ATR.SPG) 


(fiile  =  I3401DTR.SPG) 


Figure  E.6,  4.5  lbs  of  pentolite  detonated  at  3.5'  HOB 
above  center  of  arena,  7-18-91  (HE  2-1)  (6.6  lbs  TNT  equiv.) 


Overpressure,  pslg 


(file  =  P3404DTR.SPG) 


(file  =  P3404ATR.SPG) 


(file  =  I3404DTR.SPG) 


Figure  E.7,  4.5  lbs  of  pentolite  detonated  at  8.7'  HOB,  offset 
1’  away  from  center  of  arena,  7-20-91  (HE  2-2)  (6.6  lbs  TOT  equiv.) 


Overpressure,  psig 


(file  =  P3399DTR.SPG) 
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Figure  £.8, 4.5  lbs  of  pentolite  detonated  at  14'  HOB,  offset  1' 
away  from  center  of  arena,  7-16-91  (HE  2-3)  (6.6  lbs  TNT  equiv.) 


180 


Impulse,  psi-ms  Overpressure,  psig 


(ffle  =  P3248DTR.SPG) 


(file  =  P3248ATR5PG) 


Distances  per  areno  loyout,  on  Uie  ground, 
not  corrected  for  offset 


(ffle  =  B248DTR.SPG)  (file  =  I3248ATR.SPG) 


Figure  E.9,  2'  dia.  steel  vessel,  53  cu.  ft.  volume,  3250  psi,  split  at  near  center  at  3.5'  HOB, 
1'  offset  away  from  center  of  arena,  6-26-90  (P-1)  (TNT  equivs:  gas  =  31.7,  LSC  =  .07) 


Overpressure,  psig 
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Figure  E.10, 2'  dia.  steel  vessel,  53  cu.  ft.  volume,  4700  psi,  split  at  near  center  at  3.5'  HOB, 
1'  offset  away  from  center  of  arena,  6-25-90  (P-2)  (TNT  equivs:  gas  =  47.0,  LSC  =  0.0) 


Impulse,  psi-ms  Overpressure,  psig 


(ffle  =  P3251DTR.SPG) 


(file  =  P3251ATR.SPG) 


Distonces  per  arena  layout,  on  the  ground, 
not  corrected  for  offset 


(ffle  =  B251DTR.SPG) 


(ffle  =  I3251ATR.SPG) 


Figure  E.11, 2'  dia.  steel  vessel,  53  cu.  ft.  volume,  1475  psi,  split  at  center  at  3.5'  HOB, 
1 '  offset  away  from  center  of  arena,  1-25-91  (1-1)  (TNT  equivs:  gas  ~  13.5,  LSC  =  0.27) 


(file  =  P3252DTR.SPG) 


(file  =  P3252ATR.SPG) 


(file  =  I3252DTR.SPG) 


Figure  E.12, 2'  dia.  steel  vessel,  53  cu.  ft.  volume,  3450  psi,  spUt  at  center  at  3.5'  HOB, 
1'  offset  away  from  center  of  arena,  1-29-91  (1-2)  (TNT  equivs:  gas  =  33.8,  LSC  =  .07) 


Overpressure,  psig 


(ffle  =  P3253DTR.SPG) 


(ffle  =  P3253ATR.SPG) 


(file  =  B253DTR.SPG) 


Figure  E.13,  2'  dia.  steel  vessel,  53  cu.  ft.  volume,  5425  psi,  split  at  center  at  3.5'  HOB, 
1'  offset  away  from  center  of  arena,  1-26-91  (1-3)  (TNT  equivs:  gas  =  54.8,  LSC  =  .03) 


Overpressure,  psig 


(ffle  =  P3402DTR.SPG) 


(file  =  P3402ATR^PG) 
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Figure  E.14, 2'  dia.  steel  vessel,  53  cu.  ft.  volume,  7125  psi,  split  at  center  at  3.5'  HOB, 
1'  offset  away  from  center  of  arena,  7-19-91  (1-4)  (TNT  equivs:  gas  =  73.4,  LSC  =  .13) 


Overpressure*  psig 
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(ffle  =  P3403DTR.SPG) 
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Figure  E.15, 2'  dia.  steel  vessel,  53  cu.  ft.  volume,  3450  psi,  split  at  center  at  3.5'  HOB, 
V  offset  away  from  center  of  arena,  7-20-91  (2-1)  (TNT  equivs:  gas  =  33.8,  LSC  =  .03) 


Overpressure,  psig 


(ffle  =  P3400DTILSPG) 


(file  s  P3400ATILSPG) 


(file  =  B400DTR.SPG) 


Figure  E.16. 2'  dia.  steel  vessel,  53  cu.  ft.  volume,  3450  psi,  split  at  center  at  8.7'  HOB, 
1'  offset  away  fi*om  center  of  arena,  7-17-91  (2-2)  (TNT  equivs:  gas  =  33.8,  LSC  =  .03) 


Overpressure,  psig 


(ffle  =  P3398DTR.SPG) 


(file  =  P3398ATR.SPG) 


(file  =  I3398DTR.SPG) 


Figure  E.17,  2'  dia.  steel  vessel,  53  cu.  ft.  volume,  3475  psi,  split  at  center  at  14'  HOB, 
\>  offset  away  from  center  of  arena,  7-15-91  (2-3)  (TNT  equivs;  gas  =  34.1,  LSC  =  .03) 


Overpressure,  psig 


(file  =  P3573DTR.SPG) 


(file  =  P3573ATR.SPG) 


(file  =  I3573DTR.SPG) 


Figure  EI.18, 1S4 '  dia.  composite  vessel,  2.7  cu.  ft.  volume,  3975  psi,  split  at  horizontal  center 
at  3.5'  HOB,  V  offset  away  from  center  of  arena,  1-25-92  (PC)  (TNT  equivs:  gas  =  2.0,  LSC  =  . 
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Figure  E.19, 2'  dia.  steel  vessel,  53  cu.  ft.  volume,  3600  psi,  split  at  3/4  - 1/4  length 
at  3.5'  HOB  above  center  of  arena,  6-16-92  (5-1)  (TNT  equivs;  gas  =  35 J,  LSC  =  .C 


Overpressure,  psig 
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Figure  E2,0,  2*  dia.  steel  vessel,  53  cu.  ft.  volume,  3600  psi,  split  at  3/4  - 1/4  length 
at  3.5'  HOB  above  center  of  arena,  6-17-92  (5-2)  (TNT  equivs;  gas  =  35  J,  LSC  =  .O': 


(file  =  P3577DTR.SPG)  (file  =  P3577ATR.SPG) 


(file  =  B577DTR.SPG)  (file  =  I3577ATR.SPG) 


Figure  EJl,  2'  dia.  steel  vessel,  53  cu.  ft.  volume,  3600  psi,  split  at  full  length-endcap 
at  3.5'  HOB  above  center  of  arena,  6-19-92  (5-3)  (TNT  equivs:  gas  =  353,  LSC  =  .07) 


(file  =  P3576DTR.SPG) 


(file  =  P3576ATR.SPG) 


(file  =  I3576DTR.SPG)  (file  =  13576ATR.SPG) 


Figure  E^,  2*  dia.  steel  vessel,  53  cu.  ft.  volume,  3600  psi,  split  at  full  length-endcap 
at  3.5'  HOB  above  center  of  arena,  6-18-92  (5-4)  (TNT  equivs:  gas  =  353,  LSC  =  .07) 


Overpressure,  psig 
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(ffle  =  P4053DTR.SPG) 


(ffle  =  P4053ATR.SPG) 
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Figure  E^3,  3'  dia.  steel  vessel,  53  cu.  ft.  volume,  3280  psi, 
split  into  3  fragments  (2  endcaps  +  1  sidewall  piece)  at  3.5'  HOB 
above  center  of  arena,  11-2-92  (6A-1)  (TNT  equivs;  gas  =  32J1,  LSC  =  3.41) 


Overpressure,  psig 


(file  =  P4054DTR.SPG) 


(file  =  P4054ATR^PG) 


(file  =  I4054DTR.SPG) 


Figure  E24, 1.94'  dia.  composite  vessel,  2.7  cu.  ft.  volume,  4000  psi,  split  at  vertical  center 
at  3.5'  HOB  above  center  of  arena,  11-3-93  (6A-2)  (TNT  equivs:  gas  =  2.0,  LSC  =  .07) 
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Figure  E.25,  3'  dia.  steel  vessel,  53  cu.  ft.  volume,  3300  psi,  split  at  center 
at  3.5'  HOB  above  center  of  arena,  11-3-93  (6A-3)  (TNT  equivs:  gas  =  323,  LSC  =  39) 


(ffle  =  P4056DTR.SPG) 


(ffle  =  P4056ATR.SPG) 


(ffle  =  I4056DTR.SPG) 


(ffle  =  I4056ATR.SPG) 


Figure  2 '  dia.  steel  vessel,  53  cu.  ft.  volume,  3500  psi, 
split  into  14  fragments  (2  endcaps  +  12  sidewall  pieces)  at  3.5'  HOB 
above  center  of  arena,  11-4-93  (6A-4)  (TNT  equivs:  gas  =  34J,  LSC  =  2.18) 


Appendix  F 

Comments  on  Multi-fragment  ^linder  Overpressure  Measurements 


V 


4^ 


Figure  F.I,  Overpressure  versus 
Angle  for  Vessel  Burst  6A-4 


Figure  F.2,  Overpressure  versus 
Distance  for  Vessel  Burst  6A~4 


A  late  attempt  was  made  to  find  a  better  analysis  model  to  fit  the  overpressure  and 
impulse  data  firom  the  multi-fragment  cylinder,  vessel  6A-4.  Ten  multiple  linear  regressions  were 
performed^’  using  various  combinations  of  two  to  four  independent  variables.  The  two  best 
analyses,  using  two  variables,  regressed  the  log  of  overpressure  versus  the  log  of  radial  distance 
from  the  vessel  center  (as  the  first  independent  variable  for  both)  and  a  second  variable.  The 
second  independent  variable  was  for  1.)  the  axial  distance  from  the  vessel  center  (18.2%  error) 
and  for  2.)  the  angle  from  the  vessel  center  (21.8%  error).  Symmetry  about  a  line  perpendicular 
to  the  center  of  the  vessel  axis  was  assumed  for  all  analyses  and  Ae  three  points  at  five  foot 
distance  were  not  used.  Analysis  #2  uses  the  independent  variables  typically  used  for  the  other 
vessels.  Analysis  #1  is  peculiar  to  this  vessel.  The  use  of  all  three  independent  variables  reduced 
the  error  only  0.1%.  Error  used  is  the  standard  deviation  in  percent  of  reading  of  measured 
overpressure. 


Analysis  #1  is  shown  in  a  plot  of  arena  angle  versus  overpressure  in  Figure  F.l  and  in  a 
plot  of  radial  distance  versus  overpressure  in  Figure  F.2.  The  lines  of  constant  axial  distance 
versus  overpressure  are  curved  lines  in  Figure  F.l.  The  data  is  also  plotted  for  comparison  in 
Figure  F.l  using  both  filled  and  open  symbols  and  with  points  connected  using  a  spline.  Analysis 
#1  is  shown  by  the  three  solid  lines  in  Figure  F.2  labeled  O’,  25’  and  50’  axial  (distance).  The 
line  labeled  ''Max.-radial"  is  equivalent  to  the  line  labeled  "Coefs.,Asymmetric,  max.”  in  Figure 
5.26b.  The  line  labeled  "Max  -perpendicular"  in  Figure  F.2  is  the  same  line  labeled  "Max.  90® 
distance"  in  Figure  5.27.  This  line  was  found  by  regressing  overpressure  against  perpendicxilar 
distance  with  the  filled  symbol  points  in  Figure  F.l  deleted  but  with  the  three  points  at  5 
distance  retained. 


The  equations  for  overpressure  and  impulse  for  analysis  #1  are: 

log(OP)=6.478 11-1 .06481  *log(radial  distance)-0.02419*axial  distance 
log(Imp.)=5.82776-0.53925*log(radial  distance)-0.02700*axial  distance 


F-1 


